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SUMI-IARY 


In this report an analytical method is presented for predicting 
the lift coefficient, C^, the pitching moment coefficient, C^, and the 
drag coefficient, of light, twin-engine, propeller-driven airplanes. ' 

The method is applied to the Advanced Technology Light Twin-Engine 
(ATLIT) airplane. The calculated characteristics are correlated against 
full-scale wind tunnel data. 

In addition, results obtained with the "thick wing" lifting sur- 
face program of Reference 6 and the "thin wing" lifting surface pro-, 
gram of Reference 5 are used in the correlation. 

. The method described in this report predicts the drag and pitching 
moment fairly well. The lift prediction, however, is rather poor and 
needs improvement. In the case of lift prediction the lifting surface 
methods of References 5 and 6 show better agreement xjith the wind tunnel 
results of Reference 2. 
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propeller downwash, e (Section 
5.1) P 

2) Factors used in determining wing 
lift distribution (Section 4.2) 


Dimension 


m (in, ft) 


m (in, ft) 
m (in, ft) 

m (in, ft) 
m (in, ft) 

m (in, ft) 
m (in, ft) 
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LIST OF SYMBOLS (continued) 


Symbol 

Definition 

Dimension 


Factors used in determining wing 
lift distribution 


S 

Equivalent circular perimeter of 
the maximum frontal area of the 
body 

m (in, ft) 

S 

Equivalent circular perimeter 
of the base area 

m (in, ft) 


Drag coefficient 



Base drag coefficient 



1) Wing zero-lift drag coefficient 
of total wing (Section 4.9) 



2) Skin friction and pressure drag 
coefficient of the body (Section 
4.12) 


S. 

X 

Induced drag coefficient 


S 

0 

Zero-lift drag coefficient 



Zero-lift drag coefficient of the 
fuselage with fuselage interference 
accounted for 


0 

Zero-lift drag coefficient of the 
horizontal tail with interference 
effects accounted for 


^‘"D^^n(w) 

Net zero-lift drag coefficient of 
the nacelles in presence of the 
wing 


^S K 
0 

Net contribution of the vertical 
tail to the zero-lift drag coef- 
ficient with interference effects 
accounted for 
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LIST OF SYMBOLS (continued) 


Symbol 


Definition 


Dimension 


(Cjj Net zero-lift drag coefficient of 

o the wing-fuselage combination with 

interference effects accounted for 


C 


f 


C 


L 


C 


L 


f » 


^h(hf) 


Skin-friction coefficient of a flat 
plate 

Lift coefficient 

Lift coefficient of an isolated 
wing including the increments of 
lift due to the normal force of one 
propeller and due to the lift com- 
ponent of one propeller 

Lift coefficient of the horizontal 
tail with tail-fuselage interaction 
effects included j referred to wing 
area and free-stream dynamic pressure 


\(hf) 


Lift coefficient of the horizontal 
tail, referenced to the tail area 
(unless noted otherwise) , xjith tail- 
fuselage interaction effects, angle 
of attack, stabilizer deflection, — 
and tab deflection accounted for 


(C, ) ’ Net lift coefficient of the hori- 

■^(hf) — =1.0 zontal tail due to a^, i^^, and 

with tail-fuselage interaction 
effects included, referenced to 
the horizontal tail area and a. 
djnnamic pressure-ratio at the tail 
equal to one 



= 0 


Stabilizer effectiveness, 9C^/9i^, 
with tab fixed at zero setting 





Stabilizer effectiveness, 

Xh tab 0 "li 

the tab geared to the stabilizer 
to deflect in the ratio of 
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LIST OF SYMBOLS (continued) 


Symbol 

Definition 

Dimension 

\ 

Maximum lift coefficient 


max 



Cl 

w(f)+f (w) 

Lift coefficient of the wing inclu- 
ding mutual wing-fuselage interfer- 
ence effects. 


a 

Lift-curve slope 

deg rad ^ 

^h(f) 

a 

Liit-curve slope of the horizontal 
tail alone with fuselage effects 
on the tail accounted for 

deg ^ 

^^L^^h(h£) 

Lift-curve slope of the horizontal 
tail with interacting tail-fuselage- 
effects accounted for. 

deg ^ 

^Pol 

a 

Wing lift-curve slope according to 
Polhamus 

rad“^ 

\ 

tab 

Tab effectiveness, ^Cj^/96^^^ 

deg 

C 

m 

Pitching moment coefficient 


^“h(hf) 

Horizontal tail contribution to 
the pitching moment coefficient 
based oh 

ni(hf) 


\(h£) 

Contribution of C to pitching 

hi(hf) 

moment coefficient 


C 

mr 

Stabilizer effectiveness in pitch 
with the tab geared to the stabilizer 
to deflect in the ratio of 

deg~^ 

C 

m 

Zero-lift pitching moment coefficient 


o 

C 

m 

a 

Variation of the pitching moment 
coefficient with angle of attack 

deg ^ 
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LIST OF SYMBOLS (continued) 


Symbol 


^f(e)+n,(e) 

ct 

^w(D) 

a 

^^m \7(f)+f (w) 

Ott 


(c ) 

ra w 
a s 


K 

P 




c 

X 


c. 


c 



Definition 


Dimension 


Slope of the "free moment" deg 

coefficient of the body 

Sum of the "free moment" slopes deg 

of the fuselage and nacelles 


Slope of the pitching moment deg 

coefficient due to wing drag 

Slope of. the pitching moment deg 

coefficient about the leading 
edge of the wing mean aerodyn- 
amic chord due to the effective 
wing lift, including the effects 
of fuselage upwash on the wing and 
wing lift carryover onto the fuse- 
lage 

Average slope- of the wing pitching deg 
moment coefficient about the lead- 
ing edge of the wing mean aerodyn- 
amic chord in the non-linear lift 
range to stall 

Normal-force derivative of the pro- rad 
peller based on the propeller disk 
area 


Reference normal-force derivative rad 
of a propeller having a normal 
force factor, equal to 80.7 


Circumference of cross-section of m (in, ft) 
body at distance x from the nose 

Chord . m (in, ft) 

Mean aerodynamic chord m (in, ft)- 

Taper-ratio correction factor 



Steady-state crossflow drag coef-^ 
ficient 


Standard mean chord-, (c^ .+ c^)/2 


m (in, ft) 
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LIST OF SYMBOLS (continued) 


Symbol 


c . 

X 



Definition Dimension 

Mean aerodynamic chord of the m (in, ft) 

portion of the wing immersed in 
the propeller slipstream 

Chord at the intersection of the m (in, ft) 

vertical and horizontal tail 

Section lift coefficient 
"Additional" lift coefficient 

"Basic" lift coefficient 


c 


S, 


max 


Maximum section lift coefficient 



Section lift-curve slope at low 

deg 

rad 

Mach number (M _< 0.2) 



Section effectiveness of the tab. 

deg 

rad 


c 


m 


Section pitching moment coefficient 


c 

m 

o 


Section pitching moment coefficient 
at zero lift 


c 

n 


Chord of the wing at the centerline m (in, ft) 
of the nacelle 


c Wing chord m (in, ft) 

w 


d Fuselage width at the wing m (in, ft) 

dg Diameter of the equivalent circular m (in, ft) 

perimeter of the maximum frontal 
area of the body (fuselage or na- 
celle) 

Diameter' of the equivalent circular m (in, ft) 
perimeter of the base area 


-(dC /dC^) 
m L 


eg 


Static margin relative to the center 
of gravity as a fraction of the wing 
mean aerodynamic chord 
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LIST OF SYMBOLS (continued) 


Symbol 


d 

max 


e 


f 


1 

w 


(i ) 
wo 


i 


T 


K 


K’' 




Definition 


Fuselage width at the horizontal 
tail 

Maximum diameter of- ah equivalent . 
circular body 

Oswald span efficiency factor used 
in the induced drag equation 

1) Propeller inflow factor, ratio 
of the propeller normal force coef- 
ficient at power-on to power-off 

(T ’ = 0)- conditions- (section 5.1) 
c 

2) Lift distribution function 
(Section 4.2) 

Incidence of the horizontal tail 


Wing incidence at the root, angle 
between wing chord and X-body axis 

Incidence of the zero-lift line of 

the wing relative to the X-body 

axis, i - a 
w o 

w 

Incidence of the thrust axis rela- 
tive to X-body axis 

Correction factor for maximum lift 
due to power 

Correction factor for the lift effect- 
iveness of the tab at large tab de- 
flections 

Correlation parameter for additional 
wing lift due to power effects on 
the wing 

Span factor for inboard flaps or tabs 


Propeller drag factor 


^fCh) 


Ratio of the tail-lift carryover onto 
the fuselage to the tail alone 


Dimension 
m (in, ft) 

m (in, ft) 


deg 

deg 

deg , rad 
deg 
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LIST OF S'gfflOLS (contiimed) 


Symbol 


Definition 


Dimension 


^f(w) 




^ol 

^w(f) 


k 



Ratio of wing lift carryover on the 
fuselage to the wing alone 

Ratio of the lift on the tail in the 
presence of the fuselage to the tail 
alone 

Propeller normal force factor 


Error in Polhamus formula when com- 
pared with lifting surface method 

Ratio of the lift on the wing in the 
presence of the fuselage to the wing 
alone 

Correction factor to account for the 
effects of the wing planform on the 
increment of maximum lift coefficient 
due to tab position 

1) c /2 tt (Section 4.2) 

■V 

a 

2) d/b (Section 4.8) 

3) Equivalent sand roughness of a m (in, ft) 
surface (Section 4.12) 

Reduced mass factor 


Correction factor to account for 

c^ , /c other than 0.25 
tab 

k 2 Correction factor to account for 

tab deflection other than 60 degrees 

k^ Correction factor to account for 

chord extension due to tab deflection 


^f(h) 




Ratio of the lift carryover, due to 
stabilizer deflection, onto the fuse- 
lage to the lift of the stabilizer 
alone 

Ratio of the lift on the stabilizer, 
due to stabilizer deflection, in the 
presence of the fuselage to stabilizer 
alone 
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LIST OF SYMBOLS (continued) 


Symbol Definition 

LER Leading edge radius as a fraction 

of the chord 

I Reference length, for lifting sur- 

faces, equal to the mean aerodynamic 
chord of the surface, for bodies, 
equal to the length of the body 

Jl Length of body (fuselage or nacelle) 

B 




eff 



SL 


h 


t 


Z 

n 

Z 

ru 


Z 


2 


Distance, in the wing root-chord 
plane, from the tip vortex at the 
quarter chord of b to the 

eff 

quarter chord- of the horizontal- 
tail mean aerodjmamic chord 

Length of fuselage 


Distance from the nose of the fuse- 
lage to the quarter chord of the 
wing mean aerodynamic chord 

1) Distance from the wing trail- 
ing- edge to the centroid of the 
last aft Ax segment of the fuselage 
(Section 4.8) 

2) Distance from the center of. 
gravity to the quarter chord' of the 
horizontal tail mean aerodynamic 
chord (Section '5.1) 

Distance from the quarter chord of 
the mean aerodynamic chord of the 
immersed portion of the wing to the 
quarter chord of the horizontal tail 
mean aerodynamic chord 

Nacelle length 


Distance required for the complete 
rollup of the wing-tip -vortices 

Tail length in the wing-root chord 
plane from the root-chord trailing 
edge to the quarter chord of the 
horizontal tail mean aerodynamic 
chord 


Dimension 

m (in, . ft) 

m (in, ft) 
m (in, ft) 

m (in, ft) 
ra (in, ft) 

m (in, ft) 
m (in, ft) 

m. (in, ft) 

m (in, ft) 
m (in, ft) 

m (in, ft) 
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LIST OF SYMBOLS' (continued) 


Symbol 


Definition 


Dimension 


Distance from the leading edge of m (in, ft) 
the wing mean aerodynamic chord to 
the trailing edge of the wing root 
chord 

M Mach nximber 


M 

c 


M sina^. 


ef f 


Reynolds number 


n 


n 


1 


n 


2 


prop 


Number of propellers 

Number of junctures of the tail sur- 
face with the fuselage 

Number of corners in the juncture 
of the vertical tail with the hori- 
zontal tail (cruciform intersection 
would have four corners) 


Propeller 

Dynamic pressure ratio at the hori- 
zontal tail 

2 2 

Free stream dynamic pressure N/m (Ib/ft ) 


R Leading-edge suction parameter 


R 

P 

Propeller radius 

m 

(ft) 


\f 

Ratio of the wing-fuselage to the 
fuselage-alone zero-lift drag with 
the base drag omitted 




r 

Cross-section radius of equivalent 
circular body 

m 

(in, ft) 

S 

Area 

2 

m 

(in^, 

ft^) 


Maximum frontal area of the body 
(fuselage or nacelle) 

2 

m 

(in^. 

ft^) 

^f 

Planform area of the fuselage 

2 

m 

(in^, 

ft^) 
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LIST OF SYMBOLS (continued) 


Symbol 







wet 


T 

T ' 
c 


t/c 


(t/c) 


int 


V 


Definition 


Dimension 


Planform area of the fuselage for- 
ward of the quarter chord of the 
wing mean aerodynamic chord 

Area of the horizontal tail im- 
mersed in the propeller slipstreams 

Horizontal tail area in front of 
and including the tab 

Wing area immersed in the slip- 
stream of the propellers 

Area of the exposed panels of a 
lifting surface 


2 2 2 
m“ (in^, ft^) 


2 2 2 
m Cin , ft ) 


2 2 2 
m^ (in , ft ) 


2 2 2 
m (in , ft ) 


2 2 2 
(in , ft^) 


Cross-section area of an equivalent 
circular body. 

Disk area of propellers (total) 


2 2 2 

m^ (in, ftn 


m^ (in^, ft^) 


Cross-section area of an equiva- 
lent circular body at the fore- 
going station being considered 

Wetted area 


*7 9 2 

m^" (in , ft ) 


2 2 - 2 . 
m (xn , ft ) 


Thrust of the propellers (total) 


N. (lb) 




Maximum thickness ratio 

Average maximum thickness ratio of 
intersecting vertical-and horizontal 
tail surfaces 

Induced drag factor due to linear 
twist 

1) Width (diameter) of an equiva- 
lent circular body at the foregoing 
station being considered (Section 
4.7) 

<V ' ' 

2) 2ero-lift drag factor due to 
linear twist (Section- 4.12) 


w 


m (in, ft) 



LIST OF SYMBOLS (continued) 


Symbol 


w. 


B 


X 


X 

o 


X 


1 




/ 


X 

ac 


(x /c ) . . 
ac r f(w) 
6 


(x /c ) 
ac r w(£) 
e 


Definition Dimension 

Mean width of the body planform m (in, ft) 

segment, dx 

1) Distance from nose of body to 
the centroid of AS for the fore- 
body, and to the centroid of Ax 
of the afterbody (Section 4.7) 

2) Distance, parallel to the cen- 
terline of the wake, from the trail 
ing edge of the wing root chord to 
the quarter chord of the horizontal 
tail mean aerodynamic chord (Sec- 
tion 4.9) 

Distance from the nose of body to m (in, ft) 

the point where the potential flow 

ceases 


■ m (in, ft) 


m (in, .ft) 


1) Distance from the nose of body m (in, ft) 
to the point of raaximxnn negative 

rate of change of body cross-section- 
al area with body length (Section 4.3) 

2) Distance from the wing leading m- (in, ft) 
edge, to., ther^rentroid' of" the' 'forward 

Ax .segment- of' the body- planform area 
(Section. 4'j^)'- 

Distance from the wing trailing edge m (in, ft) 
to the centroid of the aft Ax seg- 
ment of the body planform area 

Length of the Ax segment of the body m (in, ft) 
planform area adjacent to and for- 
ward of the wing leading edge 

Distance from the lifting-surface m (in, ft) 
apex to the aerodynamic center of 
the surface 

Contribution to the aerodynamic cen- 
ter due to the lift carryover of the 
wing onto the fuselage, as a fraction 
of the root chord of the exposed wing 
panels 

Aerodynamic center of the wing with the 
wing in the presence of the fuselage, 
as a fraction of and about the leading 
edge of the root chord of the exposed 
wing panels 
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LIST OF SYMBOLS (continued) 


Symbol 


X 

eg 


x/c 


m 


X 

n, 


X 

P 


X ' 

P 


X 

w 




z/c 


Definition 


Distance to the center of gravity 
from the leading edge of the wing 
mean aerodynamic chord 

Section coordinate dimensions 

Distance, parallel to the X-rbody 
axis, from the quarter chord of the 
horizontal tail mean aerodynamic 
chord -to the leading edge of the 
wing mean aerodynamic chord 

Distance from the hose of body to 
the chosen moment center (leading 
edge of the mean aerodynamic chord 
in this case) 

Distance from the lifting surface 
apex to the desired reference cen- 
ter (leading edge of the mean aero- 
dynamic chord in this case) 

Distance from center of gravity to 
the propeller,, positive- forward 

Distance from quarter chord of the 
mean aerod3mamic chord of the im- 
mersed portion of the wing to the 
propeller, positive forward 

Distance from the aerodynamic cen- 
ter of the mean aerodynamic chord 
of the immersed portion of the wing 
to the center of gravity 

Lateral distance from the root chord 
to the mean aerodynamic chord 

Section coordinate dimens-ions- 


z^ Distance, parallel to the Z-body 

axis, from the X-body axis to the 
quarter chord of the horizontal 
tail mean aerodynamic chord, pos- 
itive down 


Dimension 
m (in, ft) 

m (in, ft) 

m (in, ft) 

m (in, ft) 

m (in, ft) 
m (in, ft) 

m (in, ft) 

m (in, ft) 

m (in, ft) 
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LIST OF S'jHBOLS (continued) 


Symbol 


z 


h 


I 





z 


s 


z 


T 


z 

w 


a 


Definition 


Vertical distance from the root- 
chord plane of the wing to the 
quarter chord point of the hori- 
zontal tail mean aerodynamic chord 

Vertical distance from the center- 
line of the wake to the quarter 
chord of the horizontal tail mean 
aerodynamic chord 

Effective distance, parallel to 
the Z-body axis, from the quarter 
chord of the horizontal tail mean 
aerodynamic chord to the centerline 
of the propeller slipstream, posi- 
tive down 

Distance, parallel to the Z-body 
axis, from the thrust axis to the 
quarter chord of the horizontal 
tail mean aerod3mamic chord, pos- 
itive down 

Distance, parallel to the Z-body 
axis, from the X-body axis to the 
centerline of the propeller slip- 
stream at the longitudinal station 
of the quarter chord of the mean 
aerodynamic chord of the immersed 
portion of the wing, positive down 

Distance, parallel to the Z-body 
axis , from the X-body axis • to the 
thrust axis, positive down 

Distance, parallel to the Z-body 
axis , from the X-body axis to the 
quarter chord of the mean aerodyn- 
amic chord of the immersed portion 
of the wing, positive down 

Angle of attack relative to X-body 
axis 


Limit of linearity of c 


xxxxxx 


Dimension 
m (in, ft) 


m (in, ft) 


m (in, ft) 


m (in, ft) 


m (in, ft) 


m (in, ft) 


m (in, ft) 


deg, rad 
deg 




LIST OF SYMBOLS ('continued) 


Symbol 

Definition 

Dimension 

“b 

Angle of attack of actual body, 
identical to airplane angle of 
attack, a, using X-body axis as 
reference 

deg 

“b .. 

eff 

Effective angle of attack of an 

equivalent circular body, a + 

ii' 

a 

°B 

Angle of attack of a surface at 
its 

max 

deg 

max 

deg 

^b(hf) 

■‘^max 

Angle of attack for maximum lift 
of the horizontal’ tail in the pres- 
ence of the fuselage 

deg 

max 

Angle of attack for maximum lift 
of the wing relative to the zero- 
lift line of wing 

deg 

max 

Angle of attack for maximum lift 
of wing-fuselage combination rel- 
ative to the zero-lift line of wing 

deg 

a 

^5, 

max 

Angle of attack at c^ 

max 

deg 

“h 

Local angle of attack of the hori- 
zontal tail with the stabilizer 
setting equal to zero 

deg 

a 

o 

Angle of attack for zero lift 

deg 

a 

Zero-lift angle of an equivalent 
circular body relative to the ref- 
erence X-hody axis- 

deg 

a 

o 

w 

Angle of attack of wing for zero 
lift relative to chord line 

deg 

a 

P 

Angle of attack of the propeller 
plane, includes the effect of the 
wing' upwash 

deg 

a 

s 

Average value of angle of attack in 
the nonlinear lift range to stall 

deg 
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LIST OF SYMBOLS (continued) 


Symbol 


Definition 


Dimension 


Angle of attack of the thrust deg 

axis 


a 

w 


Angle of attack of wing relative deg 

to its chord line, a + i 

w 


a * 
w 


Limit of linearity of the wing deg 

lift curve relative to chord, line 


a 

w 


abs 




(a J 
0 c. 


Wing angle of attack relative to 

the wing zero-lift line, 

a - a + i 
o w 
w 


C3C^/'36)/(3C^/9a) 


deg 


(3c^/36>/(3c^/3a) 


0 



3 ’ 


Propeller blade angle at 0.75 R 


deg 


•r 


Dihedral angle 


deg 


T 


cooling 

system 


Angle between the wing chord plane deg 
and the line connecting the trailing 
edge of the wing root chord and the 
quarter chord of the horizontal tail 
mean, aerodynamic chord 

Increment of drag coefficient due to 
cooling system 



Increment of drag due- to power effects 
on induced drag 



Increment of zero-lift drag due to 
power 


>h 

o 


Contribution of the horizontal tail 
to 


o 
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LIST OF SYMBOLS (continued) 


Symbol 


>h(f) 

o 


>n 

o 


D n(w) 
o 


(“d >t(f) 

O 


^v(f) 

o 


>v(h) 

O 


>w 

O 




AC, 


^^^L^h(fv) 






Definition 


Dimension 


Increment of zero-lift drag of the 
horizontal tail dtie to fuselage 
interference 

Contribution of nacelles (including 
nacelle-wing interference effects) 
to 4Cj, 
o 

Increment of the zero-lift drag of 
one nacelle due to wing interference 

Increment of the zero-lift tail drag 
(horizontal-or vertical tail) due to 
fuselage interference 

Increment of the zero-lift drag coef- 
ficient of the vertical tail due to 
fuselage interference 

Increment of the zero-lift drag coef- 
ficient of the vertical tail due to 
the horizontal tail interference when 
the vertical tail intersects the hori- 
zontal tail 

Contribution of the wing to 

o 

Component of the total thrust para- 
llel to the velocity vector, positive 
thrust is equal to a negative drag 
contribution 

Increment of lift 


Increment of lift coefficient of the 
horizontal tail due to the effect of 
fuselage vortices 

Increment of lift coefficient due to 
stabilizer deflection 

Increment of horizontal tail contribu- 
tion to the lift coefficient resulting 
from the power-induced change in dynam- 
ic pressure at the tail 
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LIST OF SYMBOLS (continued) 

Symbol 

Definition Dimension 

h h power 

Increment of horizontal tail con- 
tribution to the lift coefficient 
resulting from the power-induced 
change in downwash at the tail 

max 

1) Mach number- correction of the 
incompressible, maximum lift coef- 
ficient (Section 4.2) 


2) Increment of maximum lift coef- 
ficient due to power (Section 5.1) 

max tab 

Increment of maximum lift coeffic- 
ient due to the tab 

)-■ 

Ij 'POWfilT 

Difference, between predicted power- 
on lift and predicted propeller-off 
lift of the complete airplane 

(AC^ ) ’ , 

L pox^er 

Increment of the tail-off lift due 
to power at the propeller-off max- 
imum lift angle of attack 


Increment of lift coefficient due 
to the power-induced change in dyn- 
amic pressure over the portion of 
the wing immersed in the propeller 
slipstreams 

»Vt 

Increment of lift coefficient due to 
the lift component propeller thrust 
vector 

>w 
a s 

Correction to reduce (C.^ ) - to an deg rad ^ 

L^ wfn - “ 

average slope in the nonlinear lift 

range to stall 

<"Ve 

P 

Increment of lift coefficient due to 
the change in angle of attack, resul- 
ting from propeller doTmwash, e , of" 
the portions of the wing immersld in 
the propeller slipstreams 

P 

Increment of pitching moment coef- 
ficient due to propeller effects on 
body (fuselage or nacelle) 

<3 

Increment of. tail contribution to the 
pitching moment coefficient due to the 
propeller induced increments of dynamic 
pressure and downwash at the tail 
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LIST OF SYMBOLS (continued) 


Symbol 


Definition 


Dimension 




m n 


m r 
o 


(AC ),- 
ni 

o 




AC /e 

m 

o 


(AG ) 

m power 


<“m)T 


(AC ) 
m w. 


(AC ).- 
m Aq, 


(AG ).- 
m Aq 


w 


(AC ),. s 
m*^ (A^h^ 


power 


Increment of pitching moment due to 
the propeller normal force 

Increment of pitching moment due to 
the propeller effects on the nacelles 

Correction to (C )- for the non- 
m X 
o 

midwing, configuration 

Increment of zero-lift pitching mo- 
ment coefficient due. to the increase 
in dynamic pressure induced by the 
propeller slipstreams on the immersed 
portions of the wing 

Change in wing zero-lift pitching deg 
moment coefficient due to a unit 
change in linear wing twist 

Difference between predicted power-on 
pitching moment coefficient and pre- 
dicted propeller-off pitching moment 
coefficient of the complete airplane 

Increment of pitching moment coef- 
ficient due to the thrust of the 
propellers 

Increment of pitching moment coef- 
ficient due to the net change in wing 
lift resulting from propeller-slip- 
stream- induced dynamic pressure and 
angle-of-attack changes on the wing 

Increment of tail contribution to the 
pitching moment coefficient due to 
the propeller induced increment of 
dynamic pressure at the tail 

Increment of pitching moment coef- 
ficient due to the increase in 
dynamic pressure induced by the 
propeller slipstreams on the immersed 
portions of the wing 

Increment of tail contribution to the 
pitching moment coefficient due to ’ 
the propeller induced increment at 
downwash of the tail 
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LIST OF SYMBOLS (continued) 


Symbol 


Definition 


m £ 


Ac, 


max 


(Ac^ ) 
max 


(Ac^ ) 

max 




base 


tab 


(Ai/ij^ 


w n 


AS 

X 


Ax 

Ay 

^^wake 


Change in the pitching moment . coef- 
ficient of the wing due to the 
propeller slipstream induced 
change in angle of attack 


Correction of c for Reynolds 
max 


number 


Section maximum lift increment for 
the 25-percent-chord tab (flap) at 
a reference tab-deflection angle 
of 60 degrees 

Section maximum lift increment due 
to the tab 

Increment of power-induced dynamic 
pressure acting on the horizontal 
tail as a ratio of the free-stream 
dynamic pressure 

Dynamic pressure loss at the hori- 
zontal tail as a ratio of the free- 
stream dynamic pressure 


Dynamic pressure loss at the wake 
centerline as a ratio of the. free 
stream dynamic pressure 


Increment of power-induced dynamic 
pressure acting on the wing Immersed 
in the propeller slipstream as a ra- 
tio of the free stream d 3 m.amic 
pressure 

Wing area overlapped by one nacelle 


Change in the cross-section area of 
the body across the Ax segment of 
the body segment considered 

Incremental length of the body 

Section leading edge sharpness 
parameter, percent of chord 

Half-width of the wake at a distance 
X from trailing edge of the wing 
root chord 

xxxix 


Dimension 


2 2 2 
m' (in , ft ) 

m (in , ft ) 
m (in, ft) 

m (in, ft) 




LIST OF SYMBOLS (continued) 


Symbol 

Definition 

Dimension 

Aa 

max 

Angle of attack correction at 
for flow separation 

max 

deg 

''“s. 

X 

Change in angle of attack of the 
portion of the wing immersed in 
the propeller slipstream due to 
the propeller 

deg 

ha /6 
o 

Shift in a per degree of wing 
twist 


(As, ) 

b power 

Increment of downwash at the 
horizontal tail due to power 

deg 

6 

Deflection 

deg 

9S/3a 

Variation of upwash and downwash 
with angle of attack at the Ax 
segment of the body forward of the 
wing leading edge and aft of the 
wing trailing edge 


83/3a 

Variation of upwash with angle of 
attack of the Ax segment of the 
body forward of and adjacent to 
the wing leading edge 


3s, /-3a 
n. 

Average downwash gradient at and 
across the horizontal tail 


3e /3a 
P 

Downwash gradient behind the pro- 
peller 


-Oe^Sa) 

Upwash gradient at the propeller 


3e / 3a 

VC 

Rate of change of downwash, in the 
plane of symmetry at the height of 
the vortex core, with the absolute 
angle of attack 


3e /3a 
w 

Rate of change of downwash, behind 
the wing, with angle of attack 


3e /.3a 

CO 

Downwash gradient at infinity 



Twist at spanwise station p 

deg 
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LIST OF SYMBOLS (continued) 


Symbol 

^ . Definition 

Dimension 

£ 

Effective propeller downwash angle 
average over entire wing 

deg 


Downwash in the plane. of symmetry 
at the vortex sheet 

deg 

"h 

Average downwash across the hori- 
zontal tail 

deg 

£ 

P 

Downwash angle behind the propeller 

deg 

-e 

u 

Upwash angle at the propeller 

deg 

ri ■ 

1) Drag proportionality factor 
(Section 4.3) 

2) Non-dimensional span-wise station 
y/(b/2) (Section 4.2) 

5 


Distance from centerline of tail to 
the inboard tab edge as a ratio of 
the tail semispan 


% 

Distance from centerline of tail to 
the outboard edge of the tab as a 
ratio of the tail semispan 


9 

Twist of the wing tip with respect 
to the root chord (negative for 
washout) 

deg 

A 

Sweep angle 

deg 

' 

Corrected sweep angle, 
tan"^(tan(A^^^)/3] 

deg 

A 

Taper ratio, c^/c 
t r 


^te 

Section trailing-edge angle 

deg 

Subscripts 



B 

Body 


c/2 

Half -chord line 


c/4 

Quarter-chord line 
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LIST OF SYMBOLS Ccontiaued) 


Subscripts 

Definition 

eg 

Center of gravity 

cooling system 

Engine cooling system (cowl flaps 
and engine inlets) 

e 

Exposed panels 

f 

1) Fuselage 

2) Flap (Section 4.13) 

fn 

Fuselage-nacelle combination 

h 

Horizontal tail 

i 

Immersed in propeller slipstream 

o 

II 

Stabilizer not deflected 

£ e 

Leading edge 

low speed 

Low subsonic Mach number 

£s 

Lifting surface 

M 

High subsonic Mach number 

max 

Maximum 

n 

Nacelle (s) 

prop off 

Propeller (s) removed 

power on 

Propeller (s) installed and power on 

r 

Root 

s 

Nonlinear lift range to stall 

t 

Tip 

tab 

Tab 

V 

Vertical tail 

w 

Wing 

wf 

Wing-fuselage combination 
. xlii 
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^Subscripts 

wfn 


^tab 


= 0 


Definition 

Wing-fuselage-nacelle combination 
Tab deflection is zero 


xliii 
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CHAPTER 1 


INTRODUCTION 

The Advanced Technology Light Twin (ATLIT) airplane was developed 
by the University of Kansas Flight Research Laboratory as part of a 
general aviation research program sponsored by the National Aeronautics 
and Space Administration, Langley Research Center. Background infor- 
mation on this program together with the history of the ATLIT develop- 
ment and performance predictions are given in Reference- 1. The broad 
objective of this research program was to apply existing jet-transport 
wing technology to general aviation airplanes for the purpose of 
improving safety, efficiency and utility. 

The ATLIT project was managed by the University of Kansas and, 
in addition, involved the Robertson Aircraft Company for detailed 
design; the Piper Aircraft Company for modification, fabrication and 
initial flight testing; and Wichita State University for ^rind tunnel 
tests in support of the spoiler lateral-control-system development. 

The Piper PA-34-200 Seneca I was selected for modification 
as a typical example of current twin-engine general aviation aircraft 
which are virtually all limited in single-engine climb performance at 
gross weight. As pointed out in Reference 1, the conceptual study 
which led to the ATLIT design showed a potential for the much needed 
improvement in single-engine climb performance without increasing 
installed power or reducing gross-weight. Improved single-engine 
climb performance became the major goal of the ATLIT project. 

The preliminary flight test results (Reference 1) indicated that 


the ATLIT failed to achieve the predicted improvements in climb and 



cruise performance. They were not better than those of the basic 
Seneca at the same gross-weight and with the same installed power. 

A full-scale xri.nd tunnel investigation of the ATLIT airplane was 
undertaken to evaluate the various advanced aerodynamic concepts and 
to determine the cause for the lack of performance improvement. The 
results of this wind tunnel investigation are presented in Reference 2. 

In this report the first phase will be described of a project 
performed by the Flight Research Laboratory of the University of 
Kansas sponsored by Grant NSG 1574 from NASA, Langley Research Center. 
The objectives of this project are to 

1. correlate theoretical predictions of longitudinal aerodynamic 
characteristics with full-scale wind tunnel data; 

2. correlate theoretical predictions of lateral directional 
aerodynamic characteristics with full-scale wind tunnel 
data; 

3. correlate the results of point 1 and point 2 with flight 
test data. 

In this report only work toward the first objective will be 
discussed. References 3 and 4 will be used as guidelines for the 
theoretical predictions. Use will also be made of References 5 and 6, 
and the results will be compared with those obtained with References 
3 and 4 and the experimental data (Reference 2) . Also the results 
mentioned in Reference 7 will be discussed. 
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CHAPTER 2 


THE ATLIT AIRPLANE 

The airplane used in the analysis is the full-scale ATLIT 
(Advanced Technology Light Twin) which is an extensively modified 
Piper PA-34-200 Seneca I general aviation low-wing monoplane with 
an all-movable horizontal stabilizer. Pertinent physical character- 
istics are listed in Table 2.1 and a three-view drawing is presented 
in Figure 2.1. 

The advanced technology improvements implemented on the ATLIT 

were 

1. replacement of the basic untapered, aspect ratio 7.25 Seneca 

2 2 

wing having a 652415 airfoil and an area of 19.4 m (208.7 ft ) 
with a newly designed tapered (taper ratio 0,5), high aspect 
ratio (aspect ratio 10.32) wing having a 17-percent thick 
GA(W)-1 airfoil and an area of 14.4 m^ (155.0 f t^) ; 

2. installation of full-span 30-percent chord Fowl^“ flaps 
instead of the partial span 20-percent chord plain flaps 
used on the Seneca; 

3. use of a spoiler lateral control system instead of conven- 
tional ailerons. 

The planform changes were made to lower the induced drag, while the 
GA(W)-1 airfoil was chosen for its high lift-to-drag ratio and for 
its high maximum lift coefficient. Also supercritical propellers, 
designed for increased propulsive efficiency, were used. These, 
however, were not installed during the wind tunnel investigation. 
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Special thrust-torque balances were installed on the propeller 
shafts to measure the propeller characteristics. The balances made 
it necessary to extend the nacelles with 8 inches. 

As already mentioned in Chapter 1, the climb and top speed 
performance of the ATLIT fell short of predictions. This was the 
reason that part of the full-scale wind tunnel investigation was 
devoted to drag evaluation of the airplane "as built." A drag 
clean-up investigation was initiated and several factors were estab- 
lished which adversely affected the lift and drag of the airplane. 
These factors are discussed in detail in Reference 2. The result of 
the drag clean-up program was a "fully clean" configuration of the 
ATLIT airplane. Comparisons between experimental and theoretical 
results will be based on the ATLIT "fully clean." 
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Table 2,1: Specifications of the ATLIT airplane 


Wing ; 


Area, 

14.40 (155.0) 

Span, m (ft) 

12.19 (40.0) 

Aspect ratio 

10.32 

Thiclcness ratio 

0.17 

Dihedral, deg 

7 

Taper ratio 

0,5 

Incidence angle at root, deg 

0.5 

Incidence angle at tip, deg 

-2.5 

Leading edge and trailing edge sweep, deg 

3.67 

Mean aerodynamic chord, n (ft) 

1.225 (4.018) 

Root chord, m (in) 

1.575 (62.0) 

Tip chord, m (in) 

0.787 (31.0) 

Airfoil 

GA(W)-1 

Flap; 

Type 

Fowler 

Area (total), (ft^) 

3.56 (38.3) 

Span/side, m (ft) 

5.15 (16.91) 

Chord, perc. of wing chord 

30 

Inboard wing station, m (in) 

0.711 (28.0) 

Outboard wing station, m (in) 

5.867 (231.0) 

Leading edge retracted, perc. of wing chord 

70 

Maximum deflection, deg 

40 

Maximian travel, m (in) 

0.343 (13.50) 

Root chord, n (in) 

0.445 (17.50) 

Tip chord, m (in) 

0.244 (9.62) 
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Table 2.1; Continued 



Area (total), (ft^) 
Spaa/slda, m (in) 

Inboard wing station, m (in) 
Outboard wing station, m (in) 
Binge, perc. of wing chord 
Maximum deflection, deg 
Root chord, m (in) 

Tip chord, m (in) 


Tail ; 

Stabilator: 

Area (Including tab), (ft^) 
Span, m (ft) 

Chord (constant) , m (in) 
Aspect ratio 
Sweep angle, deg 
Dihedral, deg 
Hinge line, perc, of chord 
Airfoil 
Vertical : 

Area (including cab) , (ft^) 

Span, m (ft) 

Aspect ratio 
Taper ratio 
Root chord, m (in) 

Tip chord, m (in) 


Differential-Slat- 
Dppar Surface 

0.488 (5.25) - 

3,226 (127.0) 

2.565 (101.0) 

5.791 (228,0) 

70 

60 

0.177 (6.975) 

0.124 (4,880) 


3.60 (38.7) 
4.13 (13.56) 
0.371 (34.29) 
4.75 
0 
0 

26.6 

NACA 0010 

1.85 (19.9) 
1.52 (5.00) 
1.19 
0.42 

1.708 (67.23) 
0.723 (28.45) 
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Table 2,1; Continued 


Mean aerodynamic chord , m (in) 

1.282 

(50.47) 

Leading edge sweep angle, deg 

39.92 


Axrfoli 

NACA 0009 

Stabilator tab:- 

. 


Area, (ft^) 

0.539 

(5.8) 

Span, m (in) 

3.023 

(119.0) 

Chord (constant), m (in) 

0.178 

(7.0) 

Tab hinge line to stabilator hinge line, 
m (in) 

0.462 

(18.17) 

Rudder: 

- 


Area (including tab) , (ft^) 

0.706 

(7.6) 

Span, m (in) 

1.525 

(60.02) 

Chord (constant) , m- (in) 

0.449 

(17.69) 

Rudder trim tab: 



2 2 
Area,, m (ft ) 

0.158 

(1.7) 

Span, m (in) 

0.895 

(35.25) 

Chord (constant) . m (in) 

0.175 

(6.9) 


2 

Lycoming 

I0-360-CIE6 

2700 

200 

HarCzall 

HC-C2YK-ICEF/F0 7666A 
2 

1.930 (76.0) 


Propeller : 

Manufacturer 

Model 

Number of blades 
Diameter , m (in) 
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Table 2.1: Concluded 


Weights and I.oading i 
Gross weight, S (lb) 

Empty weight, N (lb) 

Useful load, N (lb) 

Wing loading (at gross weight), 
N/m^ (Ib/ft^) 


17,792 ( 4 , 000 ) 
11,018 ( 2 , 477 ) 
- 6,774 ( 1 , 523 ) 

1,236 ( 25 . 8 ) 


8 



! m ^ ^ 3.01 m'^ 


P 8.72 m 

NORMAL GROUND LINE 



Figure 2.1: Three-view drawing of the ATLIT airplane (Reference 1) 





2.1 Geometric Parameters of Wing and Tail 


Depending on which characteristic has to be determined, the total 
area or the exposed area of the wing and the tail are considered. The 
total planform is considered to extend through the nacelle and the 
fuselage, while the exposed planform terminates at the fuselage. 
Pertinent dimensions for the wing, the horizontal tail and the vertical 
tail are shown in Figures 2.1.1, 2.1.2 and 2.1.3, respectively. 

Table 2.1.1 lists the geometric parameters of the wing and the tail 
pertinent in the analysis. 
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Table 2.1.1: Pertinent wing and tail geometric parameters used in the analysis 




Wing 

Horizontal Tall 

Vertical Tall 

Symbol 

Dcscripcion 

Total 

Exposed 

Total 

Exposed 

Exposed 

S 

Area, (ft^) 

14.40 

(155.0) 

12.53 

(134.8) 

3.60 

(38.7) 

3.25 

(34.9) 

1,75 (18.8) 

b 

Spaa, ra (ft) 

12.19 

(40.0) 

10.96 

(36.0) 

4.13 

(13.56) 

3.73 

(12.23) 

1.52 (5.0) 

A 

2 

Aspect ratio 9 b /S 

10.32 


9.61 


4.75 


4.28 


1.33 


■rip chord, m (In) 

0.787 

(31.00) • 

0.787 

(31.00) 

0.871 

(34.29) 

0,871 

(34.29) 

0.723 (28.45) 


Root chord, m (in) 

1.575 

(62.00) 

1.495 

(58.87) 

0.871 

(34.29) 

0.871 

(34.29) 

1.575 (62.0) 

\ 

Taper ratio, Cj./c^ 

0.50 


0.527 


1 


1 


0.459 

C 

* Mean aerodynamic 
chord, m (In) 

1.225 

(48.22) 

1.178 

(46.38) 

0,871 

(34.29) 

0.871 

(34.29) 

1.201 (47.30) 

^3 

** Lateral position of 
mean aerodynamic chord, 
m (in) 

2.709 

(106.67) 

2.46 (96.85) 

1.033 

(40.68) 

0.932 

(36.69) 

0.334 (13.15) 

r 

Dihedral angle, deg. 

7 


7 


0 


0 



^le 

Leading-edge sweep, 
deg. 

3.67 


3.67 


0 


0 


40.0 

^c/4 

Sweep of c/4 line, , 
deg. 

1.835 


1.835 


0 


0 


34.5 

^c/2 

Sweep of c/2 line, 
deg. 

0 


0 


0 


0 


29.0 




** 


^3 


^ 1 . 1 + 2 \ b 
3 ^ 1 + A ^2 




Figure 2.1.1: Definition sketch, of wing dimensions 
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CHAPTER_3 

LISTING OF IMPORTANT RESULTS 

In this chapter the calculated longitudinal characteristics of the 
ATLIT airplane will.be compared with full-scale wind tunnel data of 
Reference 2. The results are only discussed briefly, a more detailed 
discussion of the results and how to improve the results can be found 
in Chapters 4 through 6. 

3.1 Lift Characteristics of the ATLIT Airplane 

The lift of the airplane without horizontal tail is discussed in 
Section 4.1 to 4.4. The calculated lift curve of the ATLIT without 
horizontal tail is shown in Figure 3.1.1 and compared xjith wind tunnel' 
results. In this figure results obtained with Reference 5 are also 
shown. These references both represent lifting surface prediction 
methods. Reference 5 is based on the thin wing theory, while Reference 
6 is based on the thick wing theory. The calculated lift curve shows 
fair agreement with the experimental results, except for the discrepancy 
in the, zero-lift angle of attack. 

The lift of the complete airplane is discussed in Section 4.10 and 
the calculated lift curve is shown in Figure 3.1,2. The predicted lift 
curve shows fair agreement with the wind tunnel data. No results were 
obtained with the computer program of Reference 6, 

The power-on lift of the airplane is determined in Section 5.1. 

The predicted lift curve is drawn for T^’ = 0.0915 and 0.1970 in Figure 
3.1,3 and 3.1.4, respectively. The calculated lift curves show poor 
agreement with the experimental results. This is caused, however, by 
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the difference between the predicted and experimental propeller-off 


lift cuirve. When the predicted increment in lift due to power, 

(AC^ ) -j,is added to the experimental lift curve (propellers off, 

L pred 

power 

cowl flaps and engine inlets open) good agreement is obtained with the 
wind tunnel data. 


3.2 Pitching Moment Characteristics of the ATLIT Airplane 

The pitching moment of the airplane without horizontal tail is 
described in Section 4,5 through 4.8. In Figure 3.2.1 the predicted 
horizontal-tail-off pitching moment is compared with experimental re- 
sults of Reference 2 and with results obtained with Reference 5 and 
Reference 6. The calculated results show good agreement with the wind 
tunnel data. 

The pitching moment of the complete airplane is discussed in Sec- 
tion 4.11. The calculated results are plotted in Figures 3.2.2 and 3.2.3 
and they show good agreement with experimental data. 

In Section 5.2 the power-on pitching moment is described and the 
results are shown in Figures 3.2.4 through 3.2.7 for two power settings. 
The calculated characteristics show fair agreement with the wind tunnel 
data. The predictions have been performed for a Reynolds number of 2.3 
million. No power-on wind tunnel data, however, were available for the 
ATLIT in the "fully clean" configuration at a Reynolds number of 2.3 
million. Data, were available for a Reynolds number of 3,5 million. In 
the case of the ATLIT, this increase in Reynolds number resulted in an 

increase in pitching moment (AC^)^^ — 0.03, in the linear lift region. 

^Re 

The predicted pitching moment, including the Reynolds nttmber correction. 
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is also shown in Figures 3.2.4 through 3.2.7. the predicted move- 


ment in pitching moment due to power, (AC ) added to the 

power ^ 

experimental pitching moment curve (propellers off, cowl flaps and en- 
gine inlets open) good agreement is obtained with the wind-tunnel data. 


3.3 Drag Characteristics of the ATLIT Airplane 

The drag of the airplane is discussed in Section 4.12, while the 
power-on drag is described in Section 5.3. In both cases the calculated 
drag polar shows good agreement with the wind tunnel data, as shown in 
Figures 3.3.1 through 3.3.4. 
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Figure 3.1,1: Comparison, of lift predictions with wind tunnel data 

2.3x10^, horizontal tail and propellers removed) 
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Figure 3.1.2;’ Comparison of predicted airplane lift curve with 

wind tunnel data (propellers removed, million) 
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Figure 3.1.4; Comparison o£ predicted lift curve with full-scale 
wind tunnel data (T ’ = 0,1970, no_ stabilizer 
deflection) 
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Figure 3,2,3: Comparison of predicted airplane pitching moment with 

f^l^~scale wind tunnel data (propellers removed, stabilizer 
not deflected, 2.3 million) 
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RElERODUCIBILnY OP THJJ 
ORlG^jiL PAGi<j IS POOR 



Figure 3.2,4; Comparison of predicted pitching moments with 


full-scale wind-tunnel data (T ' = 0,0915, no 
stabilizer deflection) 









M 
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Figure 3.2,7; 


Comparison of predicted pitching moments with 
experimental results (T^* = 0,1970, no stabilizer 
deflection) 




Figure 3.3.1: Comparison of predicted airplane drag with full- 

scale wind tunnel data (propellers removed, no 
^ engine cooling drag, 2.3 million) 


EEPR0DUCF*F.ITY" OF TRW 
ORIGP"^ L . IS'POOP 










Figure 3.3.2: Comparison of predicted airplane drag with full- 

scale wind tunnel data (propellers removed, no 
engine cooling drag, 2.3 million) 
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Figure- 3.3.3: Comparison of calculated and wind- tunnel-determined 

drag at different power conditions 


REPEODUCIBIL'ITY OF THE) 
ORIGINAL, TS POOR 
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Figure 3.3.4: Comparison of predicted drag, with full-scale wind 

tunnel data at different power conditions 
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CHAPTER 4 


PREDICTION OF PROPELLER-OFF AERODYNAMIC CHARACTERISTICS 

In this chapter the propeller-off aerodynamic characteristics 
will be discussed. Mainly, the method of Reference 3 will be used. 
However, Reference 3 does not discuss the influence of wing twist 
on the lift, pitching moment- and drag of the airplane. References 
4 and -8 will be used to incorporate this effect. 

4.1 Wing and Horizontal Tail Airfoil Section Characteristics 

A great deal of theoretical and experimental work has been done 
toward the development of airfoil sections. The theoretical airfoil 
design, however, is hampered by the viscous effects in the form of 
the boundary layer between the airfoil surface and the free stream. 
This boundary layer affects the section drag and maximum lift coef- 
ficient and also the slope of the lift curve, the angle of attack 
for zero lift and the section pitching-moment coefficient. Since 
the boundary layer is influenced by the surface roughness , the cur- 
vature of the surface, the pressure gradient, heat transfer between 
the surface and the boundary layer and viscous interaction with the 
free stream, it is apparent that no simple theoretical considerations 
can accurately predict ail the airfoil characteristics. For these 
reasons , experimental data are always preferable to theoretical 
calculations . 

Table 4.1.1 stimmarizes experimental data for the NACA four- and 
five-digit airfoils and for the NACA 6-series airfoils. The data, 
from Reference 3, are for sraooth-leading-edge conditions and 9' x 10®' 
Reynolds number. Information is presented on the following airfoil 


33 - 




characteristics ; 


1. angle of attack for zero-lift, a 

o 

2. lift curve slope, c^ 

X/ 

Cl 

3. angle of attack at which the lift curve deviates from 

* 

linear variation, a 

4. maximum lift coefficient, c. 

max' 

5. angle of attack for maximum lift, a 

max 

6. zero-lift pitching moment coefficient, c 

m 

o 

7 • position of the aerodynamic center as a ratio of the 

chord length, a.c. 

From the first five quantities the approximate lift-curve shape can 
be synthesized. Experimental data for a large number of additional 
airfoils are available in the literature (see Reference 4) . 

In this section a theoretical approach will not be considered. 
Section characteristics will be based on experimental data with the 
maximum lift coefficient corrected to the Reynolds number being 
considered'. The effect of Reynold.s number on the maximum lift coef- • 
ficient can be accounted for by using, Figure 4.1.1, which uses the 
leading-edge sharpness parameter. Ay, as the correlating parameter. 

The leading-edge sharpness parameter is defined in Figure 4.1.2. 

From this figure, the leading-edge sharpness parameter can be obtained 
as a function of the airfoil type and the thickness ratio. Another 
important parameter is the airfoil trailing-edge angle, <p^^, which 
can be obtained from Figure 4.1.3. 

In the case of the ATLIT airplane. Reference 9 has been used to 
obtain the section airfoil characteristics of the wing.' Figure 4.1.4 
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shows the section lift coefficient and the section pitching moment 
coefficient of the GA(W)-1 airfoil as being used in this study. 

The section airfoil characteristics of the horizontal tail are 
determined from Table 4.1.1. The leading-edge sharpness parameter 
and the trailing-edge angle of the horizontal tail (NACA 0010) air- 
foil can be obtained from Figures 4.1.2 and 4.1.3, respectively. 

For the GA(W)-1 airfoil, however, these parameters have to be ob- 
tained from the section shape of the airfoil (Figure 4.1.5). Table 
4.1.2 summarizes the airfoil section characteristics of the wing and 
horizontal-tail . 
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Table 4.1.1: Experimental low~speed airfoil section aerodynamic characteristics 
(Reference 3) 


(^) 5-<iigit airfoils, Ngg « 9 X 10^, smooth leading edge 



(b) 6-series airfoUs, Ngg = 9 x 10® , smooth leading edge 


Airfoil 

(i3-00G 

63-009 

63-206 

63-209 

63-210 

63^-012 

63^-212 

63j^-412 

639-015 
632-215 
632— tlo 

632- 615 

633- 018 
633-218 
633-418 

633- 618 

634- 021 
634-221 
634-421 

63.4- 120 

63.4- 420 a = 
63 (420) -422 
G3(420)-517 

64-006 

64-009 

64-108 

64-110 

64-206 

64-208 

64-209 

64-210 



0 -.015 

- 1 . 0 020 


.110 .255 10.0 1.10 10.0 

.110 .261 13.0 1.40 10.0 



-.040 
-. 039 
-.040 
-.040 


110 .253 12.0 
113 .257 10.5 
107 .261 13.0 
110 .238 14.0 


1.03 8.0 

1.23 3.3 

1.40 8.9 

1.45 10.3 


64 1 - 0 12 


64j-112 

-.8 

64,-212 

-1.3 

64^412 

-2.6 


0 .111 .262 

-.017 .113 .267 

-.027 .113 .262 

-.065 .112 .267 


14.5 1.45 11.0 

14.0 1.30 12.2 

15.0 1.55 11.0 

15.0 1.67 b.O 
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Table 4,1,1: Concluded 


Airfoil 



=fo,per deg 

a.o. 

^C, , deg 

^max 

‘’max 

or, deg 

642-015- 

0 

0 

0. 112 

0.267 

15.0 

1.48 

13.0 

G42"215 

-1.6 

030 

.112 

.265 

15.0 

1, 57 

10, 0 

642-410 

-2.8 

-. 070 

. 115 ■ 

.264 

15.0 

1. 65 

8.0 

643-018 

643-218 

0 

-1.3 

.004 
-. 027 

.111 

.115 

.266 

.271 

17.0 

16.0 

1.50 

1.53 

12.0 

10.0 

643-418 

-2,9 

-. 065 

.116 

.273 

14.0 

1.57 

s.o 

643-618 

-3.8 

095 

.116 

.273 

16.0 

1.58 

3.6 

644-021 

0 

.005 

.110 

.274 

14.0 

1,30 

10.3 

644-221 

-1.2 

-.029 

.117 

.271 

13.0 

1.32 

6.3 

644-421 

-2.8 

-. 068 

.120 

.276 

13.0 , 

1.42 

6.4 

65-006 

0 

0 

.105 

.258 

12.0 

.92 

7.6 

65-009 


0 

.107 

.264 

11.0 

1. 08 

9.8 

65-206 - 

-1.6 

-.031 

.105 

.257 

12.'0 

1,03 

G.O 

65-209 

-1.2 

031 

.106 

.259 

12.0 

1.30 

10.0 

65-210 

-1.6 

-.034 

.108 

.262 

13.0 

1.40 

9.6 

65-110 

-2.5 

-.067 

.112 

.262 

14.0 

1.52 

8.0 

65|-0X2 

0 

0 

.110 

.261 

14.0 

1.36 

10.0 

65,-212 

-1.0 

-. 032 

.108 

.261 

14.0 

1.47 

9.4 

a = .6 

-1.4 

, -.033 

.108 

.269 

14.0 

1.50 

9.6 

65^-4 iZ 

-3.0 

-. 070 

■ 111 

.265 

15.5 

1.66 

10.5 

602-015 

0 

0 

.111 

.257 

15.0 

1.42 

11.2 

G52“2X5 

-1.2 

-. 032 

.112 

.269 

15.5 

f.53 

10.0 

G52“4X5 

-2.6 

-.060 

■ 111 

.268 

16.0 

1. 61 

8.7 

652“4X5 a = .5 

-2.6 

-.051 

-111 

,264 

20.0- 

1. 60 

7,0 

■6o(21S)-114. 

-.7 

0 

.112 

.265 

15.0 

1.44 

. 10.5 

65(216) -415 a = .5 

-3.0 

-.057 

.106 

.267 

18.0 

1. 60 

6.0 

65,3-018 

0 

0 

.100 , 

.262 

17.0 

1.44 

10.0 

65-418 a = .8 

-3.0 

-. 081 

.112 

.266 

20.0 

1, 58 

4.4 

65-618 

-t.O 

100 

.110 

.273 

20.0 

1. 60 

4.9 

6S3-018 

0 

0 

.100 

.267 

16.0 

1.37 

10.0 

603-218 

-1.2 

-.030 

.100. 

.263 

18.0 

1.48 

8.8 

653-418 

-2.4 

-. 059 

.110 

.265 

18.0 

1.54 

4.9 

65|j-418 a = .3 

-2.8 

055 

.115 

.267 

18.0 

1.50 

6.0 

653-618 

.4.0 

-.102 

.113 

,276 

18.0 

1.64 

5.2 

6&3-6I8 a » « 5 

-4.2 

-.078 

,104 

.265 

20.0 

1.51 

5.3 

654-021 

-0 

0 

.112 

.267 

18.5 

1.40 

7.4 

604-221 

-1.3 

-.029 

.115 

.274 

20,5 

1.46 

6.0 

604-421 

-2.8 

-.066 

,116 

.272 

22.0 

1.56 

S.O 

G54“i2l a — • 5 

-2.8 

-. 052 

.116 

.272 

20.0 

1.43 

5.6 

6S(421)-420 

-2.4 

-. 061 

.116 

.276 

' 20.0 

1.52 

4,7 

66-006 

0 

0 

.100 

.252 

9.0 

.80 

6.5 

66-009 

0 

0 

.103 

.259 

10.0 

1.05 

10.0 

66-206 

-1.6 

-.038 

.108 

.257 

10.5 

1.00 

7.0 

66-209 

-1.0 

-.034 

.107 

.257 

11.0 

1.17 

9.0 

66-210 

-1.3 

-.035 

.110 

.261 

11.0 

1.27 

10. 0 

6G4-012 

0 

0 

.106 

.258 

14.0 

1.25 

11.2 

664-212 

-1.2 

-. 032 

.102 

.259 

IS.O 

1.46 

11.6 

662-015 

0 

.005 

.105 

.265 

IS. 5 

1.35 

12.0 . 

660-2 IS 

-1.3 

-.031 

.106 

.260 

16.0 

1. 50 

11.4 

662-415 

-2.6 

-. 069 

.106 

,260 

17.0 

1, 60 

lO-O 

66(215)-016 

0 ■ 

0 

.105 

.260 

14.0 • 

1,33 

10.0 

66(215)-216 

-2.0 

-.044 

.114 

.262 

16.0 

1. 55 

8.8 

66(215)-216 a= .6 

-1.2 

-. 030 

.100 

.257 

16.0 

1.46 

7.0 

66(215)-416 

-2.6 

-.068 

.100 

.265 

18.0 

1.60' 

4.0 

63A010 

0 

.005 

.105 

.254 

13.0 

1.20 

10. 0 

63A210 

-1.5 

-.040 

.103 

.257 

14. O' 

1.43 

10.0 

64AOIO 

0 

0 

.110 

.253 

12.0 

1.23 

10.0 

64A210 

-1.5 

-. 040 

.105 

.251 

3.0 

1.44 

10,0 

64A410 

-3.0 

-. 080 

.100 

.254 

15.0 

1.61 

10.0 

i644A212 

-2.0 

040 

,100 

.252 

14.0 

1.54 

11.0 

G42A.Z1£) 

-2.0 

-. 040 

.095 

.252 

15.0 

1.50 

12.0 
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Table 4.1.2: Airplane wing and horizontal tall airfoil section characteristics 







Horizontal 

Symbol 

DoscrlpCion 

Reference 

Wing 

Reference 

Tail 

- 

Airfoil oeccion 

Table 2.1 

GA(H)-i 

Table 2.1 

NACA 0010 

t/c 

Thickness ratio 

Table 2.1 

0.17 

Table 2.1 

0.10 

Ay 

Leading-edge-sharpnees parameter 

Figure 4.1.5 

4.9 

Figure 4.1.2 

2.6 

'*'te 

Trailing-adge-angle, deg 

Figure 4.1.5 

14.5 

Figure 4.1.3 

13.0 

a 

o 

Zero-lift a relative to chord line, 
dag 

Figure 4.1.4 

-3.7 

Table 4.1.1 

0 


Lift-curve slope, per deg 

Figure 4.1.4 

0.115 

Table 4.1.1 

0.109 

Ot 

per rad 


6.589 


6.245 

* 

a 

Limit of linearity of Cj^ , relative 

a 

to chord line, deg 

Figure 4.1.4 

3.6 

Table 4.1.4 

12.2 

a 

c( at (relative to chord line) , 

Figure 4.1.4 

16.0 

Table 4.1.1 

14.9 

Z 

max 

max 

deg 





base 

Reynolds number of airfoil wind- 
tunnel data 

Figure 4.1.4 

2.1 X 10® 

Table 4.1.1 

9 X 10® 


Maximum lift coefficient 

Figure 4.1'. 4 

1.59 

Table 4.1.1 

1.45 

max . 

base 







Correction of maximum c^^ to 

Figure 4.1.1 

0 

Figure 4.1.1 

-0.10 

max 

Nj^ = 2.3 X 1Q6 
e 






Maximum lift coefficient at wind- 


1.59 


1.35 

max 

tunnel test condition « 






(Cj^ ) + ACj^ 






max ; max 

base 




' 

c 

m 

o 

Zero-lift pitching-moment coefficient 

Figure 4.1.4 

-0.095 

Table 4.1.1 

0 

a. c. 

Location of aerodynamic center 



Table 4.1.1 

0.250 




Ay, % chord 


Figure 4.1.1: Effect of Reynolds number on section maximum lift coefficient 
(Reference 3) 
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% chord 


Figure 4.1, 



Variation of leading-edge sharpness parameter with airfoil 
thickness ratio (Reference 3) 
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t/c 




Figure 4.1.3; Variation of trailing-edge angle with airfoil thickness 
ratio (Reference 3) 
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4.2 Lift Characteristics of the Wing and Horizontal Tail 

The approximate wing lift-curve can be estimated when the following 
characteristics are known: 

1. lift curve slope, 

a 

2.. angle of attack for zero-lift, 

3. limit of linearity of the lift curve slope, a 

4. maximum lift coefficient, 

max 

5. angle of attack for maximum lift, a„ 

max 

In the following subsections methods' are presented for calculating 
these characteristics. 

4.2.1 Lift Curve Slope 

The lift curve slope of a tapered straight wing, in the subsonic 
region to M = 0.6, can be determined by the modified lifting line 
theory method of Polhamus . The lift curve slope is calculated as a 
function of the aspect ratio, A, the midchord sweep angle, 
ntimber, M, and the section lift curve slope, c. , by the following 

A/ 

a 

expression: 
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As compared to results from the lifting surface theory, Equation 
(4. 2.1.1) overestimates the value of by the value 


a \ /a 


(4. 2. 1.2) 


Pol 


where 


follows from Figure 4. 2. 1.1. 


4.2.2 Angle of Attack for Zero-Lift 


According to Reference 4, the zero-lift angle of attack of a wing 
may be calculated as follows ; 

Ao! 

ct = a + 0 (4. 2. 2.1) 

O O 0 

w 

where 

can be assumed to be the section zero-lift angle at low Mach'' 
number, obtained from Section 4.1. 

Aa^/0 represents the shift in the wing angle of attack for zero 
lift per degree of wing twist, obtained from Figure 4. 2. 2.1, 

0 is the twist of the wing tip with respect to the root section 
(negative for washout) . 

A Mach number correction is presented in Figure 4, 2.2. 2. This 
chart gives the ratio of the zero-lift angle of attack at any subsonic 
Mach nmnber to the corresponding value at M = 0.3. 

The upper limit of linearity of the wing lift curve slope is 
considered to be; 

A A 

a = a + 0 (4. 2. 2. 2) 

w 0 


45 



where 


a indicates the limit of linearity of the section airfoil lift 
curve slope, obtained from Section 4.1. 


4.2.3 Maximum Lift of the Wing 


The maximum lift coefficient and angle of attack for maximum 
lift of wings at subsonic conditions may be determined by the empirical 
method of Reference 4. For high-aspect-ratio (general aviation air- 
planes are concerned with high-aspect-ratio wings), untwisted, constant 
section wings: 


"L ~ c '^S, 
max Z max max 

max 


max , , . 

+ a + Aa 

-■t L- o g,. 

L L w L 

max a max 


(4. 2-3.1) 


(4. 2. 3. 2) 


The first term on the right side of Equation (4. 2. 3.1) is the maximum 
lift coefficient at M = 0.2 and the second term is the lift increment 
due to Mach effect. 

/c^ is obtained from Figure 4. 2. 3.1. 
max max 

c^ is the section maximum lift coefficient at low Mach number 
max 

obtained from Section 4.1. 

ACl is the Mach number correction obtained from Figure 4. 2. 3. 2. 
max 

is the wing lift curve slope obtained from Equation (4. 2. 1-2). 
a 

“q is the wing zero-lift angle obtained, for the appropriate 
w 

Mach number, from Equation (4. 2. 2.1). 
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&a„ is obtained from Figure 4. 2. 3. 3. The leading-edge 

max 

sharpness parameter. Ay, must be used in reading values from the 
charts. The value of Ay is expressed in percent chord and is obtained 
from Section 4.1. 

For twisted wings the calculations are more involved. The method 
used in this report is identical to the one described in Reference 8. 
The lift may be divided into additional and basic lift: 



(4. 2. 3. 3) 


The additional lift coefficient can be written as follows: 


a g g 


and the basic lift coefficient is; 


Aa 

® ^4 ^ I + “T > 

b a a 


(4. 2. 3. 4) 


(4. 2. 3. 5) 


where 

c/c is the ratio of the chord length at p = y/(b/2) to the mean 
§ 

geometric chord : 


c - 
g 


c 

r 


2 


1 + X 
'^r 2 


through follow from Figure 4. 2. 3. 4. 


(4. 2. 3. 6) 


f is the lift distribution function and can be obtained from 


Figure 4. 2. 3. 5. 

i 

e/0 is the ratio of the wing twist at p to the wing twist of the 
tip, with respect to the root section. 
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/ tan A ,r \ 

Ag = tan" ( g-- - - ' ) (4. 2. 3. 7) 

The maximum lift coefficient of the wing, , may be estimated from 

max 

the assumption that this coefficient is reached when the local section 

lift coefficient, c^, at any position along the span is equal to the 

local c^ for the corresponding section. This value may be found 
max 

by the process indicated in Figure 4. 2.3. 6. Spanwise variations of 

the local c^ and of the additional lift coefficient, c^ , for 
max a 

C, = 1 (Equation [4. 2. 3. 4]) and c„ (Equation [4. 2. 3. 5]), distributions 

are plotted. The spanwise variation of (c. - c, ) is plotted, and 

max b 

the minimum value of the ratio of (c. - c. ) to c. at C_ = 1 is 

tWm Xf Ij 

max b a 

then found. This ratio is considered to be the maximum lift coefficient 
of the wing. 

The angle of attack for maximum lift can be estimated with 
Equation (4. 2. 3. 2). 

4.2.4 Lift of the Wing and Horizontal Tail for the ATLIT Airplane 

Pertinent aspects of the calculation of the wing and horizontal 
tail of the ATLIT airplane at wind tunnel Mach conditions are summarized 
in Table 4 . 2 . 4 . 1 . 

During this study it appeared that the ATLIT wing had parabolic 

twist instead of the common linear twist. This led to a change in 

the value of Aa^/6. In Appendix A the calculations to obtain Aa^/8 

for the case of parabolic twist are shown. The calculations which 

lead to maximum lift coefficient, , of the wing are presented 

max 

in Appendix 's. 
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The results from Table 4. 2.4.1 are applied to the lift curves 
shown in Figure 4. 2.4.1. In this figure, results obtained with 
References 5 and 6 are also shown. The results obtained with Ref- 
erence 5 are in good' agreement with those obtained with Reference 6. 
However, the lift curve slope of the former is slightly less steep 
than the lift curve slope of the latter. This difference is caused 
by the thickness effect and discussed in more detail in Appendix C. 

The difference in the angle of attack for zero-lift between the results 
of References 5 and 6 and the results from Table 4. 2. 4.1 is under- 
standable. The former are lifting surface prediction methods, while 
the latter is based on the lifting line theory. 
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Table 4. 2, 4.1: 


Lift chacact eristics of wing and horizontal tail for 
the ATLIT airplane 


SysboL 

Desecipcion 

Raferenea 

Wi 

ng 


1 Horizontal Tail 

'Total 

Sxposec 

1 Reference 

Total 

1 Exposed 

M 

Hach ausber 

1 - mZ 

' 

0.081 

0.993 

0-081 

0,933 

- 

0.081 

0.993 

0.081 

0.993 

A 

Aspect ratio 

Table 2.1.1 

10.32 

9.61 

Table 2.1.1 

4.75 

4.28 


Twisc angle ac the clp wlcb re- 
spect CO rooc section, deg 

Table 2,1 

-3,0 

-3.0 

Table 2,1 

0 

0 

*c/2 

Sveep At half chord, deg 

Table 2.1.1 

0 

0 

' Table 2,1.1 

0 

0 

h. 

Leading edge sweep, deg 

Table 2.1,1 

3.67 

3.67 

Table 2.1.1 

0 

0 

Jr 

Leading-edgc-shacpness parasecer 

Table 4.1.2 

4.9 

4.9 

Table 4,1,2 

2.6 

2.6 

■=1 

a 

Section airfoil lift curve slope, 
rad'l 

Table 4.1.2 

6-589 

6.589 

Table 4,1,2 

6.245 

6.245 

tc 

/2i 

a 

- 

1.049 

1.049 

- 

0.994 

0.994 

“ Pol 

Wing lift curve slope according 
CO Polhaaus, rad*^ 

Eq. W.Z.1.1) 

5.401 

5,322 

Eq. (4.2.1.2) 

4,165 

3.993 


Correction factor 

Figure 4. 2. 1.1 

5.75 

5.95 

Figure 4,2. 1.1 

6.90 

6.85 

h. 

a 

Wing lift curve slope, rad"*^ 
deg”i 

Eq. C4. 2.1.2) 

5,090 

5-005 

Eq. (4.2. 1.2) 

3.873 

3.719 


“ 

O.OS38 

0.0374 

' 

0,0677 

0.0649 

sax 

Section airfoil oaxinum lift 
coefficient 

Table 4.1.2 

1.59 

1.59 

Table 4.1,2 

1.35 

1.35 


Haxlaus lift coefficient factor 




Figure 4, 2. 3.1 



sax aax 

at M-0,2 




0.9 

0.9 


Mach number correction to M-0.081 




Figure 4,2.3. 2' 



sax 




0 

0 

sax 

Haxlmua lift coefficient at 
H-0.081 

Appendix 3 

1.494 

1.494 

Eq. (4.2.3.1) 

1.215 

1.215 

“o 

Section cero-llft angle of 

Table 4,1.2 

-3.7 

-3.7 

Table 4.1,2 

0 

0 

ia /a 
0 

Shift in cero-llft angle of 
attack per unit twist 

Appendix A 

-0.27 

-0.27 

Figure 4.2.2, I 

0 

0 

“o 

V 

Wing rero-lift angle of attack, 
deg 

Eq. (4.2.2.1) 

-2.39 

-2.89 

Eq. (4.2. 2.1) 

0 

0 

o* 

Limit of linearity of section 
lift curve, deg 

Table 4,1.2 

3.6 

3.6 

Table 4.1.2 

12.2 

12.2 

< 

Limit of linearity of wing lift 
curve, deg 

Eq. (4, 2. 2. 2) 

4.41 

4.41 

Eq. (4. 2. 2. 2) 

12.2 

L2.2 

(C^ /Ct. ) 
sax a 

MaxiauB o for extended linearity 
conditions, deg 

- 

16.82 

17.09- 

- 

17-95 

L3.72 

sax 

Correction for flow separation, 
deg 

Figure ‘4. 2, 3. 3 

2.5 

2.5 

'Igure 4. 2. 3. 3 

0.3 ' t 

.S 

\ 

sax 

Angle of attack at at 

max 

Eq. (4. 2. 2. 2) 

16.43 

16.70 £ 

q. (4. 2. 3.2) 

18.75 I 

9.52 


M“0.031, relative to chord, deg 
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Figure 4. 2. 1.1;. Error' in Polhamus -formula' when compared with lifting 
surface method (Reference 10) 
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(a) X = 0 



T I j r I i 

0 10 20 30 40 50 60 


QUARTER-CHORD SWEEP ANGLE, (deg) 


Figure 4. 2. 2.1: Effect of linear twist on angle of attack for zero-lift 
(Reference 4) 
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(C) X = 1.0 



Figure 4, 2. 2.1: Concluded 





Figure 4. 2. 2. 2; flach number correction for zero-lift angle of attack 
for cambered airfoils (Reference 4) 
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0 


0 



Figure 4t2.3.2: Mach number correction for subsonic maximum lift of 

high— aspect-ratio, untwisted, constant airfoil section 
wings (Reference 3) 


55 














Figure 4, 2. 3. 3: Angle of attack increment for subsonic maximum Ixft 
of high-aspect-ratio, untwisted, constant airfoil 
section wings at M=0,2 - 0,6 (Reference 3) 


Figure 4. 2,3.4 ;• Coefficients for additional and basic lift distribution 
(Reference 8)’ 
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Figure 4. 2, 3. 5: Lift distribution function (Refe: 
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Figure 4. 
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2.3.6: Example of lift distribution (Reference 8) 
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Figure 4. 2. 4.1: Lift curve of wing and horizontal tail at N = 2.3x10 

■ ■ -fmr 05 * 

for the ATLIT airplane (total surface^v-rQ^ILt^^ -onO'V^ 
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4.3 Lift Due to Fuselage and Nacelles 


The method used in Reference 3 to predict the lift due to the 
fuselage and the nacelles is also used in this study. This method 
is also discussed in Reference 4. It is based on the asstimption that 
the flow is potential over the forward part of the body and has no 
viscous contribution in this region. At the position the potential 
flow over the forward portion of the body changes into a viscous flow 
for the remainder of the body. 

The expression for the lift coefficient of a body of revolution, 
based on the wing area, S^, is: 



2(k2-k)ag nc^ Jj 

eff eff c r 

( 57 . 3 ) 2 s J 

W X 


(4.3.1) 


where the first part on the right hand side of Equation (4.3.1) repre- 
sents the potential flow contribution ani the second part represents 
the viscous crossflow contribution and where 

(k 2 ~k^) is the apparent mass factor obtained from Figure 4.3.1 as 

a function of body fineness ratio, d /&„• 

max 5 

S is the body cross-sectional area at x . In this study the 
o o 

maximum cross-sectional area of the equivalent circular body will be 
used, which results in slightly optimistic contributions of bodies 
(see Reference 3) . 

x^ is the body station where the flow ceases to be potential. 

This is a function of x^ , the body station where the rate of change 
of the cross-sectional area with x reaches its maximum negative value, 
x^ and x^ are correlated in Figure 4.3.2. 
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°^Be£f effective angle of attack of the equivalent; circular 

body in degrees , or : 

“b „ * “b %, 

eff B 


a is the angle of attack of the actual body, identical to the 

D 

airplane angle of attack, a, using X-body axis as a reference. 

a is the zero-lift angle of the equivalent circular body rela- 

tive to the reference X-body axis. 

c^ is the steady-state cross-flow drag coefficient of a circular 
c 

cylinder of infinite length, obtained from Figure 4.3.3. 

q is the ratio of the drag of a finite cylinder to the drag of 
an infinite cylinder, obtained from Figure 4.3.4. 


S, is - the body length. I 

h '’MGINAr. PAGE is /o™® 

J" rdx represents half of the projected area of the equivalent 


circular body from x to the end of the body. Using the Simpson inte- 

O-f 


gration method: 


^B 

J r dx = 2 r Ax 


r is the body radius at any longitudinal station, 
dx, Ax is the increment length of the body. 

In many cases it will be possible to determine the location of 
Xj^ by inspection. For cases that are doubtful, the area distribution 
should be plotted and examined to determine the location where dS^/dx 
first reaches its maximum negative value. 
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Because the equation for lift of bodies is based on bodies of 
revolution, it is necessary to replace the actual body by an approxi- 
mate equivalent body of revolution to serve as a mathematical model 
for analysis. 

For the AIXIT Figure 4.3.5(a) shows the estimated equivalent 
.circular fuselage in relation to the actual fuselage. The equivalent 
fuselage has a zero-lift angle of 2.5' degrees below the reference 
X-body axis. The nacelle does not lend itself to such a simple esti- 
mate of equivalence because of its wide rectangular shape- and irregu- 
lar profile. However, in Figure 4.3.5(b) an approximation for 
equivalence is shown. The axis of the equivalent nacelle is parallel 
to the X-body axis. 

In Figure 4.3.5(b) the actual nacelle is shown.- Not included is 
the extension of 8 inches due to installation of the thrust/torque 
balance. This factor has not been taken into account in the calcu- 
lations. However, the effect of the extension on the lift of the 
nacelle will be negligible. 

Table 4.3.1 summarizes the calculations for the lift contributions 

of the fuselage and the nacelles. The lift coefficient of the fuselage 

based on the wing, S , is: 

w 

C == 0.002574(a^-2.5)+0.0000402(a^-2.5)2 (4.3.4) 

Ijf D 15 

and the total lift coefficient of the' nacelles based on the reference 
wing area is; 

2 

= 0.002031 + 0.0000201 (4.3.5) 

B B 

where 

is the angle of attack, in degrees , of the actual body and 
identical to the airplane angle of attack. 


62 



» K> 


Table 4.3.1: Fuselage and nacelle lift contribution 





Description 

Reference 

Fuselage 

Nacelle 

Mach number 

- 

0.081 

0.081 

Reference wing area, (ft^) 

Table 2.1.1 

14.40 (155.0) 

14.40 (155.0) 

Maximum diameter of equivalent 
circular body, m (in) 

Figure 4.3.5 

i.238‘ (48.75) 

0.889 (35.0) 

Maximum cross sectional area of 
equivalent circular body, (ft^) 

- 

1.204 (12.96) 

0.621 (6.68) 

Location of from nose of body, 

max 

Figure 4.3.5 

3.97 (156.3) 

0.76 (30.0) 

■Hi (in) 




Body length, m (in) 

Figure 4.3.5 

8.35 (328.8) 

2.43 (95.6) 

Fineness ratio 

- 

6.745 

2.731 

Reduced mass factor 

Figure 4.3.1 

0.882 

0.675 

Ratio of drag coefficient of finite 
to infinite length cylinder 

Figure 4.3.4 

0.650 

0.575 

Location from nose of S 

0 

max 

- 

0.475 

0.314 

Location from nose where potential 
flow ceases 

Figure 4.3.2 

0.627 

0.541 

® (in) 


5.24 (206.2) 

1.31 -(51.7) 

Half projected area of equivalent 

circular body from x to fft^) 

0 B 

Figure 4.3.5 

1.22 (13.1) 

0.34 (3.7) 

Zero-lift angle of equivalent circu- 
lar body, deg 

Figure 4.3.5 

-2.5 

0 

Angle of attack of equivalent 
circular body, deg 

Eq. (4.3.2) 

Ub-2.5 

°3 

M sin a 

®e£f 

- 

0.081 sin(a,-2.5) 

0.081 sin a 

15 

Crossflow drag coefficient of 
infinite length cylinder 

Figure 4,3.3 

1.2 

1.2 


Fuselage = 0.002574 Ca,-2- 5) +0.0000402 (a„-2,5)2 

D D 


Nacelles (2) C. = 0.002031 + 0.0000201 a_^ 

L b o 


63 







Me = M sin 




Figure 4.3.3: 


Steady state drag coefficient for 
infxnxte cylinder (Reference 3 ) 


Fxgure 4.3.4; 


Ratio of drag coefficient of a cylinder 
of finite length to that of a cylinder 
infinite length (Reference 3) 


L»I»j 


Equivalent circular fuselage 



Figure 4.3.5: Pertinent dimensions of approximated 

equivalent circular fuselage and nacelle 



Figure 4.3.5; Concluded 


4.4 Lift Due to Combined Wing-Fuselage-Macelle 


The addition of a body to a wing results in mutual interference 
effects. Lift of the wing-body combination is influenced by the body 
upwash effect on wing lift and the lift carry-over of wing panel 
loading -onto the body. Net wing upwash’ and downwash effects on the 
body influence body pitching moment primarily. Symmetrical body 
vortices which result from flow separation just behind or above the 
area of minimum pressure along the side of the body near the nose 
are normally negligible for most airplane types of wing-body config- 
urations . 


4.4.1 Wing-Fuselage-Nacelle Lift in the Linear Lift Range 

The lift of the wing-fuselage-nacelle combination accounting for 
the mutual interference effects of wing and fuselage may be estimated 
from:- 


w 

+ C., + [K + K.. xKC., ) a 

L „ L, L w(f) f(w)’’^ L V7 , S 

wfn f n a abs w 

w 

e 


(4.4. 1.1) 


where 

C, is the fuselage lift obtained from Section 4.3 
^f 

is the lift from the nacelles from Section 4.3 
n 

^w(f) ratio of the lift on the wing in the presence of 

the body to the lift on an isolated wing, obtained from Figure 4. 4. 1,1 

Kf(^) is the ratio of wing lift carry-over onto the body to wing 

lift alone, obtained from Figure 4. 4. 1-1 

(C^ ) is the lift-curve slope of the exposed wing panels,, obtained 
a 

w 


from Section 4.2 
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is the absolute angle of attack of the wing: 


w 


abs 


a • = a + X - a 
w , wo 

abs w 


(4. 4. 1.2) 


i^ is the incidence angle of the wing root with respect to the 


X-body axis 


is the sero-lift angle of attack of the wing, obtained from 
w 

Section 4.2 

is the exposed wing area; see Section 2.1. 
e 

Because of the lack of suitable data, the interference effects 

of the nacelles are not accounted for. 

The use of the interference factors, K and is restricted 

w(f) f(w) 

to wings which do not have sweptback trailing edges or sweptforward 
leading edges. The factors were obtained for wings mounted as midwings 
on bodies of revolution but have been used for other configurations. 


4.4.2 Maximum Lift of Wing-Fuselage-Nacelle Combination 


The addition of a body of revolution to a wing at high angles of 
attack increases the wing-induced angle of attack at all spanwise stations. 
The increase is greatest at the root and falls off in an exponential 
manner with increasing distance from the body. This effective increase ■ 
in angle of attack tends to make the wing in the presence of the body 
stall at a lower geometric angle of attack than that corresponding to 
the wing alone. However, this tendency to stall at a lower angle of 
attack may be modified by changes in the wing stall pattern. The 
changes are the result of nonlinear spanwise variation of body-induced 
flow and also of the partial coverage of the wing by the body. The 
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relative magnitudes of these latter two effects are largely dependent 
on specific wing planform shape. 

In absence of theoretical methods which predict the effect of 
the addition of a body to a wing on the maximum lift coefficient and 
corresponding angle of attack, Reference 4 presents empirical relations 
to predict this effect. The following, expressions are defined: 


and: 


where 


> 

max 


wf 


■(c^ ) 


max 


wf 


max 


w J 


(C. ) 


(4. 4.2.1) 


max 


w 


(a„ ) 

max 


wf L 


(“p ) 

max 


wf 


) 

max 


(ctp ) 
max 


(4. 4. 2. 2) 


w 


w 


(C^ ) ^ empirical correction factor, obtained 

max £ max 
wf w 

from Figure 4.4. 2. 2. 

(a ) / (a ) can be obtained from Figure 4. 4. 2. 3. Both 

max max 

wf w 

factors are presented as a function of the ratio of the fuselage diameter 
to the wing span, d/b, and the factor which follows from Figure 4. 4. 2.1. 


4.4.3 Wing-Fuselage-Nacelle Lift for the ATLIT 

/■ 

For the ATLIT the lift of the wing in the presence of the body 

and the carry-over of the wing lift onto the body, C , is calcu— 

w(f)+f(w) 

lated in Table 4. 4. 3.1. The net lift of the wing-fuselage-nacelle combi- 
nation in the linear lift range is: 
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(4. 4. 3.1) 


or; 


= C^+C^+C 
wfn f n w(f)+f(w) 


c: = 0.002574(0-2. 5) -f0.0000402(a-2. 5)2+ O.Q02031a + 0. 0000201 a2 + 0.0920(a+3. 4) (4.4, 3. 2) 

wfn 


Pertinent aspects of the calculations for 


> 

max 


and (a ) 
L 

max 


wf 


are listed in Table 4.4.3. 2. 

The net lift of the wing- fuselage-nacelle combination for the 
ATLIT airplane is summarized in Table 4. 4. 3. 3, while in Figure 4.4. 3.1 
the results are plotted and compared with full-scale wind tunnel data 
from Reference 2. The "fully clean" configuration of the ATLIT air- 
plane has not been tested with the horizontal tail removed. However, 
in Appendix D a procedure is shown with which the lift curve of the 
"fully clean" ATLIT, without horizontal tail, can be determined. 

Figure 4. 4. 3.1 also shows results obtained with References 5 and 6, 
respectively. 

The results obtained with Reference 5 do not include the lift 
due to the nacelles. However, the data obtained with Reference 6 
show that the lift contribution of the nacelles is substantial. Addition 
of the lift due to nacelles to the results of Reference 5 will result in 
a lift curve for the wing-fuselage-nacelle combination which shows good 
agreement with the experimental results in the. linear lift region. 

The results from Table 4.4^3. 3 show poor agreement with the full- 
scale wind tunnel data. The angle of attack for zero-lift is 1.5 degrees 
off. Part of the discrepancy is caused by the wing lift prediction as 
is shown in Figure 4. 2. 4.1. An additional factor is the omission of 
wing-nacelle interference effects. 
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The results obtained with Reference 6 show, in the linear lift 
region, good agreement with the experimental data. However, the 
predicted lift curve slope is too steep. This is caused by the 
thickness effect as is explained in Appendix C. 
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Table 4. 4. 3.1: Wing lift in linear range including mutual wing-fuselage 

interference 
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Table 4.4. 3. 2: Maximum lift of wing including mutual wing-fuselage 

interference 


Symbol 

X 

=2 

A, 

le 


(c,+l)A Can A, 
i- le 


d/b 


wf 


max 


^ ■■ 

max j- 
^ 

^ 

> max 

w 


> 


(a„ ) 

''L 

max 


(C, ) 


wfn 




) 

max 


wfn 


Description 

Reference 

Magnitu 

Wing caper ratio 

Table 2.1.1 

0.5 

Taper ratio correction factor 

Figure 4. 4. 2.1 

1.07 

Leading edge sweep of wing, deg 

Table 2.1.1 

3.67 

Wing aspect ratio 

Table 2.1.1 

10.32 

- 

- 

1.37 

- 

Table 4. 4. 3.1 

0.10 

! Ratio of of wing-fuselage to wing alone 

max 

Figure 4. 4, 2. 2 

1.0 

Ratio of stall angle of wing fuselage to wing 
alone 

Figure 4.4-2-3 

1.0 

Maximum' lift coefficient of wing alone 

Table 4. 2, 4.1 

1.494 

Stall angle of wing alone relaCaive to 
zero-lift, deg 

1 

; Table 4. 2. 4.1 

16.4 + 

Maximuni lift coefficient of wing— fuselage*' 
nacelle combination 

Eq. (4.4. 2.1) 

1.494 

stall angle of wing-fuselage-nacelle 
combination relative to zero lift, deg 

Eq. (4. 4. 2. 2) 

19.3 
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Table 4. 4. 3, 3: Summary of wing-fuselage-naeelle lift of the ATLIT 

airplane 


a, deg 

. deg 

abs 

Eq. (4. 4. 1.2) 

% 

Sq. C4.3.4) 

n 

Eq. (4.3.5) 

% 

w(f) + f (w) 
Table (4. 4. 3.1) 

\fn 

Eq. (4.4. 3.1) 

-4 

-0.6 

-0.01673 + 0.00170 

-0.00812 + 0.00032 

-0.05520 

-0.07803 

-2 

1.4 

-0.01158 + 0.00081 

-0.00406 + 0.00008 

0.12880 

0.11405 

0 

3.4 

-0.00644 + 0.00025 

0 

0.31280 

0.30661 

2 

5.4 

-0.00129 + 0.00001 

0.00406 + 0.00008 

0.49680 

0.49966 

4 

7.4 

0.00386 + 0.00009 

0.00812 + 0.00032 

0.68080 

0.69319 

15.9 

*19.3 

0.03449 + 0.00722 

0.03229 + 0.00508 

*1.49400 

1.57308 


i: 

from Table 4*4.3. 2 
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2.0 



Figure 4. 4. 1.1: 


Lift ratios K 


w(f) 


and K 


f(w) 


based on slender-body 


theory (Reference 3) 
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1.5 



Figure 4.4. 2.1: Taper ratio correction, factor (Reference 3) 
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Figure 4. 4. 3.1: Comparison. of lift predictions with wind tunnel data 

^6 

horizontal tail and propellers removed) 


2,3x10' 
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4.5 Wing Zero-Lift Pitching-Moment and Aerodynamic Center of Win 


and Horizontal Tail 


Subsonic zero-lift pitching moment coefficients for constant 
section airfoil lifting surfaces can be predicted with the following 
expression: 


A cos^ A / , 

C = c 

m m A + 2 cosA , , 
o o c/4 


where 


c is the section zero-lift pitching moment obtained from 
o 

Section 4.1 

AC /e is the change in wing zero-lift pitching moment coefficient 
o 

due to a unit change in linear wing twist. This parameter is obtained 
from Figure 4.5.1. 

9 is the twist of the wing tip with respect to the root section, 
in degrees (negative for washout) 

The aerodynamic center, the point about which the lifting surface 
pitching-moment coefficient is invariant with lift, may be determined 
relative to a desired reference center on and as ratio of the mean 
aerodynamic chord of the lifting surface by using Figure 4.5.2 and 
the expression: 


X X \ c 
n ac \ r 


where 


r / c 


-dC /dC_ is the static margin, the distance from the reference 

m Li 

center on the mean aerodynamic chord of the lifting surface to the 
aerodynamic center of the surface as a ratio of the mean aerodynamic 


chord 
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x^/c^ is the chordwise distance from the wing apex to the point 
about which the pitching moment is desired, measured in root chords, 
positive aft or: 



y- tan 
c le 


(4.5.3) 


when the leading edge of the mean aerodynamic chord is the moment 


center . 


X / c is the chordwise distance from the wing apex to the aero- 
dynamic center, measured in root chords, positive aft, obtained from 
Figure 4.5.2. 

labels 4.5.1 and 4.5.2 summarize the calculations made to deter- 


mine the zero-lift pitching moment of the wing and the horizontal tail 
of the ATLIT and the location of the aerodynamic centers of the surfaces 
relative to the leading edges of the mean aerodynamic chords of the 
surfaces . 


82 



Table 4.5.1: Zero~lift pitching moment coefficient of wing and horizontal 
tail 


Symbol 

Description 

Reference 

Wing 

Horizontal Tail 

Total 

Exposed 

Total 

Exposed 

A 

Aspect ratio 

Table 2.1.1 

10.32 

9.61 . 

4.75 

4.28 

8 

Twist angle at wing tip, deg 

Table 2.1 

-3.0 

-3.0 

0 

0 

^c/4 

Sweep of quarter chord line, dag 

Table 2.1.1 

1.835 

1.835 

0 

0, 

X 

Taper ratio 

Table 2.1.1 

0.50 

0,527 

1 

1 

AC /0 
m 
o 

Shift in per degree of wing twist 

o 

Figure 4.5.1' 

-0.0004 

-0.0004 

0 

0 

c 

Section airfoil zero-lift pitching 

Table- 4.1. 2 

-0.095 

-.095 

0 

0 


moment coefficient 






C 

m 

o 

Zero-lift pitching moment coefficient 

Eq. (4.5.1) 

-0,0783 

-0.0774 

0 

0 


Table 4.5.2; Aerodynamic center of wing and horizontal tail of the ATLIT 


Symbol 

Description 

Reference 

Wing 

Horizontal Tail 

Total 

Exposed 

Total 

Exposed 

M 

Mach number 

- 

0.081 

0.081 

0.081 

0.081 

8 

Vl - 

- 

0.9967 

0.9967 

0.9967 

0.9967 


Root chord of surface > m (in) 

Table 2.1.1 

1.575 

(62.00) 

1.495 ■ 
(58.87) 

.0.871 

(34.29) 

0.871 

(34.29) 

c 

Mean aerodynamic chord, m (In) 

Table 2.1,1 

1.225 

(48.22) 

1.178 

(46.38) 

0.871 

(34.29) 

0.871 

(34.29) 


Lateral position of c from root 
chord, m (in) 

Table 2.1,1 

2.709 

(106.67) 

2.460 

(96.85) 

1.033 

(40.68) 

1.033 

(40.68) 

le 

Sweep of leading edge» deg 

Table 2.1.1 

3.67 

3.67 

0 

0 

A tan A, 
le 

' 

- 

0,6619 

0.6164 

0 

0 

tan A^^/8 

* 

- 

- 

0.0644 

0.0644 

0 

0 

X /c 
ac r 

Distance from apex of surface 
to a.c as ratio of c 

r 

Figure 4.5.2 

0.294 

0.297 

0.250' 

0.250 

X /c 
n r 

Distance from apex of surface to 
leading edge of m.a.c. as ratio 
of c 

r 

Eq. (4.5.3) 

■0.110 

0.106 

0 

0 

a« c. 

Aerodynamic center relative to 
leading edge of mean aerodynamic 
chord as ratio of c 

Eq. (4-5.2) 

0.236 

0.243 

0.250 

0.250 
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(a) X = 0 



QUARTER-CHORD SWEEP ANGLE, Ac/4 (deg) 


(b) X = 0.5 



QUARTER-CHORD SWEEP ANGLE, Ac/4 (deg) 


Figure 4,5.1. Effect of linear twist on the zero— lift pitching 
moment of the lifting surface (Reference 4) 
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1.6 



P tan Ajg 


(e) A =0.5. 


'ac 


2.4 

2.0 

1.6 

1.2 

.8 

.4 


A tan Ajg 


6 






5 4 


3 




2 

1 

— 


— - 




tan A|e 


tan A 


le 


(f) A=l. 


Figure 4.5.2: Concluded 
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4.6 Wing-Fuselage Pitching Moment at Zero-Lift 


The addition of a fuselage to a wing results in a fuselage contri- 
bution, (C^ )^, to the pitching moment at zero lift. This contribution 
o 

may be estimated from Figure 4.6.1, which is based on streamline bodies 

of circular or near circular cross section for midwing conditions. For 

high- or low-wing configurations a- positive or negative increment, 

(AC )^, of 0.004 has to be added, respectively, to the value obtained 
o 

from Figure 4.6.1. In the absence of suitable data, the effects of 

the nacelles on C are considered to be negligible, 

m 

o 

The effect of Mach number on the wing-fuselage zero-lift pitching 
moment is presented in Figure 4.6.2. When using this chart, no correc- 
tion should be made to the section c value and the wing pitching 

• m 

o 

moment, (C ) , at zero-lift. The wing-fuselage-nacelle pitching- 
m- w 
o 

moment coefficient is : 




(C ) ^ 
m wfn 
o 


(C^ + (C^ ). + (AC ) 

m w m r m r 

O O 0-1 


o 


(4.6.1) 


where 

(C ) is the wing zero-lift pitching moment coefficient uncorrected 
o 

for Mach number effects, obtained from Section 4.5. 


(C is the fuselage zero-lift pitching moment, obtained from 
o 


Figure 4.6.1. 


(AC is the correction for high- or low-wing configurations, 
m t 
o 

(C ),./ (C is the Mach number correction factor obtained 

m M m M=0 
o o 

from Figure 4.6.2. ' 
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In Table 4.6.1 the calculations are summarized which lead to the 
wing-fuselage-nacelle pitching moment coefficient at zero-lift of the 
ATLIT airplane. The final result is: 



'wfn 


o 


-0.1072 


(4.6.2) 
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Table 4.6.1: Wing-fuselage-nacelle pitching moment coefficient at zero 
lift of the ATLIT airplane 


Symbol 

Description 

Reference 

Magnitude 

w 

Width o£ fuselage at the wing, m (ft) 

Figure 2.1.1 

1.219 (4.0) 

Sf 

Planform area of fuselage, (ft^) 

Figure 4.6.3 

7.72 (83.1) 


Planform area of fuselage forward 
of c/4 of wing, (ft^) 

Figure 4.6.3 

3.17 (34.1) 


Length of fuselage, m (ft) 

Figure 4.6.3 

8.352 (27.40) 


Distance from nose of fuselage to 
c/4 of wing, m (ft) 

Figure 4.6,3 

3.127 (10.26) 

S 

V 

Reference wing area, (ft^) 

Table 2.1.1 

14.4 (155.0) 

c 

w 

Wing mean- aerodynamic chord 

Table 2.1.1 

1.225 (4.013) 


■ Incidence of zero-lift line of 
wing = > rad 

w 

Table 4. 4. 3.1 

0.0593 

- w2/S^ , 

- 

- 

0.193 


- 

- 

- 0.154 

o 


Figure 4.6.1 

-0.115 

0 

C of fuselage 

m 

O 

- 

-0,0249 

)f 
m r 

o 

Correction for low-wing configuration 
of airplane 

- 

-0.004' 

0 

Zero— lift pitching moment coefficient 
of wing 

.Table 4.5.1 

-0.0783 

(c ) . 

IS wfn 
Q 

Zero-lift pitching moment of wing- 
fuselage-nacelle combination 

Eq. (4.6.1) 

-0.1072 
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Figure 4.6.1; Effect of fuselage on G 
(Reference 3) 


Midwing configuration 
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Figure 4.6.2: Mach effect on zero— lift pitching 

moment coefficient (Reference 4) 
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c/4 of ulng 





Figure 4.6.3; Top view of the ATLIT airplane 


4.7 Fuselage and Nacelle Pitching Moments 


The slope of the pitching moment curve of the fuselage and the 
nacelles at subsonic Mach numbers may be determined from the following 


expression; 

a 


where 

(C^ )g is the slope, of the body pitching moment curve, based on 


j r(x^-x)dx (4.7.1) 




4ri c 


>B = - 
a 57 


‘2 1 ^ J 

.3 S c J 


w w 


(x -x)dS + 

"" (57. 3)^3 


a 

the reference wing area, about a moment center at a longitudinal dis- 
tance, X , from the nose of the body. 

Above expression is based on potential-flow lift effects on the 
forebody and on viscous-flow lift effects on the afterbody, which are 
discussed in more detail in Section 4.3. 

Table 4.7.1 summarizes the calculation of the slope of the pitching 
moment curve of the fuselage and the nacelles of the ATLIT airplane 
about the leading edge of the total wing mean aerodynamic chord. Table 
4.7.2 shows the tabular integration procedure used to obtain the values 

X 

o B 

of the integrals, / (x^-x)dS^ and f r(x -x)dx, for the fuselage. The 

o . X 

o 

same procedure is used for the nacelles as is shown in Table 4.7.3. 

The slope of the pitching-moment curve of the fuselage and the 
nacelles of the ATLIT about the leading edge of the total wing mean 
aerodynamic chord is: 


(C^ 5fn " 0-00533 - 0.000252a , (4.7.2) 

a 

where 

a is the angle of attack with respect to the X-body axis in 
degrees . 
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Not included in the calculation of the pitching moment curve 
slope of the nacelles is .the extension of the nacelles hy- 8- Inches 
due to installation of the thrust/torque balances. However, the 
effect of the extension on the pitching moment of the nacelles is 
negligible. 
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Table 4.7.1: 


Fuselage and nacelle pitching moments 
of the ATLIT airplane 


Symbol 

Description 

Reference 

Fuselage 

Nacelle 

M 

Mach number 

- 

0.081 

0.081 

s 

w 

Reference wing area, (ft^) 

Table 2.1.1 

14.4 (155.0) 

14.4 (155.0) 

s 

X 

Equivalent body cross-section, (ft^X 

Figure 4.3.5 

Variable 

Variable 

s 

Body length, m (in) 

Table 4.3.1 

8.351 (323.8) 

2.428 (95.6) 


Location from body nose where 
potential flow ceases, m (in) 

Table 4.3.1 

5.237 (206.2) 

1.313 (51.7) 


Distance from body nose to leading 
edge of wing mean aerodynamic chord, 
m (in) 

Figure 4.6.3 

2,819 (111.0) 

1.334 (52.5) 

X 

Distance from body nose to centroid 
of d quantity,, m (in) 

Figure 4.3.5 

Variable 

Variable 

r. 

Effective body radius of Ax segment, 
m (in) 

Figure 4.3.5 

Variable 

Variable 

c 

w 

Wing mean aerodynamic chord, m (ft) 

Table 2.1.1 

1.225 (4.018) 

1.225 (4.018) 

^2-^1 

Reduced mass factor 

Table 4.3.1 

0.882 

0.675 

n 

Ratio of drag coefficient of finite 
to infinite length cylinder 

Table 4.3.1 

0*650 

0.575 

c 

Cross flow drag coefficient of 
Infinite length cylinder 

Table 4.3.1 

1.2 

1,2 

X 

o 

S (x ^x)dS 

m X 

o 

X 

o 

5 (s^-x)AS^, m3 (ft 3) 
o 

Table 4. 7. 2/3 

2.10 (74.2) 

0.62 (22.0) 

S' 

/ r(x^-x)dx 

X 

0 

2 r(x -x)Ax, m3 (ft3) 
01 
X 

o 

Table 4. 7. 2/3 

-4.38 (-154,7) 

-0.17 (-6.0) 

a 

Slope of fuselage pitching moment 
curve, deg~3 

Eq. (4.7.1) 

0. 00367-0. 000236O 


(C. )'n 

a 

Slope of nacelle (one) pitching 
moment curve, deg~^ 

Eq. (4.7.1) 


0.00083-0. 0000081a 

Surnmary ; (C 

), = 0.00533 - 0.000252a 

a 
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Table 4.7.2: Tabular integration of fuselage pitching 

moment parameters 



Bistance from 


nose CO S , 
x’ 


in 

X, in 

0 

10.9 

16.3 


! 

26.5 

32.5 i 

i 

41.0 

49.4 

j 

56.6 

63.3 

71.6 

79.4 

87.5 

95.6 

103.5 

111.3 

119.1 

■ 126.9 

134.7 

142.5 

150.0 

157.5 

165.7 

173.8 

181.6 

189.4 

197.8 

206.2 



/ (.x-x) d s 
o “ ^ 



w, in 

S , 1q 2 
x' 

4S , 

X 

(x -x)AS 
in3 

0 

0 





209 

20921 

16.3 

209 





334 

28223 

26.3 

543 





354 

24780 

33.3 

897 





285 1 

15504 

38.3 

1182 





197 

7762 

41.9 

1379 





J.69 

3972 

44.4 

1548 





136 

1020 

46.3 

1684 





88 

-713 

47.5 

' 1772 





98 

-2323 

48.8 

1870 





0 

0 

48.8 

1870 





0 

0 

48.8 

1870 





-124 

8754 

47.5 

1746 





-239 

20745 ■ 

43.3 

1507 




• 


128645 


(x -rx) is = 128645 iti^ = 74.2 

01 1C 
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Table h.l .2'. Concluded 






.. . 

Distance 
from nose. 


■I 


r(x -x)Ax, 

in. 

Ax, in. 

B 

X, in. 

ln3 

206,2 

16.3 

21.6 

214.4 

-34603 

222.5 

16.3 

19,7 

230.7 

-38434 

238.8 

14.9 

17.5 

246.3 

-35280 

253,7 

16.3 

15.6 

261.9 

-38376 

270.0 

15.6 

13.8 

277.8 

-35914 

285.6 

15.6 

11.6 

293.4 

-33005 

301.2 

16.3 

9.7 

309.4 

-31363 

317.5 

11.2 

7.3 

323.1 

-18533 

328.7 









-267308 


^f 




/ r(x -x)dx = 2 r(x -x)Ax - -267308 in^ 

m Q 

O ‘ o 

= -154.7 ft^ 
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Table 4.7.3: Tabular integration of nacelle pitching 

moment parameters 








4-8 Wing-Fuselage-Nacelle Pitching Moments 

The wlng-fuselage-nacelle pitching moment characteristics are 
considered in terms- of pitching moment slopes , aerodynamic center and 
pitching moment coefficient. A first-order approximation of the 
pitching moment coefficient beyond the limit of linearity of the lift 
curve slope up to the stall is also considered. The prediction method 
described in this section is similar to- the one in Reference 3. 

4.8.1 Factors Contributing to Wing-Fuselage-Nacelle Pitching Moments 

a) Zero-lift pitching moments; The zero-lift pitching moments of the 
wing, fuselage, nacelles and wing-fusel-age interference were accounted 
for in Section 4.6. For the ATLIT: 

’wfn - (4.6.2) 

o 

b) Fuselage and nacelle pitching moments;. The. fuselage and nacelle 

pitching moments due to potential- and viscous-flow lift effects were 

calculated in Section 4.7. For the ATLIT, with the moment center about 

the leading edge of the total wing mean aerodynamic (geometric) chord: 

(G^ ^fn ~ 0*00533 - 0.000252a per degree ^(4.7.2) 

a 

c) Wing pitching moment due to effective wing lift; This coefficient 
includes the effects of body upwash on the wing and wing lift carry- 
over onto the fuselage. The following expression calculates the pitching 
moment slope about the leading edge of the total wing mean aerodynamic 
(geometric) chord; 



(4. 8. 1.1) ■ 
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where 


(x /c ) is the aerodynaroic center of the wing in the presence 

SC r w( i ) 
e 

of the fuselage as a fraction of and about the leading edge of the root 

chord of the exposed firing panels, obtained from Figure 4. 5. ’2. 

^w(f) ratio of the lift on the wing in the presence of the 

body to the lift on the isolated ^'rLng, obtained from Figure 4. 4. 1.1. 

'K-.^ V is the ratio of wing lift carry-over onto the body to wing 
r {.w; 

lift alone, obtained from Figure 4-. 4. 1.1. 

(x^^/c^ ^f(w) contribution to the aerodynamic center due 

to the lift carry-over of the wing on the fuselage. For general aviation 
airplanes (A^ /l-M^‘ ^ 4 , d/b < 0.5) this contribution is: 



1 b-d 
4 2c 

r 


tanA if 
c/4 




(l-k)+ 

Jlk k 


(l-k)2 

-f- (1-k) 


k 


(4. 8. 1.2) 


The wing pitching moment slope, about the leading edge of the mean 
aerodynamic chord due to the effective Tjing lift in presence of the 
body for the ATLIT' airplane is summarized in Table 4. 8. 1.1: 

^^m ^w(f)+f(w) “ -0*02464 per degree • (4. 8. 1.3) 

d) Wing pitching moment due to wing drag: The wing pitching moment 

due to wing drag can be predicted as follows : 


^^m ^w(D) 
a 


= - C. 


wfn 


^ 

«L>. 

a 


w 


>w 

a 

z 

w 

wfn 


ttA e 
w 

c 

w 


where 


(4. 8. 1.4) 


is the vertical distance from the X-body axis to the quarter 
chord of the mean aerodynamic chord of the wing, positive down, 
e is the wing efficiency factor for induced drag. 
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For the ATLIT the distance, z^, is negligible. In that case the 
wing pitching moment slope due to wing drag is zero : 


^^m ^w(D) ^ 

a 


(4. 8, 1.5) 


e) Fuselage and nacelle "free moments": The fuselage and nacelle "free 

moments" due to induced flow from the wing can be estimated with the 
Multhopp method. It is indicated that, in considering wing lift carry- 
over onto the body, there remains a free moment of the body due to 


wing upwash ahead of the wing and downwash behind the wing or: 


^^m ^B(e) 
ct 


2 


36.5 S c 
■w w 


2 


where 

w is the mean width of the body planform segment Ax 

D 

33/ 3a is the variation of the local flow with the angle of attack, a. 
Curves of 33 / 3a are shown in Figure 4. 8. 1.1 as a function of the Ax seg- 
ment position ahead of the wing leading edge, x^^/c^, where c^ is the 
exposed root chord of the wing for the fuselage, and the -chord at the 
centerline of the nacelle for the nacelle. For Ax segments immediately 
ahead of the wing leading edge, 33' / 9a rises so abruptly that integrated 
values, 33/ 9a, are given based on the length, x^, of the segment 
adjacent to the wing leading edge. For segments aft of the trailing 


edge of the wing, 33/3d is assumed to vary linearly, or; 


93 _ _1 / -[ w 

9a \ 3a 


where 

9e /3a can be obtained from Section 4.9 and is considered to be 

W 

similar to (3e^/3a)j^ 

is the distance from the wing trailing edge to the centroid of 
the last aft Ax segment 
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is the distance from the wing trailing edge to the centroid 
of the Ax segment. 

Fuselage and nacelle "free moments" for the ATLIT are stimmarized 
in Table 4. 8. 1.2 or; 

^^m^^f(s) +n(e) " ^“^m^^fCe) 

= 0.00966 + 0.00737 (4. 8.1.8) 

= 0.01703 per degree 


4.8.2 Static Margin of Wing-Fuselage-Nacelle Combination 


The pitching-moment characteristics in terms of static margin, 
the distance from the center of gravity to the aerodynamic center, 
are obtained from the following expression: 


^d C. 


m. 


tdC 


h 'eg 


eg 


w 


>le 

a 

2 C, 


(4. 8. 2.1) 


where 

Xcg/ c^' is the distance from the leading edge of the total wing 
mean aerodynamic chord to the center of gravity as a ratio of the mean 
aerodjmamic chord. 

(2C^ )j^g is the pitching moment slope about the leading edge of 
a 

the mean aerodynamic chord 

2Cj^ is the lift curve slope of the wing-fuselage-nacelle combi- 
• a 

nation, which can be obtained from Section 4:4. 

For the wing-fuselage-nacelle configuration: 


^ ^ -i-f (w) 

^wfn 
a 


^ (s) +n(e)j 


le 


(4. 8.2.2) 


104 



where 


is the pitching moment slope of the fuselage and nacelle 
a 

obtained from Section 4.7 

^wff ) + f fw") wing pitching moment slope due to effective 

wing lift, obtained from Equation (4. 8. 1.1) 

(C ) , , is the pitching moment slope due to wing drag , obtained 

a ^ ^ 

from Equation (4. 8. 1.4) 

(C pitching moment due to the fuselage and nacelle 

"free moments" 

^wfn lift curve slope of the wing-fuselage-nacelle 

a 

combination, obtained from Section 4.4 

To express the static margin as a function of the lift coefficient 
of the wing-fuselage-nacelle configuration, replace a in Equation 
(4. 8. 2. 2) by: 


wfn 


\<rfn 

a 


+ (a ) ^ 

0 wfn 


(4. 8. 2. 3) 


where 


represents the angle of attack at zero-lift of the wing- 
fuselage-nacelle combination, which can be obtained from Section 4.4. 
For the ATLIT airnlane: 


d C 


m 


d C. 


= - 0.25 - 


(0.00533 - 0.Q00252a) - 0.02464 + 0. + 0.01703 


L^0.25c 


0.0964 


= -0.2263 + 0.00261a 


Substitution of the following expression: 


(4. 8. 2.4) 


_ wfn 
“ 0.0964 

into Equation (4. 8. 2. 4) results in: 
/d C 


- 3.2 


ra 


.d C. 


= -0.2347 + 0.02707 C, 


0.25c 


wfn 


(4. 8.2. 5) 


(4.8. 2.6) 
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4.8.3 Pitching Moment Coefficient of Wing-Fuselage-Nacelle Combination 


Linear lift range-: The pitching moment coefficient of the air- 
plane without horizontal tail in the linear lift range can be deter- 
mined as follows : 


(Q 


m 


wfn 


0.25c J 


^dC 


m 


Kd C. 


L '0.25c 




wfn 


+ (C ) . 
m wfn 


(4. 8. 3.1) 


where 

- (d C /dC^)rt follows from Equation' (4. 8. 2.5) and (C ) , 
m L 0.25c ^ m 'wfn 

o 

from Equation (4.6.2). 

For the ATLIT airplane the pitching moment of the airplane without 
horizontal tail about quarter chord of the mean aerodynamic chord is 
given by : 


(C, 


)_ ^ 0.2347 C^ - 0.01354 C, ^ 

m 0.25c L , L ^ 

wfn wfn wfn 


- 0.1072 


(4. 8. 3. 2) 


Non-linear "lift range: In Reference 3 a method is given which 

estimates the pitching moment coefficient in the lift region between 
the upper limit of linearity of the lift curve slope and stall. 

The average pitching moment slope in the non-linear lift region 
can be obtained as follows: 

a) Calculate the average value of the lift curve slope of the wing 
in the non-linear range approaching stall by: 

(C, 


(<1 ’ w , 
a s 


1 

2 




-T ) - (C_ ) {a - (a ) } 

L w L w w 0 w 
max a 


a 


(a. 


* 

) “ “ 
w w 


(4. 8. 3. 3) 


max 


For the ATLIT, the wing lift data can be obtained from Section 4.2. 
The average wing lift slope is ; 

(C^ ) = 0.0796 per degree (4. 8.3.4) 

Li W 

a s 
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b) Calculate the difference in linear and non-linear lift curve slopes 
from: 

'“e ’w - >» - 'h, (4.S.3.5) 

as a as 

For the ATLIT airplane: 


(AC^ ) = 0.0888 - 0.0796 

L w 


(4. 8, 3. 6) 


= 0.0092 per degree 

c) Calculate the average slope of the wing pitching moment coefficient 
about the leading edge of the \^rf.ng mean aerodynamic chord in the non- 
linear lift range to stall by: 


m w 
a s 


- (ac) (C, 

w L w 
s as 


(4. 8, 3. 7) 


where 

(ac) is the average value of the aerodynamic center in the non- 
s 

linear range of the wing lift curve slope to stall expressed as a 
fraction of the wing mean aerodynamic chord, assumed to be 0.375. 

For the ATLIT: 


(C ) = -0.375 (0.0796) 

m w ' 

a s 


(4. 8. 3. 8) 


= - 0.0299 per degree 

d) Calculate the average value of the angle of attack in the non-linear 
lift range from: 

a^ = -I- (a* + ) (4. 8.3. 9) 

^ax 


For the ATLIT the average value of a in the non-linear range can be 


determined from Figure 4. 4. 3.1: 


cis = f + 15.9) 
- 10 degrees 


(4.8.3.10) 
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e) Calculate the pitching moment slope in the- non-linear lift range 
as follows: 


'dC 


m 


.dC, 


eg. 


(C ), + (C ) + (G . 

: m fn m w m “^fCe) +n(e) 

eg a g s a 

(C^ ) ^r - (AC^ ) 
c - L wfn L w 

a g s 


(4.8.3.11) 


where 


(C )_ follows from Sectxon 4.7 
m fn 
a 

(C ) follows from Equation (4. 8. 3. 7) 
m w ■ 

g s 

CC )r/ \ ^ is obtained from Equation (4. 8; 1.8). 

m r ( e ) + n( e) 
g ' 

For the ATLIT airplane the pitching moment slope is the non-linear 


lift range is: 



J = -0.25 
0.25c 

s 


(0.00533 - 0.000252g ) - 0.0299 + 0.01703’ 

; S 

0.0964 - 0.0092 


-0.1346 


(4.8.3.12) 


4. 8.-4 - Pitching Moment Characteristics- of the ATLIT 

The pitching moment characteristics of the ATLIT airplane, including 
the non-linear region, are summarized in Table 4. 8. 4.1, while in Figure 
4. 8.4.1 the results are shown and compared with the full-scale wind 
tunnel data. These results are obtained from Appendix D, because no 
wind tunnel data were obtained with the ATLIT in the "fully clean" 
configuration and horizontal tail removed. Figure 4- 8.4.1 also shows 
results, obtained with References 5. and 6. 

The pitching moment curve obtained with Reference 5 shows poor 
agreement with the experimental results. The results obtained with 
Reference 5 do not include the pitching moment due to the nacelles. 
However, the data obtained with Reference 6 show that the contribution 
of the nacelles is substantial. 
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The results obtained with Reference 6 show good agreement .x^ith 
the experimental data. This- is also true for the results from Table 
•4. 8. 4.1. Both methods predict the pitching moment slope, fairly well. 

In Section 4.13 the stabilizer effectiveness in lift and pitch, 

C and C" , respectively, is discussed. For the ATLIT the stabilizer 

Xj . Ill"; ’ 

^ ^h 

effectiveness in pitch is C -0.08 per degree of stabilizer deflec- 

^h ■ 

tion. This means that the discrepancy betx^een the predicted pitching 
moment curve (from Table 4. 8. 4.1) and the experimental curve is identical 
to a stabilizer deflection, i^^, of approximately one degree. 

Note; 


At the time of finishing this report it appeared that the pitching 
moment coefficient, obtained with Reference 6, is defined as; 


C _ Pitching Moment 


m 


S b /2 

“ w w 


(4,8.4.1) 


while the normal definition is: 


Q _ Pitching Moment 

m “ ri “ 
q S c 
w w 


(4. 8.4. 2) 


The pitching moment coefficients obtained with Reference 6 and shown 
in Figure 4. 8. 4.1 have tO'be corrected in the follox<;ing manner: 


b 11 

C = c ~^=4.98C (4. 8. 4.3) 

” Vf.6 c “Ref. 6- 

w 

This correction will result in poor agreement of the results obtained 
with Reference 6 with the X7ind tunnel data of Reference 2. 
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Table 4, 8. 1.1: Wing pitching moment due to wing lift including 
mutual wing-fuselage interference 


Symbol 

DescripCion 

Reference 

Ifagnitude 


Total wing mean aerodynamic' chord, m (in) 

Table 2-1.1 

1.221. (43.22) 

' c 

r 

e 

Root chord of exposed wing panels, m (in)' 

Table 2.1.1 

1.495 (58.87) 

S 

w 

' Reference wing area, (ft^) 

Table 2.1,1 

14.40 (155.0) 

S 

V 

e 

Exposed panel wing area, (ft^) 

Table 2.1.1 

12,53 (134,8) 

a e 

Lift curve slope of exposed wing panels, deg~^ 

table 4. 2. 4.1 

0.0874 

‘'w(f) 

Ratio of lift of wing in presence of fuselage 
to wing alone, 

' Table 4. 4. 3,1 

1.08 


Ratio of wing lift carry-over onto fuselage 
to wing alone 

Table 4. 4. 3.1 

0.13 


Aerodynamic center of wing in presence of fuse- 
lage, as fraction of and about leading edge of 

*^r 

Figure 4.5.2 

0,297 

d 

Fuselage width at wing, m (ia) 

Table 4.6.1 

1.219 (43.0) 

b 

Wing span, m (in) 

Table 2.1.1 

12.19 (480.0) 

k 

d/b 

- 

0,10 

^c/4 

Sweep of wing quarter chord line, deg 

Table 2.1.1 

1.835 


Contribution to the aerodynasaic center due to 
lift carry-over of wing onto fuselage, as 
fraction of c 

r 

Eq. (4.-8. 1.2) 

0.268 

ys 

Lateral distance from root chord to total wing 
mean aerodynamic chord, m (in) 

Table 2.1.1 

2.709 (106.67) 

A, 

le 

Leading edge sweep angle of wing, deg 

Table 2.1.1 

3.67 

+ f (w) 

Pitching moment slope, about leading edge 
of wing mean aerodynamic chord, due to effective 
wing lift, deg“l 

Eq. (4. 8.1.1) 

-0.02464 


no 









Table 4. 8. 1,2: "Free moments" of fuselage and nacelle 


* 

KEPRODUCIBILrrY OF THE 
ORIGINAL PAGE IS POOR 


Symbol 

Description 

Reference 

"Magnitude 

"b 

Mean, width of body planform segment ix, m (in) 

Figure 4. 8. 1.2 

Variable 


Distance from wing leading edge Co centroid 
of forward Ax segments, m (in) 

Figure 4. 8. 1.2 

Variable 

-l' 

Distance from wing trailing edge to centroid 
of aft ix segments, m (in) 

Figure 4. 3. 1.2 

Variable 

C 

r 

e 

Root chord of exposed wing panels, m (in) 

Table 2.1.1 

1.495 (58.87) 

c 

n 

Chord of wing at centerline of nacelle, m (in) 

Figure 4. 8. 1.2 

1.331 (52.4) 

\ 

Distance from wing trailing edge to centroid 
of last aft Ax segment, m (in) 

Figure 4. 8. 1.2 

3.581 (141.0) 

9e /3a 

V 

Rate of change of downwash behind wing 

Table 3.9. 3. 2 

= 0.4 

36/3a 

•For Ax segments forward of wing leading edge: 
the variation of upwash at segment with angle 
of attack 

Figure 4. 8. 1.1 

Variable 


•For Ax segments aft of wing trailing edge: 
the variation of downwash at segment with angle 
of attack 

Eq. (4. 8.1.7) 

Variable 

s 

w 

Reference wing area, m^ (ft^) 

Table 2.1.1 

14.4 (155.0) 

c 

w 

Reference mean aerodynamic chord, m (ft) 

Table 2.1.1 

1.225 (4.018) 

a 

Lift curve slope of wing, deg~^ 

Table 4. 2. 4.1 

0.0883 


111 



Table 4. 8. 1.2; Concluded 


Figure Figure 
4. 8. 1.2 4. 8. 1.2 



Summary: (C^ 


13.554 

22.759 

28.952 

31.221 

44.840 

30.287" 

18.001 

9.303 

2.872 


36.5 S c. 

w w o 


Xj_, in 

Xj^', in ' 

Figure 
4. 8. 1.2 

Figure 
4. 8. 1.2 

88.8 

- 

70.0 

- 

51.3 

- 

30.8 

- 

20,5** 

- 

- 

15.9 

- 

47.3 

- 

78.5 

- 

109.8 

- 

141.0 

, . 38 

3^r = 

219.586 

36.5 S E 
w 


JlL 


C 

r 

e 


1,503 

- 

1.139 

- 

0.871 

- 

' 0.523 

- 

0.348 

- 

- 

0,113 


0.335 

- 

0.557 

- 

0.779 

- 

1.0 



38. 

3a 

Eq. 

(4. 8. 1.7) 

2. 38 
ft3 

- 

4.144 

- 

17.092 

- 

30.634 

- 

43.423 

- 

103.061 

0.068 

3.049 

0.201 

6.088 

0,334 

6.012 

0.467 

4.344 

0.600 

1.723 

219.586 


= 0.00966 


Segment 4x, in w^, in in x^^’, in 

Figure Figure , Figure 

4. 8. 1.2 4. 8. 1.2 “ 4. 8.1. 2 




36.* ^ 

3a 3a 

Figure Eq. 

4. 3.1.1 (4. 8. 1.7) 



2« ou 
w ■‘Axr — 
n 3a 


7.828 

12.001 

12.870 

14.280 

36.889 


83.868 


• Summary: (C ) . = 2 4-— 5 w 2 a x = 2 = 0.00737 

“a 36.5 Sc " 36.5 S 5 

w w o w w 


* Including (Cj ) correction: 


** For segment 5: = Xj^ 


M = fa.) 

3a \ 3a /, 


(Cj^ 0,0785 
a 


g 

0.0785 
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Table 4. 8. 4.1: Pitching moments of the ATLIT airplane with horizontal 
tail removed 


a, deg. 

wfn 

Table 4 . 4 . 3 . 3 

wfn 

^^m ^ ^0.25c 
wfn 

Eq. (4.8.3. 2) 

-4 

-0.07803 

0.00609 

-0.12560 

-2 

0.11405 

0.01301 

-0.08061 

0 

0.30661 

0.09401 

-0,03651 

’ 2 

0.49966 

0.24966 

0.00669 

4 

0.69319 

0.48051 

0.04899 

15.9 

1.57308 

- 

*0.16742 


* With Equation (4.8.3.12) 










Figure 4, 8, 4.1: 


Comparison of predicted pitching moments with wind 


tunnel data (horizontal tail off, propellers removed 
2,3 million) 


See note on Page 109 







4.9 Downwash and Dynamic Pressure at the Horizontal Tail 


The method presented- 'for' predicting the subsonic downwash and dy- 
namic pressure in the region of the tail plane was developed for the 
linear lift region for swept- and unswept-wing airplanes. This 
method, however, provides a reasonable approximation for the doxrawash 
and dynamic pressure in the nonlinear lift range below stall. 


4.9.1 Downwash at the Horizontal Tail liEPRODUCffilLITY OP THS 

ORIGINAL PAGE IS POOR 

The method discussed is also in References 3 and 4. The method 

is applicable to configurations in which the span of the \-7ing is at 

least 1.5 times as large as that of the horizontal tail (b /b, > 1.5). 

w h — 

The basic approach is as follows; 

1) Detemnine the do^mwash in the plane of symmetry at the height 
of the vortex cores and at the longitudinal station of the quarter 
chord point of the horizontal tail mean aerodynamic chord. 

2) Correct this value for the horizontal tail height above or below 
the trailing vortices. 

3) Evaluate the effect of horizontal tall span by relating the 
average do^TOwash at the tail to the downwash determined in Step 2. 

The downwash gradient, 9e/9a, at the trailing edge of the wing 
is unity. The value at a distance infinitely for downstream is given 
by: 


9e 


CO 


3a 


360 

tt^A 

w 


a 


(4. 9. 1.1) 


where 

(C^ )^ is the wing lift curve slope per degree. If these two 
a 

values are known, the downwash gradient for any intermediate longitudinal 
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position can be found by means of the lifting line theory. 

For straight wings that have tip stall or thin swept wings that 
shed the leading edge vortices inboard of the wing tips, the effective 
wing aspect ratio is considerably less than the geometric aspect ratio. 
An effective aspect ratio is presented for these wings. 

Because of the spanwise variation of downwash, the effective 
downwash acting on the horizontal tail is different from that at the 
plane of symmetry. A, correction for tail-span effect is presented 
and it is based on the assumption that the vortices are essentially 
rolled up at the longitudinal tail station. 

The subsonic downwash is obtained from the following precedure 
(see also Figure 4. 9. 1.1): 

1) Determine the effective wing aspect ratio, A , and the effective 

^eff 

wing span, b , from Figure 4. 9. 1.2 as a function of the wing angle 
eff 

of attack parameter: 


a , - a 
w o 


w 




) - Ct 

w o 


w 


max 

2) Determine the downwash gradient, t^^yc^^^^low speed plane 

of symmetry at the height of the vortex core from Figure 4. 9. 1.3. 

This figure is entered with 2i„/b and A , where is the distance 

z w w z 

eff 

measured parallel to the wing root chord between the aft end of the 
wing root and the quarter chord point of the mean aerodynamic chord' 
of the horizontal tail. 

3) Determine the vertical position, a, of the quarter chord point of 
the mean aerodynamic chord of the horizontal tail relative to the 
vortex core. This distance depends upon the type of wing flow sepa- 
ration as determined from Figure 4. 9. 1.4. For trailing edge separation: 
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a IT 

^ " ^h' ■ ^eff Vl^ 


0.41 C. 


w 


w 


eff 


IT A 


tan r 


(4. 9. 1.2) 


w 


For leading edge separation: 


^ ^3^ I 180 


0.41 C. 


w 


TT A 

W 


eff 


where 


(4, 9, 1.3) 


is the distance measured parallel to the wing root, between 
the effective wing tip quarter point and the quarter chord point of 
the horizontal tail mean aerodynamic chord, 
r is the dihedral angle of the wing 

z^’ is the height of the horizontal tail mean aerodynamic chord 
quarter chord point above or below the plane of the wing root chord, 
measured in the plane of symmetry and normal to the extended wint 
root chord, positive for the tail above the plane of the mng root 
chord. 


is the distance measured parallel to the plane of symmetry, 
between the leading edge of the wing mean aerodynamic chord to the 
trailing edge of the wing root chord. 

5,2 is defined in Step" 2. " ' ' 

4) Calculate the span of the vortices at the longitudinal location 
of the quarter chord point of the horizontal tail mean aerod 3 mamic 
chord by: 


b = b 
vor w 


eff 


- (b, 


- b ) 

"eff 


' 2 % _ 
eff 

b 5, 
w ru 


where 


(4. 9.1.4) 


b^^, the span of the completely rolled up vortices, is obtained 

from: 
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(4. 9.1.5) 


= [0.78 + 0.10 a - 0.4) + 0.003 (A ,,) ] b 
ru w c/4 w w 


eff 


where 


is the quarter' chord sweep angle of the wing in degrees 


and: 


0.56 A 


^ru n 


w 


(4. 9. 1.6) 


w 


5) Determine the average downwash gradient acting on the tail by: 


9e. 

( — 

3 a low speed 


r 3s^/3a 1 


3e 


3 3a j low speed ^ 3a "^lowspeed 


(4.9.1. 7) 


where the [ ] quantity is obtained from Figure 4. 9.1.5 by using the 
parameters calculated in the above steps. 

6) For high subsonic Mach number the do\fnwash gradient is given by: 

3e. ai ^w 

" ^"^low speed (C^ )^ (4.9.1. 8) 

a low speed 


7) Determine the average downwash acting on the horizontal tail by 
integrating the average do^mwash gradient from Step 5 or 6 as follox<rs; 

(4.9.1.9) 

o 


4.9.2 Dynamic-Pressure Ratio at the Horizontal Tail 

In the linear angle of attack range, the ratio of the dynamic pres- 
sure in the plane of symmetry at some distance x aft of the wing root 
chord trailing edge to the free stream dynamic pressure, q^^/q^, is ob- 
tained from the procedure outlined in the following steps (see also Fig- 
ure 4. 9. 2.1). 
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1) Calculate the half width of the wing wake by: 


= 0-68 c VC^ + 0.15) 

wake w V ]) 

f w 


(4. 9.2.1) 


where C^_ is the wing zero-lift drag coefficient, 
f 

2) Calculate the downwash in the plane of symmetry at the vortex sheet 

(assumed to be the same location as the wake centerline) by: 

1.62 C, 




180 

ir 


w 


ttA 


(degrees) 


( 4 . 9 . 2 . 2 ) 


w 


3) Determine the vertical distance, z. ", from the vortex sheet to the 

ti 

point of interest (usually the quarter chord point of the mean aerodyn- 
amic chord of the horizontal tail) by; 


z, " = X tan (y- + s, - a ) 
h ' h w 


(4. 9. 2. 3) 


where y is defined in Figure 4.9. 2.1. 

4) Determine the dynamic pressure loss in the wake center by: 


^co 0 


2.42 C. 


0.5 




X 


(4. 9.2.4) 


+ 0.30 


w 


5) Determine the dynamic pressure loss ratio for points not on the 
wake center' line by: 


_ ; . 4 q . 2 ^ tt Tr. 

” h “> n wake 


) 


(4. 9. 2. 5) 


where ( — -r- ) is expressed in radians. When- 2, "/A2 , is greater 

2 Az - ' ^ h wake “ 

wake 

than one, the dynamic-pressure ratio loss is zero. 
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6) Determine the dynamic-pressure ratio in the plane of symmetry at an 
arbitrary distance x (usually the quarter chord point of the mean aero- 
dynamic chord of the horizontal tail) aft of the wing root chord trail- 
ing edge by: 

^ = 1 - (M ) (4. 9. 2. 6) 

4.9.3 Downwash and Dynami.c-Pressure Ratio at the Horizontal Tail of 
the ATLIT Airplane 

Tables 4. 9. 3.1 and 4. 9. 3. 2 show the calculations to obtain the down- 
wash angle at the horizontal tail as function of the angle of attack 
while the result is shown in Figure 4. 9. 3.1. 

In Figure 4. 9. 3.1 also the downwash angle at the tail obtained 
from the experimental wind tunnel data is shown. Appendix E demonstrates 
in what manner these results were obtained from the wind tunnel data of 
Reference 2. The slope of the predicted curve shows fair agreement- 
with the slope fo the experimental downwash curve. However, in the lin- 
ear lift range the predicted value of is approximately 1.5 degrees 
too small. 

A summary of the calculations which lead to the dynamic pressure 
ratio at the horizontal tail is listed in Table 4. 9. 3. 3 and 4. 9. 3. 4. 

The results indicate that the horizontal tail, except at high angles of 
attack, stays out of the wake. Therefore, the dynamic pressure ratio 
at the horizontal tail will be assumed to be constant and equal to one. 

The result of the calculations is shoim in Figure 4.9. 3. 2. 
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In, Section 4.13- the effect of deflection of the stabilizer on 
the lift and pitching moment of the ATLIT is predicted. From the re- 
sults of that section it seems that the assumption: 

is too optimistic. In Reference 8 the following values are listed: 
may be assumed equal to 0.85 for a fuselage mounted stabilizer 
and 0.95. for a fin mounted stabilizer, except in the case of a T-tail 
(with- q^/ q^ - 1) . For the ATLIT airplane this would mean that the 
dynamic pressure ratio, > is equal to 0.85 instead of equal to 

1 (one) . 
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Table 4. 9. 3.1: Pertinent parameters for computing average downwash 
at the horizontal , tail of the ATLIT airplane. “• 


j Symbol 

DescrxpCion 

Reference 

Magnitude 

(a„)„ 

o w 

Wing zero-lift angle relative to chord, deg* 

Table 4. 2. 4.1 

-2.9 

^ C* 'w 
ciax 

Wing angle of attack at , relative to wing 

chord, deg. max 

Table 4.2. 4.1 

16.4 


Wing sweep angle along c/4 line, deg» 

Table 2.1.1 

1.835 

Ay 

Wing leading edge sharpness parameter, perc* of 
chord 

Table 4.1.2 

4.9 

Xw 

Wing taper ratio 

Table 2.1.1 

0.5 

A 

V 

Wing aspect ratio 

Table 2.1.1 

10.32 

b 

w 

Wing span, a(ft) 

Table 2.1.1 

12.19(40.0) 

r 

Wing dihedral, deg. 

Table 2.1.1 

7.0 

h.'* 

Perpendicular distance froa'wing-root chord plane 
to c/4 of horizontal tail,- m(ft) 

Rigure 4.9.1.1 

0.832(2.73) 

"2 

Tall length in wing-root chord plane from wing- 
root, trailing edge to c/4 of horizontal tail, m(fc) 

Figure 4, 9. 1.1 

3. 716 a2. 19) 

s 

Distance from leading edge of wing, mean aerodynasdc 
to trailing edge of wing root chord, m(ft) 

Figure 4. 9. 1.1 

1.414(4.64) 

A /A 

w w 

eff 

Ratio of effective to geometric wing aspect ratio 

Figure 4. 9. 1.2 

1.0 

b , /b 
erf 

Ratio of effective to geometric wing span 

Figure 4. 9. 1.2 

1.0 


Tail length in root chord plane from vortex tip 
of c/4 of b to c/4 of horizontal tail, raCft) 

eff 

Figure 4. 9. 1.1 

4.700a5.42) 


Lift curve slope of wing, (deg) ^ 

Table 4. 2. 4.1 

0.0888 

3c^/3o 

Downwash gradient at infinity 

Eq.(4.9.1.1) 

0.314 


— 

— 

0.6095 

(3e^^/3a)l^„ speed 

Downwash gradient in plane of symmetry at height 
of vortex core 

Figure 4.9.1.3 

0.395 

b 

ru 

Span of complete rolled up wing tip vortxces, m(ft) 

Eq. (4.9.1.5) 

7.714(31.82) 

£ 

ru 

Distance required for complete rollup of wing 
tip vortices in chord plane, semispans 

Eq,C4-9.1.6) 

5,7792/Cj^ 

w 

b 

vor 

Span of vortices at longitudinal station^ of c/4 
of horizontal tail, m(ft) 

Eq.C4.9.1.4) 

Variable 

a 

Vertical distance from horizontal tall root chord 
to vortex core if trailing edge separation, o(ft) 

Eq.(4.9.1.2) 

Variable 

isjia 

C ^ 1 

low speed 

Ratio of average downwash acting on horizontal tail 
to downwash at vortex core height 

Figure 4. 9, 1.5 

Variable . 

-( speed 

Downwash gradient 'at horizontal taxi at low speed 

Eq.(4.9.1.7) 

Variable 

( 3Sj^/3a)jj 

Downwash gradient at Mach number 

Eq.(4;9.1.8) 

Variable 

"h 

Downwash at horizontal tail, deg. 

Eq.C4.9.1.9) 

Variable 
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Table 4, 9. 3. 2; Calculation of downwash at the horizontal tail of the 
ATLIT 


a, deg 
© 

. 

CE^, relative to chord 

abs 

deg 

% 

Figure 4.2. 4.1 
© . 

a, ft . 
Eq.(4.9.1.2) 

© ■ 

1 

ru 

Eq.(4,9.1.6) 

© 

deg 

0 

rad 

0 

-4 

-3,5 

-0.0611 

-0.6 

-0.0533 

1.206 

-108.43 

-2 

-1.5 

-0.0262 

1.4 

0.1243 

0.702 

46.49 

0 

0.5 

0.0037 

3.4 

0,3019 

0.199 

19.14' 

2 

2.5 

0.0436 

5.4 

0.4795 

-0.305 

12.05 

4 

4.5 

0.0785 

7.4 

0.6571 

-0.808 

8.80 

15.9 

16.4 

0.2862 

19.3 

1.494 

-3.848 

3.87 


b 

vor 

Eq.(4.9,1.4) 

© 

2a/b 

vor 

© 

h vor 
© 

3e^/3a 

speed 
Eq. (4.9.1. 7) 

0 

to 

01 

11 

0 

3e /3a 

VC 

Figure 4. 9.1.5 
0 

38.95 

0.036 

0.348 

.■1.019 

0.403 

0.564 

38.36 

0.010 

0.353 

1.025 

0.405 

1.377 

37.93 

-0.016 

0.358 

1.027 

0.406 

2.192 

37.58 

-0.043 

0.361 

1.023 

0.404 

2.990 

36.35 

-0.212 

0.373 

0.950 

0.375 

7.238 
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Table 4. 9, 3. 3: Dynamic— pressure ratio at the horizontal tail of the 


ATLIT 


Syrabol 

Description 

Reference 


a 

w 

Wing angle o£ attack relative to root chord = 

a + i , deg. 

W o 

Table 4. 4. 3.1 

a + 0,5 

Y 

Angle between wing chord plane and line connecting 
trailing edge of wing root chord and c/4 of hori- 
zontal tail, deg. 

figure 4. 9.2.1 

12.9 

c 

w 

Wing mean aerodynamic chord, m(ft) 

Table 2.1.1 

1.225 (4.018) 

A 

w 

Wing aspect ratio 

Table 2.1.1 

10.32 

X 

Distance from trailing edge of wing root chord 
to c/4 of horizontal tail measured along center- 
line wake, m(ft) 

Figure 4.9.2.X 

3.65(11.98) 

% 

Wing zero-lift drag coefficient of total wing per 
procedure of Section 4.12 

— 

0.0111 


Downwash in plana of symmetry at vortex sheet, deg. 

Eq.(4.9.2.2) 

2,863 

w 


Vertical distance from vortex sheet to c/4 of hori- 
zontal tail, m(ft) 

Eq.(4.9.2.3) 

Variable 

&z . 
wake 

Half width of wake, m(ft) 

Eq.(4.9.2.1) 

0.155(0.5094) 


Dynamic pressure loss in the wake centerline 

Eq.(4.9.2.4) 

0.0777 


Dynamic pressure loss at the horizontal tail 

Eq.(4.9.2.5) 

Variable 


Dynamic pressure ratio at the horizontal tail 

Eq.(4.9.2.6) 

Variable 


a, deg 

a ,deg 

W 

% 

Figure 4. 2. 4.1 

"h’ 

Eq.(4.9.2.2) 

2 *• ft* 

Eq.(4.g.2.3 

"h”^^"wake 

Eq.(4.9.2.5) 

Eq.(4.9.2.6) 

-4 

-3.5 

-0.0533 

-0.1526 

3.49 

6.85 . 

0. 

1. 

-2 

-1.5 

0.1243 

0.3559 

3.04 

5.97 

0. 

1. 

0 

0.5 

0.3019 

0.8643 

2.83 

5.56 

0. 

1. 

2 

2.5 

0.4795 

1.3728 

2.49 

4.89 

0. 

1. 

4 

4.5 

0.6571 

1.8813 

2.17 

4.26 

0. 

1. 

14 

14.5 

1.465 

4.194 

0.54 

1.07 

0. 

1. 

15 

15.5 

1.490 

4.266 

0.35 

0.68 

0.018 

0.982 

15.9 

16.4 

1.494 

4.277 

0.17 

0.33 

0.059 

. 0.941 
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Figure 4. 9. 1.2: Effective wing aspect ratio and span for low speeds 
(Reference 3) 
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Predominant trailing-edge 
separation • 





Predominant leading-edge 
separation 


0 10 20 30 40 50 

(Ac/ 4)^^. deg 

Figure 4. 9,1,4: Type of flow sep.aration as a function of airfoil 
and wing sweep (Reference 3) 


Figure 4, 9,1, 5: Average downwash acting on aft lifting surface 
for low speeds (Reference 3) 






131 





Figure 4. 9. 3.1: Comparison of predicted downwash at horizontal tail 
with full-scale wind tunnel data (propellers removed, 
Nj, = 2,3 million) . ’ 
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Figure 4. 9. 3. 2; Prediction of dynamic- pres sure ratio at horizontal tail 
(propellers removed) 


133 



4.10 Lift of the Complete Airplane 


In this section the lift of the complete airplane will be dis- 
cussed. Lift due to flap-and-elevator-def lection will not be included. 
The method discussed is similar to the method used in Reference 3. 


4.10.1 Lift in the Linear Lift Range 


The net lift of the complete airplane in the linear lift range 
may be written as follows: 


C. = C, 


+ C 


■“win 


^i(hf) 


^^^L^h(fv) 


(4.10.1.1) 


where 


C is the tail-off lift coefficient considered in Section 4.4. 
wfn 

C is the lift contribution of the horizontal tail including 

\(hf) 

tail-fuselage interactions, wing downwash and dynamic-pressure effects. 

is the effect of fuselage vortices on the lift coeffic- 
ient of the horizontal tail. 

The lift contribution of the horizontal tail in the presence of 
the fuselage due to angle of attack at the tail can be estimated as 
follows : 


^(hf) 


a e 




+ K 


f(h) 


a. 



S q 
w ~ 


(4.10.1.2)' 


where 


«L>h 
a e 


is the lift curve slope of the exposed panels of the tail. 
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is the ratio of the lift on the horizontal tail in the 
presence of the body to tail alone, obtained from Figure 4.4.1. 1. 

is the ratio of the tail-lift carryover onto the body to 

f (h) 

tail alone, obtained from Figure 4. 4. 1.1 

a, is the angle of attack of the horizontal tail: 
h 

= a - (4.10.1.3) 

where 

£, is the average downwash acting on the horizontal tail, obtained 
h 

from Section 4.9. 

i^ is the incidence angle of the horizontal tail. 

S, is the exposed area of the horizontal tail, 
h 

e 

(q^/q^) represents the dynamic-pressure ratio at the horizontal 
tail, obtained from Section 4.9. 

According to Reference 3, the effect of body vortices on the .lift 

of the horizontal tail can be considered negligible when the horizontal 

tail span is greater than three times the body diameter at the tail 

(b._ > 3(d^).^). This ratio is exceeded by almost all general aviation 
ti r £1 

aircraft. 

4.10.2 Maximum Lift of the Complete Airplane 


The T Tiaximnm lift coefficient of the horizontal tail in the pres- 


ence' of the fuselage, (C^ ^h(hf)’ corresponding angle of 

max 

attack, (a_ ), , are determined by the methods used in Section 

h(,nr; 

L 


max 

4.4 to obtain the maximiim lift characteristics of the wing in presence 
of the fuselage. The method uses an empirical taper ratio correction 
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factor, which can be obtained from Figure 4. 4. 2.1 and may be used 
up to M = 0.6: 


^h(hf) 

max 


^h(hf) 
max 

^h 

max 


(Cl (.-f-) (4.10.2.1) 


max 


w 


r (ct. 




^h(hf) 


max 


^h(hf) 


max 


(“r 


(“r 


(4.10.2.2) 


max 


max 


where 

(C^ ), is the maximum lift coefficient of the horizontal tall 

L h 
max 

based on the horizontal tail area, obtained from Section 4.2. 

(a„ ), is the horizontal tall stall angle relative to the 

max 

chord, obtained from Section 4.2. 

>h(h£)^<‘^L >h >h(hf)'('‘c, >h 

max max L L 

max max 

correction factors obtained from Figure 4. 4. 2. 3 and 4. 4. 2. 2, respectively, 
as a function of (c^ + ^*^f^h^^h* 


4.10.3 Lift Characteristics of the ATLIT Airplane 


In Table 4.10.3.1 the svimmary calculations for the tail-lift in 
the linear lift range are presented. A s umm ary of the calculations of 
the maximum tail-lift is listed in Table 4.10.3.2, while in Table 
4.10.3.3 the total lift predictions of the ATLIT are shorn. The result- 
ing lift curve is shown in Figure 4.10.3.1 and compared with the ex- 
perimental results of Reference 2. In this figure also, lift curve 
obtained with Reference 5 is shown. 

In the linear lift range the slope of the hand-calculated lift ' 
curve shows good agreement with the slope of the experimental lift^curve. 
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The angle of attack at zero-lift prediction, however, is more than one 
degree off. 

The lift curve obtained with Reference 5 shows fair agreement 
with the full-scale wind tunnel data, while no results were obtained, 
with the computer program of Reference 6 . 

In Figure 4.10.3.2 the effect of the engine cooling system on the 
lift coefficient is shown. In this report no attempt will be made to 
predict this effect because of uncertainties in the determination. 
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Table 4.10.3.1; Lift of horizontal tail in linear lift range in presence 
■ of fuselage 


Symbol 

Description 

Reference 

Magnitude 

“h* 

Limit of linearity of horizontal tail lift curve, 
relative to tali chord, deg. 

Table 4. 2. 4.1 

12; 2 


Fuselage width at horizontal tall, m(ft) 

Figure 2.1.2 

0.405 (1.33) 


Horizontal call span, m(ft) 

Table 2.1.1 

4.133 (13.56) 


— 

— 

0.098 

'Si(f) 

Ratio of lift on tall In presence of fuselage to 
tail alone 

Figure 4. 4. 1.1 

1.082 

' ^f(h) 

Ratio of tail-lift carryover on fuselage to tail 
alone 

Figure 4. 4. 1.1 

0.130 

>h 

Lift-curve slope of exposed horizontal tail 

Table 4. 2. 4.1 

0.0649 

a e 

panels, deg"*^ 



e 

2 2 

Area of exposed horizontal tail panels, m (ft ) 

Table 2.1.1 

3.24 (34.9) 

S 

w 

2 2 

Reference wing area, ni (ft ) 

Table 2.1.1 

14.40 (155.0) 


Horizontal tall incidence angle, deg. 

— 

0 


Dynamlc-pressur.e ratio at the horizontal tail 

Figure 4.10.3.2 

1.0 

Summary; C = 0.0177 (a - e ) based on S and up Co a._* = 12.2 

hi(hf) ^ w h 

degrees 
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REPRODUCIBILITY' OF THE 
ORIGINAL PAGE IS POOR 

Table 4.10.3.2: Maximum lift- characteristics of horizontal tail in 

presence of the fuselage 


Symbol 

Description 

Reference 

Magnitude 


Horizontal tail leading-edge sweep angle, deg. 

Table 2.1.1 

0. 


Horizontal tail aspect ratio 

Table 2.1.1 

4.75 


Horizontal tail taper ratio 

Table 2.1.1 ' 

1 


Ratio of body width to tall span at tail 

Table 4.10.1 

0.098 

max 

Horizontal tail maximum lift coefficient, based 

on S, 
h 

Table 4. 2. 4.1 

1.215 


Horizontal tail stall angle, deg. 

Table 4. 2. 4.1 

18.75 

max 





Taper ratio correction factor 

Figure 4. 4. 2.1 

0.0852 

(c2+l)Aj^Can(A^^)j^ 

— 

— 

0 


^h(hf) 

max 






L max J 



Figure 4.4. 2. 3 

1.0 


^“c ^h(hf)' 

• max 






t‘\ . 

max 



Figure 4. 4. 2. 2 

1.030 


2 2 

Horizontal tail area, m (ft ) 

Table 2.1.1 

3.60 (38.7) 

s 

V 

2 2 

Reference wing area, m (ft j 

Table 2.1.1 

14.40 (^5,0) 

max 

Maximum lift coefficient of the horizontal tail 
in the presence of the fuselage based on S 

w 

Eq. (4.10.2,1) 

0.303 

^h(hf) 

max 

Angle of attack for zero horizontal tail lift, 
with respect to chord of horizontal tail, deg. 

Eq. (4.10.2.2) 

19.3 
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Table 4.10,3.3: Lift of the ATLIT airplane (no flap and stabilizer 

deflection) 


ot, deg. 

mmmmm 

du. deg. 

C 

C. 

c 



h 

^(hf) 

^fn 

L 



Eq. (4.10,1.3) 

Table -4.10.3.1 

Table 4. 4. 3, 3 

Eq. (4.10.1.1) 

-4 

-0.26 

-3.74 

-0.0662 

-0.0780 

-0.1442 

-2 

0.56 

-2,56 

-0.0453 

0,1141 

0.0688 

0 

1.38 

-1.38 

-0.0244 

0.3066 

0.2822 

2 

2.19 

-0.19 

-0.0034 

0.4997 

0.4963 

4 

2.99 

1.01 

0.0179 

0.6932 

0.7111 
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7.24 
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0.1533 
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Figure 4,10.3.1: Comparison of predicted airplane lift curve with 

wind tunnel data (propellers removed, 2.3 million) 
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Figure 4.10.3.2: Effect of engine cooling on airplane lift coefficient 

3.5 million) 


142 







4.11 Pitching Moment of the Complete Airplane 


The pitching moment of the complete airplane may be estimated 
as follows : 


C = C + C 
™ ”wfn %(hf) 


(4.11.1) 


where 


C is the tail-off pitching moment coefficient considered in 
wfn 


Section 4.8. 


C is the contribution of the horizontal tail, including 

“h(hf) 

tail-fuselage interference effects , to the pitching moment coefficient 
of the airplane. The pitching moment due to the horizontal tail is 
determined from: 

C = (-^ - ~ Vt (4.11.2) 

^(hf) V c c / n(hf) 

' w w ' 

where 

X is the distance, parallel to the X-body axis, from the leading 
eg 

edge of the wing mean aerodynamic chord to the airplane center of 
gravity (positive aft) 

Xj^ is the distance, parallel to the X-body axis, from the leading 

edge of the wing mean aerodynamic chord to the quarter chord of the 

horizontal tail mean aerodynamic chord- (positive aft) 

C is the lift of the horizontal tail in the presence of the 

Hi(hf) 

fuselage, obtained from Section 4.10 

Summary calculations for the pitching moment of the ATLIT airplane 
relative to .the quarter chord of the wing mean aerodynamic chord are 
presented in Tables 4.11.1 and 4.11.2. The resulting pitching moment 
curves are shown in Figure 4.11.2 and 4.11.3 as function of the angle 
of attack and the total lift coefficient, respectively. 
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The results from Table 4.11.2 show good agreement with the ex- 
perimental results of Reference 2. 

In Figures 4.11.4 and 4.11.5 the pitching moment contribution 
due to the engine^ coolihg.' system is shown. These curves are obtained 
from Reference 2, and they demonstrate that the contribution of the 
cooling system to the pitching moment is not negligible. An analytical 
treatment of lift (pitching moment) due to the cooling system is 
beyond the scope of this report because of uncertainties in the 
determination . 
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Table 4.11,1; Pitching moment contribution of the horizontal tail 
(stabilizer not deflected) 


Symbol 

Description 

Reference 

Magnitude 

X /c 
eg w 

Airplane center of gravity location from leading 
edge of wing mean aerodynamic chord 

- 

0.25 

X. 

a 

Distance_f Torn 'Wing mean aerodynamic chord leading 
edge to c/4 of horizontal tail, ra (ft) 

Figure 4.11.1 

5.127 (16.82) 

c 

w 

Wing mean. aerodynarndLc (geometric) chord, m (ft) 

Table 2.1.1 

1.225 (4.018) 

V=w 

- 

- 

4.186 


Lift coefficient of horizontal tail on basis of S 

w 

Table 4.10.3.1 

Variable 

”h(hf) 

Pitching moment coefficient of total airplane 

Eq. (4.11.2) 

-3.936 (1 

^(hf) 


Table 4,11,2; Pitching moment of complete airplane- 


o, deg 

Table 4.10.3.3 

^^(hf) 

Table 4.10.3.1 

"n(hf) 
Eq. (4.11.2) 

C 

®wfn 

Table 4. 8. 4.1 

C 

m 

Eq. (4.11.1) 

-4 

-0.1442 

-0.0662 

0.2606 

-0.1256 

0.1350 

-2 

0.0688 

-0.0453 

0.1783 

-0.0806 

0.0977 


0.2822 

-0.0244 

0.0960 

-0.0365 

0.0595 


0.4963 

-0.0034 

0.0134 

0.0067- 

0.0201 


0.7111 

0.0179 

-0.0705 

0.0490 

-0.0215 


1.7264 

0;1533 

-0.6034 

0.1674 

-0.4360 
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Figure 4.11.3: Comparison of predicted airplane pitching moment with 

full-scale wind tunnel data (propellers removed, stabilizer 
not deflected, 2.3 million) 
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Figure 4.11,5: Effect of engine cooling on airplane pitching moments 
(N^^= 3.5 million) 
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4.12 Drag of the Complete Airplane 


The drag of an airplane can be split up into the following compo- 
nents : 

(1) Zero-lift drag of the wing, horizontal tail and vertical tail 

(2) Zero-lift drag of the fuselage and nacelles 

(3) Zero-lift interference drag of the wing-fuselage, tail- fuselage, 
vertical tail-horizontal tail and wing-nacelle 

(4) Lift induced drag of the wing and horizontal tail 

(5) Lift induced drag of the fuselage and nacelles 

(6) Lift induced wing-fuselage and wing-nacelle interference drag 

(7) Cooling drag due to nacelle inlets and cooling flaps. 

4.12.1 Zero-Lift Drag of Wing, Horizontal Tail and Vertical Tail 

The zero-lift drag of a lifting surface is composed of skin- 

friction drag and pressure drag. The pressure drag is caused by the 

boundary layer which prevents complete pressure recovery at the 

trailing edge of the lifting surface. The skin friction drag is 

a function of the roughness of the surface of the lifting surface 

and the type of boundary layer flow. In the case of general aviation 

aircraft the boundary layer may be considered to be turbulent. 

The following expression is used in Reference 3 to predict the 

zero-lift drag of a lifting surface based on the reference wing area, 

S : 
w 

p if -j J2.S 

hs - 2 Cj [l + 2 (i) + 120 (f> I -5-^ (4;12.1.1) 

o w 
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where 


Cj is the skin-friction coefficient of a flat plate, obtained from ■ 

Figure 4.12.1.1 as a function of Reynolds number, N_ , and the parameter 

R0 

5,/k 

i is the reference length of the lifting surface, in this case 

the mean aerodynamic chord of the exposed area of the lifting surface 

k is the surface roughness height, estimated from Table 4.12.1.1 

on the basis of surface finish 

t/c is the thickness ratio of the lifting surface 

represents the exposed area of the lifting surface 
e 

The summaa^^ calculations of the zero-lift drag of wing, horizontal 

tail and vertical tail are listed in Table 4.12.1.2, and the final 

results based on the wing area are : . 

wing (Cj^ = 0.00970 

o 

horizontal tail = 0.00221 (4.12.1.2) 

o 

vertical tail (Cj^ ) = 0.00111 

o 

4.12.2 Zero-Lift Drag of Fuselage and Nacelles 

The zero-lift drag of an isolated body is composed of skin-friction 
drag, pressure drag and base drag. Equation (4.12.2.1), which estimates 
the zero— lift drag of an isolated body, is only valid for axixymmetric 
bodies of revolution. Therefore the fuselage and nacelles have to be 
treated as equivalent bodies of revolution having an axially distrib- 
uted circumferential area similar to that of the actual body. For 
subsonic conditions (M ^ 0.6) and on the basis of the reference wing 
area, S^, the zero-lift drag coefficient of an isolated body is: 
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(4.12.2.1) 


O' f b 


o 


1 + 


60 




400 


S ^ V 
w b 


where 


is the skin friction coefficient of a flat plate, obtained from 

Figure 4.12.1.1 as a function of Reynolds number, based on the 

actual body length, and the parameter A„/k 

D i5 

is the actual body length 
li 

k is the surface roughness height, obtained from Table 4.12.1.1 

dg is the diameter of a circle having the same circumference as 

the circumference of the maximtam frontal area ^ 

S ^ is the net wetted surface area of the body and can be esti- 
wet 

mated from Figure 4.12.2.1 as a function of the body fineness ratio, 
£_/d_, and the parameter, d, /d„ 

D iJ D D 

d^ is the diameter of the equivalent circular circumference of 
the base area of the body 

is the maximum frontal area of an axisymmetric body having 

D 

a diameter, d„, or: 

B' , 


S 


B 


- d 2 

4 ^B 


(4.12.2.2) 


is the base drag coefficient which can be written as follows 
b 

(based on S ) : 
w 

= 0.029 (^) l/c (4.12.2 

^ \dg/ r 

In Table 4.12.2.1 the calculations to obtain the zero-lift drag of 
the fuselage and nacelles are summarized. The net wetted area of the 
fuselage and nacelles are. summarized and the resultsc^based on the refer- 
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ence wing area, S , are: 
w 

fuselage (C )- = 0.00616 

o (4.12.2.4) 

nacelle (C^ ) = 0.00614. 

iJ n 
o 

4.12.3 Zero-Lift Interference Drag of Wing-Fuselage, Tail-Fuselage, 
Vertical Tail -Horizontal Tail and Wing-Nacelle 

The correlation factor, will be used to predict the wing- 

fuselage zero-lift drag coefficient. The zero-lift drag of the fuselage 
in the presence of the wing relative to the reference wing area, may be 
written as • follows ; 

^f(w) " ^f \f S (4.12.3.1) 

o f b 

where 

(Cj^ )f is the skin friction and pressure drag coefficient of the 
fuselage obtained from Section 4.12.2 

is the ratio of the wing-fuselage to the fuselage alone zero- 
lift drag with the base drag omitted obtained from Figure 4.12.3.1 

The net zero-lift wing-fuselage coefficient based on the reference 
wing area is : 

^wf " ^^D ^f(w) (4.12.3.2) 

o o o 

where 

(Cp )^ is the zero-lift wing drag obtained from Section 4.12.1 
o 

For the ATLIT airplane the calculations are summarized in Table 

4.12.3.1, and the result is: 

(^)wf = 0.01632/ (4.12.3.3) 

o 

On the basis of the reference wing area, S , the increment of tail 

w 

drag due- to fuselage drag is approximately: 
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2 


1^0.8 (-^)^ - 0.0005 j 


(4.12.3.4) 


where 

is the nximber of junctures of the tail surface with the fuselage 

c^ is the root chord of the exposed panels 
e 

t/c is the thickness ratio of the root chord of the exposed surface 
When the vertical tail intersects with the horizontal tail instead 
of the fuselage’^, the interference drag on the basis of the reference 


wing area , , is ; 


^v(h) 

o 




17 - 0.05 (■ 

c xnt 


-)^ 1 
c-^int J 


(4.12.3.5) 


where 


n 2 is the number of corners at the intersection (normally 1 I 2 = 4) 
(t/c) _.^ ^ is the average thickness ratio of the intersecting sur- 
faces at the intersection 

c. ^ is the chord at the intersection 
int 

For the ATLIT, the horizontal and vertical tail intersect with 
the- fuselage. In Table 4.12.3.2 the calculations are summarized and 
the net zero-lift fuselage-horizontal tail drag based on the reference 


wing area is : 


>h ■ 'S >h + >h(f) 

o o o 


(C^ ), = 0.00224 
D h 
o 


(4.12.3.6) 


while the net zero— lift fuselage-vertical tail drag based on the ref- 
erence wing area is: 


^v ^v ^v(f) 

o o o 


(Cjj )^ = 0.00112 


(4.12.3.7) 
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where 


(Cd represents the zero-lift drag coefficient of the horizontal 
o 

tail , obtained from Section 4.12.1 

(Cjj is the zero-lift drag coefficient of the vertical tail 
o 

from Section 4.12.1 

According to Reference 3 the interference drag of a nacelle 
faired into the firing may be roughly accounted for by the increment 


of zero-lift wing drag due to wing area covered, by the nacelle, or ■ 
for one nacelle: 


(AC^ ^ = (C„. )• 

D n(w) D w 
o o 


w n 


w 


(4.12.3.8) 


where 

(AS^)^ is the wing area overlapped by one nacelle 
The summary calculations of Table 4.12.3.3 show the net zero-lift 
drag of the two nacelles on the basis of the reference wing area, S : 


^n(w) ^ ^n ^n(w) j 

o '■ o o ) 

(^) . . = 0.01384 
n(w) 


(4,12.3.9) 


where. 

(Cjj )^ is the zero-lift drag of one nacelle, obtained- from 
o 

Section 4.12.2 


4.12.4 Lift Induced Drag of Wing and Horizontal Tail 

In Reference 3 the drag of twisted wings due to lift has not 
been discussed. The method for straight- tapered wings of Reference 4 
‘will be used. 

The drag due to lift. of a twisted, swept-back straight wing 
(and horizontal tail) is given by: ’ • 
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(4.12.4.1) 


-■n / A 9 <^0 V + (0 C„ )^W 

D.'^w TrAe L Z I 

X w a a 


(C. )„ = 


where 


is the- wing lift coefficient obtained from Section 4.2 
w 

c.^ is the airfoil section lift-curve slope from Section 4.1 
a 

6 is the wing t^^ist (negative for washout) 

V is the induced drag factor due to linear twist obtained from 
Figures 4.12.4.1 

w is the zero-lift drag factor due to linear twist obtained from 
Figures 4.12.4.2 

e is the span-efficiency factor determined by: 

a 


R (C^ )^ + (1 - R) ttA 
a 


(4.12.4.2) 


where 


(C, ) is the wing-lift curve slope in radians, obtained from 

Ci 

Section 4.2 

R is the leading-edge-suction parameter defined as the ratio of 

leading-edge suction actually obtained to that theoretically possible. 

The parameter is presented in Figure 4.12.4.3 as a function of Mach 

number, wing aspect ratio, wing sweepback and leading-edge-radius 

Reynolds number, (N_ )-. 177,3 is based on the leading edge radius of the 

airfoil at the wing mean aerodynamic (geometric) chord: 

c. 


<\e>LER * • 




w 


(N ) 
Re^MAC 


(4.12.4.3) 


where 


LER is the leading-edge radius of the airfoil as a ratio of the 


chord 

In Table 4.12.4.1 the induced drag calculations of the wing and 
horizontal tail are summarized. The drag of the wing at angle of 
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attack based on the reference wing area, is: 


(Cr, )„ = 


C 2 
w 


"D ^w 29 487 “ 0* 000210- + 0.000264 

^i * w 


(4.12.4.4) 


while the drag of the horizontal tail based on the horizontal tail area. 


S^, IS! 


(C. ). = 




D/h 14.6987 

X 


(4.12.4.5) 


Above expressions are applied in Table 4.12.4.2 to determine the lift 

drag contributions of the two surfaces as function of angle of attack, 

a, and based on the reference wing area, S . 

w 

4.12.5 Lift Induced Drag of Fuselage and Nacelles 

According to Reference 3 the drag of a body at angle of attack 
may be written 

(4.12.5.1) 


“%£f 


where 


C is the lift of the body, obtained from Section 4.3 
Cg ^ is the angle of attack (in radians) of the equivalent 


'ef f 

circular body, which can be determined as follows: 


= a + ct 

Bgff 0 


(4.12. 5.2) 


where 


a is the zero-lift angle of the equivalent circular body relative 
B 

to the reference X-body axis of the airplane. 

The drag of the fuselage of the ATLIT at angle of attack, based 

on the reference wing area, S , is: 

w 


<V>f = \ 


TT (g - 2.5) 
180 


(4.12.5.3) 
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while the drag of the nacelles at angle of attack is: 


^S.^n 180 

1 n 


(4.12.5.4) 


where 

C and C, (total lift of both nacelles) can be obtained from 

J-i Lr 

f n 

Section 4.3. In Table 4.12.5.1 a summary is listed of the calculations 
of the drag of the fuselage and nacelles of the ATLIT airplane. 


4.12.6 Lift Induced Wing-Fuselage and Wing-Nacelle Interference Drag 

In the absence of applicable representative data no attempt is 
made to account for wing-fuselage and wing-nacelle Interference drag 
at angle of attack. 

4.12.7 Cooling Drag 


The discussion of nacelle drag in the previous section did not take 
into account the effect of drag due to the cooling system. An analyt- 
ical treatment of cooling drag is beyond the scope of this paper 
because of the complexity and uncertainties in its determination. 
However, in Reference 2 drag curves are listed of the airplane with 
inlets and cowl flaps open and closed. From these data the increment 
of drag due to the cooling system could be determined and the cooling 
drag is shown in Figure 4.12,7,1. 

4.12.8 Drag Characteristics of the ATLIT Airplane 


The zero— lift drag of the components plus the interference drag 
between components is : 
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(4.12.8.1) 


Ctx = (C^ ) ^ + (C,. ), + (C„ ) + (C_ ) . . 

D D wf D h D V D n(w) 

o o o o o 

= 0.01634 + 0.00224 + 0.00112 + 0.01384 

= 0.03352 

The net drag of the ATLIT is summarized in Table 4.12.8.1 and 
may be written as follows : 

S “ '=D + ^ <S.>£ + (4.12.8.2) 

O 1 X X X 

and the result is plotted .in Figures 4.12.8.1 and 4.12.8.2. In these 
figures the experimental drag data are also presented. The predicted 
drag polar shows good agreement with the experimental results. 
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Table 4.12.1.1: Surface roughness height (Reference 3) 


Type of surface 

k, in. 

Aerodynamically smooth 

0 

Polished metal or wood 

0.02 to 0,08 X 10“^ 

Natural sheet metal 

0.16 X 10“^ 

Smooth matte paint, carefully applied 

0.25 X 10“^ 

Standard camouflage paint, average application 

0.40 X 10"^ 

Camouflage paint, mass production spray 

1.20 X 10"^ 

Dip galvanized metal surface 

6.0 X 10~^ 

Natural surface of cast iron 

10.0 X 10"^ 
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Table 4.12.1.2: Zero-lift drag of wing, horizontal and vertical tail 
of the ATLIT airplane 


Symbol 

Description 

Reference 

Wing 

Horizontal 

Tail 

Vertical 

Tail 


Surface roughness height, m (in) 

Table 4.12.1.1 

0.635x10-2 

(0.25x10-2) 

0.635x10-® 

(0.25x10-2) 

0.635x10-® 

(0.25x10-2) 


Mean aerodynamic chord of exposed 
area of surface, m (in) 

Table 2.1.1 

1.178- 

(46.38) 

0.871 

(34.29) 

1.201 

(47.30) 


- 

- 

1.86x102 

1,37x10® 

1.89x10® 

c 

w 

Mean aerodynamic chord of total 
wing, m(ln) 

Table 2.1.1 

1.225 

(48,22) 

1.225 

(48.22) 

1,225 
• (48.22) 

^^Re^MAC 

Reynolds number based on c^ 

- 

2.3x10® 

2.3x10® 

2.3x10® 


Reynolds number based on 

- 

2.21x10® 

1.64x10® 

2.26x10® 

% 

Skin friction coefficient of 
flat plate 

Figure 4.12.1.1 

3.87x10-2 

4‘. 05x10-2 

3.86x10-2 

(c/c)^ 

Thickness ratio of surface 

Table 4,1.2 

0,17 

0.10 

0.09 

Si 

e 

Exposed area of lifting surface, 
m2 (ft2) 

Table 2.1.1 

12,53 

(134.8) 

3.25 

(34.9) 

1.75 

(18.8) 

S 

w 

Reference wing area, (ft^) 

Table 2.1.1 

14.40 

(155.0) 

14.40 

(155.0) 

14.40' 

(155.0) 

0 

Zero-lift drag, of lifting 
surface on basis of S 

w 

Eq, (4.12.1.1) 

0.00970 

0.00221 

0.00111 








Table 4«12.2,X! ZeirO“lift drag of fuselage and nacelle 


Sycibol 

Description 

Reference 

Fuselage 

Nacelle 


Circumference of the maximum 
frontal area of body m (ft) 

Figure 2.1 

4.69 (15.4) 

3.14 (10.3) 


Diameter of equivalent perimeter 
of oaximufii frontal area = 

Cg^/r, m (ft) 


1.49 (4.9) 

1.01 03. 3) 


Frontal area of equivalent 
perimeter of body, (ft^) 

Sq. (4.12.2.2) 

1.75 (18.86) 

0.79 (8.55) 

i 

i 

Length of body, m (ft) 

Figure 4 . 8 . 1 . 2 

8.35 (27,4) 

2.65 (8.7)* 


Circumference of base area of 
body 

Figure 2.1 

= 0 

= 0 


Diameter of equivalent perimeter 
of base area = Cb^/f, m (ft) 

- 

0 

0 

^^wet^i 

Hetted surface area, m^ (ft^) 

Appendix F 

23.3 (251.0) 

6.08 (65.4)* 

*^i 

Surface roughness height, m (in) 

Table 4.12.1.1 

0.635x10"^ (0.25x10-3) 

0.635x10-3 (0.25x10-3) 

(i/w^ 

Ratio in common units 

- 

1.32x10® 

4.18x10® 

c 

V 

Mean aerodynamic chord of wing^ 
m (ft) 

Table 2.1.1 

1.225 (4.018) 

1.225 (4.018) 

^^Re^MAC 

Reynolds number based on c^ 

- 

2.3x10® 

2.3x10® 


Reynolds number based on 

- 

1,57x107 

4.98x10® _ ’ 


Skin friction coefficient of 
flat plate 

Figure 4.12.1.1 

2.3x10-3 

3,4x10-3 


Reference wing area, m^ (ft^) 

Table 2(1. 1' 

14.40 (155.0) 

14,40 (155.0) 

o 

Zero-lift drag coefficient 
of body based on S^^ 

Eq. (4.12.2.1) 

0.00616 

0,00614 


♦Including 8-inch extension of nacelle 'due to installation of thrust/torque balance. 









Table 4,12.3.1: Net zero— lift drag of wing-fuselage combination 


Symbol 

Description 

Reference 

Magnitude 

(S >w 
o 

Zero-lift drag of isolated exposed wing panels 

Table 4.12,1.2 

0.00970 


Base drag of fuselage 

Table. 4.12.2.1 

0 


Zero-lift drag of isolated fuselage with base 
drag omitted 

Table 4.12.2.1 

0.00616 

M 

Mach number 

- 

0.081 


Length of fuselage, m (ft) 

Table 4.12.2.1 

8.35 (27.4) 

c 

w 

Wing mean aerodynamic chord, m (ft) 

Table 2.1.1 

1.225 (4.01S) 


Reynolds number based on c 

w 

- 

2.3x106 


Reynolds number based on 

- 

1.57x106 

^wf 

Wing-body interference correlation factor 

Figures 4.12.3.1 

1.075 

^wf 

0 

Net zero-lift drag of wing-fuselage combination 

Eq. (4.12.3.2) 

0.01632 


164 




sss,««* 


Table 4.12.3.2: Net zero-lift drag of tail surfaces in presence of 
the fuselage 


Symbol 

Description 

Reference 

Horizontal 

Tail 

Vertical 

Tall 

>h 

0 

Zero-lift drag of isolated exposed horizontal 
tail panels 

Table 4.12.1.2 

0.00221 

- 

>V 

0 

Zero-lift drag of isolated' exposed vertical 
tall panels 

Table 4.12.1.2 

- 

0.00111 

Hi- 

Number of junctures of tail with fuselage 

Figure 2.1 

2 

1 

t/c 

Thickness ratio of tail at juncture 

Table 4.1.2 

0.10 

0.09 

c 

Root chord of exposed tail, surface, m (ft) 

Table 2.1.1 

0.871 (2.858) 

,1.575 (5.167) 

S 

w 

Reference wing area, (ft^) 

Table 2.1.1 

14.40 (155.0) 

14.40 (155.0) 


Interference drag of horizontal tail 
surface due to fuselage 

Eq,. (4.12.3.4) 

0.00003 

- 


Interference drag of vertical tail surface 
due to fuselage 

Eq. (4.12.3.4) 

- 

0.00001 

o 

Net zero-lift drag of horizontal tail la 
presence of fuselage 

Eq. (4.12.3.6) 

0.00224 

- 

o 

Net zero-lift drag of vertical tail in 
presence of fuselage . 

Eq. (4.12.3.7) 


0.00112 


Table 4.12.3.3: Net zero-lift drag of nacelles in presence of wing 


Symbol 

Description 

Reference 

Magnitude 

^n 

0 

Zero-lift drag of one isolated nacelle 

Table 4.12,2.1 

0.00614 

<S)w 

o 

Zero-lift drag of isolated exposed wing panels 

Table 4.12.1,2 

0.00970 

w n 

Wing area overlapped by one nacelle, ta^ (ft^) 

Figure 2.1.1 

1.157 (12.45) 


Reference wing area, (ft^) 

Table 2.1.1 

14.40 C155.0) 

^n(w) 

o 

Net zero- lift drag of nacelles in presence 
of wing 

Eq. (4.12.3,8) 

0.01384 


165 








Table 4.12.4.1: Drag of wing and horizontal tail due to lift 


Symbol 

Description 

Reference 

Wing 

Horizontal 

Tail 

A 

Aspect ratio 

Table 2.1,1, 

10.32 

4.75 

6 

Twist angle, deg 

Table 2.1 

-3 

0. 

a 

Airfoil section lift curve slope, 
deg“^ 

Table 4.12 

0.115 

- 

a 

Lift-curve slope of lifting surface, 
rad-1 

Table 4.2.1 

5.090 

3.878 

X 

Taper ratio 

Table 2.1.1 ' 

0.5 

- 

M 

Mach number 

- 

0.081 

0.081 

B 

V 1 - 

- 

0.9967 

0.9967 

BA 

- 

- 

10.29 

4.73 

le 

Leading edge sweep angle, deg 

Table 2.1.1 

3.67 

0, 

^c/4 

Quarter chord sweep angle, deg 

Table 2.1.1 

1.835 

0 

(">‘) 

- 

- 

1.841 

0 

V 1 

Induced drag factor due to twist 

Figure 4.12.4.1 

0-00061 

1 


Bw 

Zero-lift drag factor due to twist 

Figure 4.12.4.1 

! 0.00221 

- 

w 

Bw/B 

- 

0.00222 

- 

®Re^MAC 

Reynolds lumber based on c 

w 

- 

2.3x10® 

2.3x10® 

LER. 

j Leading-edge-radius of airfoil, 
j ratio of chord 

- 

0.06 

0.1102 

c/c 

w 

Ratio of mean aerodynamic chords 

Table 2.1.1 

1. 

0.711 

La-H^cos^A, 

(N ) ^ — 

'•Re'' LER tanA- 

le - 

- 

- 

2.14x10® 

CB 

AX/cosA^^ 

- 

- 

5.171 

4.75 

R 

Leading-edge-suction parameter 

Figure 4.12.4.3 

0.961 

0.959 - 

a 

Span-efficiency factor 

Eq. (4.12.4.2) 

0.9095 

0.9850 


>w = ® 2087 " "'000210 C + 0.000264 based on 

i w a o 


= 

1 


*c ^ 
hi 

14.6987 \ 


Summary: 
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Table 4,12.4.2: Summary induced drag of wing and horizontal • tail 


a, deg 

a “ a + i 
w w 

deg 

. W 

Figure 4 . 2.4 . 1 


ei,> deg 

Figure 4. 9, 3.1 

a. => o - E, 
h h 

deg 

Figure 4. 2. 4,1 

Eq. (4.12.4.5) 

<'D,>h 

based on S 
> w 

-4 

-3,5 

-0.0533 

0.00037 

-0.26 

-3.74 

-0.2531 

0.00436 

0.00109 

-2 

-1.5 

0.1243 

0.00075 

0.56 

-2.56 

-0.1733 

0,00204 

0.00051 

0 

0.5 

0.3019 

0.00329 

1.38 

-1.38 

-0.0934 

0.00059 

0.00015 

2 

2.5 

0,4795 

0.00796 

2.19 

-0.19 

-0.0129 

0.00001 

“ 0 

4 

4,5 

0.6571 

0.01476 

2.99 

1.01 

0.0634 

0.00032 

0.00008 

15.9 

16.4 



1.494 

0.07565 

7.24 

■ 

8.66 

0.5861 

0.02337 

0.00583 


Table 4.12.5.1: Drag due to lift of fuselage and nacelles 


a, deg 


Eq. (4.12.5.3) 

% 

table 4. 4. 3. 3 

("Oi>n 

•Eq. (4.12.5.4) 

-4 

-0.01503 

0.00171 

-0.00780 

0.00055 

-2 

-0.01077 

0.00085 

-0.00398 

0.00014 

0 

-0.00619 

0.00027 ■ 

0 

0 

2 

-0.00128 

0.00001 

0.00414 

0.00015 

4 

0.00395 

0.00010 

0.00 844 

0.00059 


15.9 


O.OU71 


0.00975 


0.03737 


0.01037 






























Table 4.12.8,1; Drag of the complete ATLIT airplane 


0, deg 

Figure 4. 4. 3.1 

S 

o 

Eq. (4.12.3.1) 

Table 4,12.4.2 

■ <\>h 

Table 1.12.4.2 

i 

Table 4.12.5.1 

3. 

Table 4.12.5.1 

S 

Eq. (4,12.8.2) 

-4 

-0.1442 

0.03352 

0,00037 

0.00109 

0.00171 

0.00055 

0.03724 

-2 

0.0688 

0.03352 

0,00075 

0,00051 

0,00085 

0.00014 

0.03577 

0 

0,2822 

0.03352 

0,00329 

0.00015 

0.00027 

0 

0.03723 

2 

0.4963 

0.03352 

0.00796 

0 

0.00001 

0.00015 

0.04164 

4 

0.7111 

0.03352 

0.01476 

0.00008 

0.00010 

0.00059 

0,04906 

15.9 

1.7264 

0.03352 

0,07565 

0.00583 

0.00975 

0.01037 

0.13512 





Figure 4,12,1,1; Incompressible skin friction coefficient of a flat 
plate (Reference 3) 







Figure 4.12,3.1: Wing-body interference correlation factor 
(Reference 3) 
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Figure 4,12«4.1; Lift dependent drag factor due to linear twist 
(Reference 4) 


172 









Figure 4.1 '2. 4.1: Continued 
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. L.E, SUCTION . 
FULL L.E. SUCTION 

R *6 


AX 
cos A, 



1(b) 

j[| 

P 

B 



■ 

IQ 

01 

B 

B^ 

n 

■ 

01 

■ 

■ 

1 

LER 
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le 

> 1.3X10^ 


8 10 


(N ) cotA. 
LER 




2 2 
•M COS' A 


AX 


le 


cos A 


LE 


Figure 4.12.4.3: Leading-edge suction parameter at subsonic speeds 
(Reference 4) 





Figure 4.12.7.1: Full-scale wind tunnel data of increment in airplane 
drag due to open inlets and open cowl flaps 
3.5 million) 
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Figure 4,12.8.1; Comparison of predicted airplane drag with full- 
scale wind tunnel data (propellers removed, no 
engine cooling drag, 2.3 million) ■ 
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Figure 4.12.8.2: 


Comparison of predicted airplane drag with full- 
scale wind tunnel data (propellers removed, no 


engine cooling drag, N 


Re 


2.3 million) 

TT rrV OF 
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4.13 Effect of Horizontal Tail and Tab Deflection on Lift and 


Pitching Moments 


The contributions of the horizontal tail to the lift and pitching 
moment were considered in Section 4.10 and 4.11, respectively, on the 
basis of a fixed tail at zero Incidence setting- In this section the 
tail is considered as an all-moving surface with a geared tabi; 

4.13.1 Lift of the Horizontal- Tail in the Linear Lift Range 

The lift of an all-moving horizontal tail equipped with a tab is 
attributed to three superimposed sources: 

(1) Lift due to angle of attack of the tail with the tail at zero 
incidence, 

(2) Lift due to stabilizer deflection, i^, from the zero incidence 
position. 

(3) Lift due to tab deflection. 

These three sources can also be found in the following expression which 
describes the tail lift in the presence of the fuselage; including 
carryover effects onto the fuselage, referenced to the reference wing 
area: 


"^(hf) 


(C. ) 


(a - 


■“a h(hf) 








tab=0 







C4. 13.1.1) 


Lift due to angle of attack of the tail with the tail at zero incidence 
This lift can be represented as follows: 
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0 ^ =0 ^ 0 ° 
•tab 


(4.13.1.2) 


where 


(Cl horizontal tail lift curve slope in the presence 

of the fuselage, including the lift on the fuselage due to lift carry- 
over of the tail onto the fuselage, based on S^* 

(a - e^) represents the angle of attack of the horizontal tail. 

The calculations for the ATLIT airplane are summarized in Table 
4.13.1.1 and the result is (referenced to S^) : 


tab 


(4.13.1.3) 


where 

e^ can be obtained from Section 4.9 as a function of the angle of 
attack. 

Lift due to stabilizer deflection from zero incidence position: 

The lift due to the stabilizer deflection, i^, with ^ degrees 

can be obtained in the same manner as the lift due to angle of attack 
of the tail in Section 4.10 or: 

>h <\(f) + ‘=£(h)> S. (4.13.1.4) 

n ae q w 

■*■00 

where 

(C ) 

L^ h^ is the lift curve slope of the exposed tail panels based 

on S, . 
h 

e 

is the area of the exposed tail panels, 
e 

k^^^^ is the ratio, due to the stabilizer deflection, i^j of the 
lift on the stabilizer in the presence of the fuselage to the stabilizer 
alone obtained from Figure 4.13.1.1. 


f(h) 


is the ratio, due to the stabilizer-deflection, i^, of the 
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stabilizer lift carryover onto the fuselage to the fuselage alone, 
obtained from Figure 4.13.1.1. 

Applied to the ATLIT airplane and referenced to the reference wing 


area-, S , Table 4,13.1.2 shows that: 
w 


- 0-0156 1 ^ 


(4.13.1.5) 


where i^ represents the stabilizer deflection in degrees. <■- 

Lift due to tab deflection: In Reference 3 the following expression is 

given for the lift contribution of the tab: 


= c. 


tab 


b(f) 

tab a h 


(a 




) 

tab C, 


(a. 


) 

^tab c 




(4.13.1.6) 


I 


where 


(C.^ ) is the lift-curve slope of the horizontal tail surface 
a h'(f) 

alone in the presence of the fuselage based on 


^ \ \(f) S 

a h(f) ah h 

e 


(4.13.1.7) 


where 


(C ) is the lift-curve slope of the exposed horizontal tail 
a h 

panels based on . 

e 

is the ratio, due to the stabilizer deflection, i^, of the 

lift on the stabilizer in the presence of the fuselage to the stabilizer 

alone, obtained from Figure 4.13.1.1. 

Cc ) is the section lift-curve slope of the untabbed tail (5^ 
a h. 

Q degrees) obtained from Section 4.1. 

4.13.1.2 as a'^function of the aspect ratio of the horizontal tail, A^^, 


tab 


.. ) / (a. ) is the tab-chord factor obtained from Figure 

tab G °tab c. 

JLj f* m ^ f 


and the parameter, Cct. ) , which can be obtained from the insert in 

°tab c 
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Figure 4.13.1.2 as a function of When (a^ varies along 

“tab 

the span, an average value of (a^ ) , based on an average value of 

■^tab 

Cfab/'^j may be used. 

is the tab-span factor, obtained from Figure 4.13.1.2 as a 

function of the horizontal tail taper ratio, X^, and the span ratio, 

Ap, as defined in Figure 4.13.1.3. 

c„ is the section lift effectiveness of the tab and may be 
'tab 

obtained as follows; 


^theory ^ 

\ab 


^ ^theory 

\ab 


where 


8, is equal to Vi - and 

(c„ ) , is the theoretical lift effectiveness of the tab, 

A- theory 

“tab 

obtained from Figure 4.13.1.4 as a function of thickness 

ratio of the horizontal tail, (t/c), 

h 

c„ /(c„ is an empirical correction factor based on 

Z^ Z~ theory 

tab tab 

experimental data obtained from Figure 4.. 13. 1.5 as a function of 


“tab'"!. ^theory 


(c„ is the section lift curve of the untabbed tail: 

^ theory 

>theory ' ifo ' *te> J 

a 

K’ is an empirical correction for lift effectiveness of the tab 
at large deflections, obtained from Figure 4.13.1,6 

In Tables 4.13.1.3 and- 4.13.1.4_^the calculations are summarized 
and the lift effectiveness of the tab, referenced to the horizontal 


tail area, S^, is; 
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= 0.0257 per degree for = 6j 0, -6 and -12 degrees 

6, 


tab 


(4,. 13.1. 10) 


= 0.0239 per degree for ~ ~18 degrees 


tab 

Applied to the ATLIT airplane and referenced to the wing area, S^: 

^ ’h, =h- 


tab 0 ^ , n q w 

tab 


(4.13.1.11) 


= 0.3745 C. 


tab 


This in the case of a tab-to-stabilizer gear ratio, ^tab'^\’ 
and a dynamic pressure ratio at the tail, q^/q^, of 1.0. 

In Table 4.13.1.5 the tail lift in the presence of the fuselage, 
referenced to the wing area, is presented as a function of the stabi- 
lizer deflection, i^, and the angle of attack of the horizontal tail. 


“h = 


\(hf)^%hf)\ 


(4.13.1.12) 


T,= 0 

tab 


4.13.2 Maximum Lift of the Horizontal Tail 

The increment" of the maximtim lift coefficient due to trailing 

edge flaps can be determined with the following expression (based on 

the horizontal tail area, S, ) : 

» h 


fAn'\ - (K \ ^^h^tab ^ 

^ ^'5 ^tab S \ 

max tab max h 


(4.13.2.1) 


where 

^^h^tab horizontal tail area in front of and including 

the tab 

is an empirically derived correction factor to account for the 
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effects of wing planform, obtained from Figure 4.13.2.1 


(A ^tab increment in airfoil maximum lift 

^max 

coefficient due to the tab '"-and can be^ calculated as follows; 


^tab ^1 ^2 S *^5, ^base 
max max 


(4.13.2.2) 


where 

(A c^ ^base section maximum lift increment for 25-percent 

max 

chord flaps at a reference flap deflection angle of 60 degrees for plain 
flaps or tabs, obtained from Figure 4.13.2.2. 

is the factor accounting for other than 0.25, obtained 

from Figure 4.13.2.3 

is the factor accounting for a tab deflection angle other than 
60 degrees, obtained from Figure 4.13.2.4 

k^ is the factor accounting for tab motion as a function of 

^tab^^^tab^ reference’ ^ 

The maximum lift coefficient for the horizontal tail may now 

be determined as follows; 


(C, 


(4.13.2.3) 


where 

[ (Ct is the maximum lift coefficient of the 

L h(hr) 1 . = 0 
max tab 

untabbed tail based on the horizontal tail area and obtained from 
Section 4.10 

The summary calculations for the maximum lift coefficient are 
presented in Tables 4.13.2.1 and 4.13.2,2. 


4.13.3 Lift Curves of the Horizontal Tail through Stall 

Because of net lift and pitching moment coefficients of an airplane 
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for different stabilizer (elevator) positions are dependent upon the 
tail lift characteristics and could involve the stall region of the 
tail, operational tail lift curves for the ATLIT airplane are plotted 
in Figure 4.13.3.1. 

The following procedure was used in constructing the lift curve 
plots; 

(1) Using the information from Table 4.13.1.5, draw the slope of 

the basic lift curve (i^ = 0) up to the limit of linearity can be 

obtained from Table 4.13.1.1). 

(2) Spot the stall point for i^ = 0 degrees using the information of 
Table 4.13.2.2. Fair a curve similar to the fairing for the isolated 
tail in Figure 4. 2. 4.1 from the limit of linearity through the stall 
point. 

(3) Using the information of Table 4.13.1.5, draw the slope of the 
other lift curves (i^ = 4, -4, -8 and -12 degrees) parallel to the 
basic lift curve. 


(4) Using the maximum lift coefficients,, (C,. ). ,. of Table 

t hi hr) 

max 

4.13.2.2, draw horizontal lines to denote the different maximum 


lift values. 


(5) Make a plot, to be used as an underlay in tracing of the non- 
linear part of the basic lift curve. Translate this underlay plot 
relative to the basic lift curve of the selected stabilizer setting 
and the corresponding maximum lift coefficient. 


4.13.4 Lift and Pitching Moment Curves Including the Effect of 
Stabilizer Position 


The lift and pitching moment characteristics of the ATLIT airplane 
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may now be determined as a function of angle of attack and stabilizer 
deflection as follows ; 


where 




^ \fn \(hf) ^w 


C = C + 
m m ^ - 

wfn c 




w 


\(hf) ®w 


(4.13.4.1) 


(4.13.4.2) 


C is the tail-off lift coefficient which can be obtained 
wfn 

from Section 4.4. 

C is the tail-off pitching moment coefficient obtained from 
wfn 

Section 4.8 


C is the lift coefficient of the horizontal tail, referenced 

^(hf) 

to the tail area, with tail-fuselage interaction effects, angle of 
attack, stabilizer deflection and tab deflection accounted for 

(x^g - Xj^)/c^ is the distance from the center of gravity of the 
airplane to the quarter chord of the horizontal tail mean aerodynamic 
chord in chord lengths of the wing mean aerodynamic chord 

The lift coefficient as a function of the angle of attack, a, 
and the stabilizer deflection, i^, is calculated in. Table 4.13.4.1. 

The pitching moment coefficient calculations are shown in Table 4.13.4.2. 
The results are compared with the full-scale wind tunnel data in 
Figures 4.13.4.1 through 4.13.4.3. 

However, no wind tunnel data were available with the ATLIT 
airplane in the "fully clean" configuration. The effect of the 
stabilizer deflection was obtained- in the same manner as the effect 
of the horizontal tail. (See Appendix D.) The incremental lift 
and pitching moment due to the stabilizer deflection are obtained 
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for the "power-off" condition for the airplane "as built" and then- 
added to the results of the ATLIT airplane in the "fully clean" con- 
figuration (propellers removed) and a stabilizer setting of zero. 

The predicted lift and pitching moment curves show poor agreement 
with the experimental data, especially in the cases ij^ = -8 and -12 
degrees. At low angles of attack, the experimental results reflect' 
the tail stall when i^ = —8 and —12 degrees . The predicted curves , 
however, do. not, show this effect. 

In Figure 3.13.4.4 the calculated incremental lift due to the 
stabilizer deflection has been added up to the experimental lift curve 
with a stabilizer setting of zero. This graph gi-ves. a better indi- 
cation of the accuracy of the calculated stabilizer effectiveness . 

The comparison with the experimental results shows that the calculated 

lift effectiveness, , is higher than indicated by the wind tunnel 

^h 

data. This has also its effect on the pitch control effectiveness, 

, which is also higher than the experimental value as shown in 

Figure 4.14.4.2. 

According to Reference 3, the factor j the ratio of the lift 

carry-over onto the fuselage to the lift of the stabilizer alone, 

causes the discrepancy. The value of this ratio, k-.,.= 0.10, is too 

high for the type of tail-body configuration of the ATLIT airplane. 

In Reference 3 it is suggested that be considered negligible for 

this type of tail-body configurations. By neglecting this factor, 

the calculated tail effectiveness, C , decreases by approximately 

6 percent. ^ 

A second factor, which affects the change in lift due to stabilizer 

deflection, is the dynamic pressure ratio at the horizontal tail, q, /q , 

h 
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This ratio has been assumed to be equal to one. However, for the 

type of tail-body configuration of the ATLIT a value q^/q^=0.85 

seems to be more accurate (see Section 4.9). The decrease in d 3 raamic 

pressure ratio will decrease the calculated lift effectiveness, 

C, , and pitch-control effectiveness, C 
L-t — m-: — 

The above two changes will result in an improved a.greement between 

the predicted results and the full-scale wind tunnel data. 
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Table 4.13.1.1; 


Lift contribution of the horizontal tail due^ to 
angle of attack (i^ = <5^^^ = 0°) 


Symbol 

Description 

Reference 

Magnitude 


Lift curve slope of tall with tail-fuselage. 
Intersection effects, referenced to S^, deg 

Table 4.10.3.1 

0.0709- 


Limit of linearity, deg. 

Table 4.10.3.1 

12.2 


Average downwash at the horizontal tail, deg. 

Figure 4. 9. 3.1 

Variable 


Dynamic-pressure ratio at tall 

Figure 4. 9. 3. 2 

1.0 


2 2 

Area of horizontal tall, m Cft ) 

Table 2.1.1 

3.60 (38.7) 

s 

w 

2 2 

Reference wing area, m- (ft ) 

Table 2.1.1 

14,. 40 (155.0) 

%hnV“ ■■ ' 

- 



“ 0.0177 (a - 8^) based on and up 

12.2 degrees 
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Table 4.13.1.2: 


Lift contribution due to stabilizer deflection 
"tab - °°> ' 


Symbol 

Description 

Reference 

Magnitude 




Fuselage width at horizontal tail, mCft) 

Figure 2.1.2 

0.405 (1.33) 




Span of horizontal tail, aCft) 

Table 2.1.1 

4.133 (13,56) 

<VMh 



— 

— 

0.098 

^(f) 



Ratio of lift on movable tail in presence of 
body to tall alone 

Figure 4.13.1.1 

0.97 

’"f(h) 



Ratio of movable tail lift carryover on 
body to tail alone 

Figure 4.13.1.1 

0.10 

«t'h 
a e 



Lift curve slope of exposed horizontal ra->i 
panels referenced to Sj^ , deg“^ 
a 

Table 4. 2. 4.1 

. 

' 0.0649 

e 



Area of exposed horizontal tail panels 

Figure 2.1.1 

3,25 (34.9) 




Reference, wing area 

Table 2.1,1 

14.40 (155.0) 




Dynamic-pressure ratio at horizontal ra-tT 

Figure 4. 9. 3.2 

1.0 

Summary: 

h 







“ 0.0156 based on S 

n V 
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Table 4.13.1.3; Tab effectiveness of the ATLIT- 


Descriptlon 


Magnitude 


Mach noabec 


Thickness ratio of horizontal tail 


Table 4.1.2 


Trailing edge angle of horizontal tail, deg. 


Table 4.1.2 


Aspect ratio of horizontal tail 


Table 2.1.1 


Taper ratio of horizontal tall 


Table 2.1.1 


2 2 

Area of horizontal call, a (ft ) 


Table 2.1.1 


3.60 (38.7) 


2 2 

Area of exposed horizontal tall panels, m (ft ) 


Table 2.1.1 


3.23 (34.9) 


• 2 2 
Reference wing area, m (ft ) 


Table 2.1.1 


14.40 (155.0) 


Ratio of tab chord to tail chord 


Distance from root chord of tail to inboard 
edge of tab as a fraction of call semi-span 

Distance from root chord of tail to outboard 
edge of tab as a fraction of tail semi-span 


Figure 2.1.2 


Figure 2.1.2 


Section lift— curve slope of tail, deg 


Table 4.1.2 


a e 


Lift-curve slope of exposed call panels, based 
on 

e 

Ratio of lift on. movable tail In presence of 
body to tall alone 


Table 4. 2.4.1 


Table 4.13.1.2 0.97 


^*^1 ^ h ( f ) 


Lift curve slope of tail surface alone in 
presence of fuselage, based on deg'^ 


Eq. (4.13.1.7) 0.0563 


*tab Z 


Section lift parameter 


Figure 4.13.1.2 -0.556 


“tab 



Ratio of finite and section lift parameter 


Span factor 

Theoretical section lift-cutve slope of tail, 
deg”^ 


Figure 4.13.1.2 j 1.071 


Figure 4,13.1.3 0.83 

Eq, (4.13.1.9) 0.1182 


^theory 


/ , Function of. c . /c. and (c ). /(c, )„. 

>thcory ^ ^ 

“tab 


Figure 4.13.1.5 0.873 


(c > -1 

1. theory Theoretical section effectiveness, deg 

^tab 


Figure 4.13.1.4 0.0634 
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Table 4.13.1.4s Lift contribution due to tab deflection 


\ab 

K' 

Figure 4.13.1.6 

C, based on S. , deg ^ 

L. n 

Cab 

Fq, (4.13.1.6) 

6 

1.0 

0.0257 

0 

1,0 

0,0257 

-6 

1.0 

0,0257 

-12 

1.0 

0.0257 

-18 

0.93 

0,0239 


Table 4.13.1.5; 


Lift contribution of the horizontal tail with 
^tab^^h " 


1,^. deg. 

^ab> ‘‘^8- 

„ p^^l, =0 

li{h£) h 0 
Cab 

Eq. (4.13,1.3) 

Eq, (4,13.1.5) 

Cab 

Eq. (4.13.1,11) 

c, 

^h(hf) 
based on 
Eq. (4.13.1.12) 

hi(Uf) 

based on S. 

n 

4 

6 







0 

0 







-4 

-6 


0.0177 (a - Cj^) 

0.01561. 

n 

0.00961j^ 

0.0177 (a - Cjj) 

0.0709(a- Ejj) 

-8 

-12 





+ 0.02521. 

n 

+ 0.10091, 
a 

-12 

-18 

0.0177(a - Cj^) 

0.01561j^ 

0.00901j^ 

0.0177(« - Cjj) 

0.0709 (a - Ejj) 






+ 0.02461. 

a 

4- 0.09851. 

h 





















Table 4.13-2-1: Change in. max^imum lift, cbeffic-leat-of the— hor-izorital 
tall due- to- tab deflection.- 


Symbol 

Description 

Reference 

Magnitude 

<^c/4>h 

Sweep of horizontal tail along c/4 line, deg. 

Table 2.1.1 

0 

(t/c)h 

Section thickness ratio of horizontal tail 

Table 4.1.2 

0.10 


Tab chord as ratio of tail chord 

Figure 2.1.2 

0.204 

^®h^tab 

Area of horizontal t^l in front of and 
Including tab, ni‘^(ft^) 

Figure 2.1.2 

3.40 (36.6) 

Sh 

2 2 

Area of horizoucal tail, m (ft ) 

Table 2.1.1 

3.60 (38.7) 

5 

w 

2 2 

Reference wing area, m (ft ) 

Table 2.1.1 

14.40 (155.0) 


Correction factor for wing planform 

Figure 4.13.2.1 

0.92 

^base 

may 

Section ma-xiTmim lift increment for 25-percent 
chord tab 

Figure 4.13.2.2 

0.81 


Factor accounting for other than 0.25 

Figure 4.13.2.3 

0.94 

^2 

Factor accounting for tab deflection other than 
60 degrees 

Figure 4.13.2.4 

Variable 

k 

Factor accounting for tab motion as a function 

of 5^ , /(6^ . ) - 

tab tab reference 

— 

1.0 


Horizontal^tail maximum lift coefficient with 
0 ^tab” presence of fuselage based on 

Table 4.10.3.2 

+ 1.214 

^h(hf)\ 
max ca 

Horizontal tail angle of attack, deg. 

3=0 



Table 4,10.3.2 

± 19-3 

summary: C^C > = k^ k 3 (Ac^ ^ 

max tao max h 




=» 0.6625 k- based on S, 
c n 
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Table 4.13.2.2: Maximum lift coefficient of the horizontal tail 


— 

S’ 

Sab* 

Figure 4.13.2.4 

max cab 
Eq. (4.13.2.1) 

max 
Ea. (4 

S(hf) 

.13,2.3) 

(a) 

(b) 

4 

6 

0.22 

0.146 

-1.068 

1.360 

0 

0 

0 

0 

-1.214 

1.214 

-4- 

-6 

-0.22 

-0.146 

-1.360 

1-068 

-8 

-12 

-0.40 

-0.265 

-1.479 

0.949 

-12 

-18 

-0.53 

-0.325 

-1.566 

0.862 


(a) for negative 

(b) for positive ol 
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Table 4,13.4.1: Lift coefficient of the ATLIT airplane 


a, deg e,^, deg deg 

Figure 9. 3.1 " St, 


from Figure 4.13.4.1 based on 


-0.26 

0.364 

1,377 

2.192 

2.990 

7.238 


wfn 

Table 4. 4. 3. 3 


-3.74 0.1384 -0.2632 -0,6688 -1,0724 -1.448 -0.0780 

-2.564 0.2218 -0,1818 -0.5854 -0.9890 -1,363 0,1141 

-1.377 0.3060 -0.0976 -0,5012 -0.9048 -1.280 0.3066 

-0.192 0,3900 -0.0136 -0.4172 -0.8208 -1.1956 0,4997 

1.010 0.4752 0.0716 -0,3320 -0.7356 -1.1104 0.6932 

8.662 1.0177 0.6141 0.2105 -0.1931 -0.5679 1.5731 


C, with Eg. (4.13.4.1) baaed on S 

L w 

1 

1 

1 

1 

I 

-0.0434 

-0.1442 

-0.2450 

-0,3458 

-0.4395 

0.1695 

0,0687 

-0.0321 

-0.1328 

-0.2262 

0.3830 

0.2822 

0.1815 

0.0807 

-0.0130 

0,5971 

0.4963 

0.3955 

0.2948 

0.2012 

0.8118 

0.7111 

0.6103 

0.5095 

0.4160 

1.8272 

1.7264 

1.6257 

1.5249 

1.4313 


Table 4.13.4.2; Pitching moment of the ATLIT airplane (6^^^/ij^=1.5) 


a, deg 


-0.26 

0.564 

1.377 

2.192 

2.990 

7.238 



% < 
S(lif) 

X - X, 

eg h 

* 

• ) from 

tu. deg 

c 

w 

a-Sh 

Era 

ih = 0“ 

1 0-4® 

-3.74 

-0,5447 

1.0438 

2.6324 

-2.564 

-0,8730 

0.7156 

2.3041 

-1.377 

-1.2044 

0.3842 

1.9727 

-0.192 

-1,5350 

0.0535 

1.6421 

1.010 

-1.8704 

-0.2818 

1.3068 

8.662 

-4.0057 

-2.4171 

-0,8285 


C vith Eg. (4.13.4.2) baaed on S, 
Ift w 


-0.1256 

-0.0806 

-0.0365 

0.0067 

0,0490 

0.1674 


lh»0" 


V-8- 

0.1350 

0.5317 

0.9283 

0.0981 

0.4947 

0,8913 

0.0594 

0.4560 

0.8527 

0.0201 

0.4167 

0.8133 

-0,0214 

0.3753 

0.7719 

-0.4361 

-0.0395 

0.3572 


3.936 obtained from Table 4.11.1 



















































Figure 4.13.1.2: Tab chord factor with c_ = c , (Reference 3) 















Figure 4.13.1.5; Empirical correction for lift effectiveness of 

plain trailing edge flaps ; c, = c^' , (Reference 3) 

t tab 
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Figure 4,13.2.1: Correction factor for wing planfona (Reference 3) 
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/Split and plain 
flap 













0 2 4 6 8 10 12 14 16 18 20 

- Airfoil thickness, percent chord 


Figure 4.13.2.2: Section maximum lift increment for 25-percent chord 

tab at reference angle of 60 degrees (Reference 3) 



1 . 

1 . 

*^1 


0 4 8 12 16 20 24 28 32 

Flap chord, percent chord 

Figure 4.13.2.3; Factor accounting for tab chord other than 25-percent 

(Reference 3) 

1.0 

.8 
.6 
,4 
.2 

0 

Figure 4.13.2.4: Factor accounting for tab deflection other than 

60 degrees (Reference 3) 
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Figure 4.13.3.1: Calculated lift curves of the horizontal tail with 

tab geared to elevator, based on ” l*0j 

^tab^\ = \e " niillion) 
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Figure 4.13.4.1: Comparison of predicted lift coefficient with full- 

scale wind tunnel data (propellers removed, N = 

2 . 3 million) 
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Figure 4.13.4.2: Comparison of predicted pitching moments with 

experimental results (propellers removed, N 
2.3 million) • Re 
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Figure 4.13.4.3: Comparison of calculated pitching moments with full- 

scale wind tunnel data (propellers removed, N = 

2.3 million) 
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Fxgure 4.13.4.4: Comparison of predicted stabilizer effectiveness 

with experimental results (propellers removed, N = 
2.3 million) Fe 
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CHAPTER 5 


PREDICTION OF POWER-ON AERODYNAMIC CHARACTERISTICS 

The effects of power from propeller operation are generally sig- 
nificant on the stability and control characteristics of an airplane. 
Unfortunately, because the propeller slipstream usually interacts with 

r 

the flow around several airplane components, a number of separate effects 
must be accounted for. Although some of the effects have been accounted 
for by theoretical analysis, many are usually estimated by empirical 
methods . 

In the following sections, the effects of power on lift, pitching 
moments, and drag are considered by the methods presented, in Reference 
3. In Figures 5,1 and 5.2 the parameters are defined, which are nec- 
essary to discuss power effects. 


5.1 Propeller Power Effects on Lift 

The effects of the propeller on the lift forces acting on the air- 
plane may be divided into two groups, those due to the propeller forces 
and those due to the propeller slipstream. On this basis the lift of 
the airplane may be written as follows: 


‘=L-=L 

prop off p ^ p 




(5.1.1) 


where 

C is the lift coefficient of the airplane with the 

prop off 

propellers removed, obtained from Section 4.13. 
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(AC ) is the lift component of the thrust. 

(ACl)jj is the lift component of the propeller normal force. 

(AC ) . - is the change in lift due to the change in dynamic pres~ 
sure over the section of the wing in the propeller slipstream. 

(AC,. ) is the change in lift of the section of the wing in the 

li £ 

p . 

propeller slipstream due to a change in angle of attack from the pro- 
peller' downwash. 

(AC^ ).- is the change in lift contribution of the horizontal tail 

hi 

resulting from change in dynamic pressure at the tail due to power. 

(AC ) / * % is change in lift contribution of the horizon- 

n h' power 

tal tail resulting from change in doimwash at the tail due to power. 

In the following discussion of power effects on lift, the airplane 
will be considered initially on the basis of tail-off lift character- 
istics, followed by horizontal tail contributions to lift with tail- 
fuselage interaction effects included, or: 


C, = C„ + C. 


L L 


wfn 


hi(hf) 


(5.1.2) 


where 




+ (AC^).- + (AC^) 

L Aq^ 


(5.1.3) 


5.1.1 Propeller Power Effects in Linear Lift Region 


The above mentioned increments in lift can be determined by the 
following steps. 
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The lift component of the thrust, (AC^ ) , can be calculated as 
follows : 

a) Calculate the angle of attack of the thrust axis measured from the 
free stream direction by: 

Ut = i^ + a (5. 1.1.1) 

where 

i^ is the incidence of the thrust line referenced to the X-body 
axis (Figure 5.2). 

b) Calculate the thrust coefficient due to one propeller (if not given 
or assumed) by: 

T^’/prop = 

where the thrust par propeller is a given or chosen quantity. 

c) Calculate the lift component of the thrust as follows: 

(AC^)^ = n(T^Vprop) sina^ = sina^ (5.1.1. 3) 


Thrust/propeller 


q S 
- “>■ w 


(5.1.1. 2) 


where 

n represents the number of propellers 


The lift component of the propeller normal force, 

P 


obtained as follows: 

a) Calculate the empirical normal force factor of the propeller by: 




= 262 


b 

(-^) 

'^0.3R 


+ 262 


b 


+ 135 


b 

(-^) 

^R '^0,9R 


(5. 1.1. 4) 


where the subscripts indicates the radial position from the center of 
the propeller where the blade width, b , is to be evaluated. 
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b) From Figure 5. 1.1.1 obtain a value for the propeller normal force 


coefficient, [(C„ ) ]„ based on the propeller normal force factor, 

\ p 

= 80.7, as a function of propeller blade angle, and type of pro- 
peller. 

c) Calciilate the true propeller normal force coefficient, (C„ ) , per 

N p 
a ^ 

radian by; 




^Sl^^p ^^Sjp^Kjj=80.7 ^80.7 “ 


(5. 1.1. 5) 


d) Determine the propeller correlation parameter defined as; 


S^(T^Vprop) 


8 R 
P 


2 


(5. 1.1.6) 


e) From Figure 5. 1.1. 2 obtain a value for the propeller inflow factor, 
f, as a function of the propeller correlation parameter, which can be 
obtained from Equation (5. 1.1. 6) 

f) Calculate the propeller disc area by; 


(S /prop) = irR- 
P P 


(5.1.1. 7) 


g) Determine the upwash gradient at the propeller, -(3£^/3a), from 

Figure 5. 1.1. 3 as a function of the wing aspect ratio, A , and the para- 

w 

meter, where is defined in Figure 5.2 and in Figure 5.1 

h) Calculate the local angle of attack of the propeller plane as follows: 


a 

P 




3e 
u 

3a 


(a^ - a ) 
W o 

w 


(5. 1.1. 8) 
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where 


a is the angle of attack of the wing obtained from: 
w 


a = a + i 
w w 


(5.1.1. 9) 


i is the wing incidence, obtained from Section 4.4. 
w 


a 


is the wing zero-lift angle obtained from Section 4.2. 


w 


i) Calculate the lift component of the propeller normal force by; 


(ACl)jj n f (Cj^ ) 
p a 


ctp (S^/prop) 


cosa„ 


(5.1.1.10) 


w 


The increment of lift due to the change in angle of attack on the 

wing induced by the propeller flow field, (AC^) , can be obtained as 

L C 

p 

follows : 

a) Calculate the downwash gradient,8e^/3a-p,by: 


ae 


to - h + S ’p 

p a 


(5.1.1.11) 


where 


C^ and C^ are obtained from Figure 5. 1.1. 4. 

(C„ ) is the true propeller normal force coefficient obtained 
N p 
a 

from Equation (5. 1.1. 5). 


b) Calculate the downwash behind the propeller, by: 

ae 


g = -r-*- a 
p da p 

P 


(5.1.1.12) 


where 

Op is the angle of attack of the propeller plane, obtained from 
Equation (5. 1.1. 8). 

c) Calculate the change in angle of attack of the portion of the wing 
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as follows; 


immersed in the propeller slipstream, » 

i 


(Aa) 


S. 

1 


3e 


(5.1.1.13) 


1 - 


u 


8a 


d) Calculate the upwash angle at the propeller, -e^, by: 


8e 

-e = - (a - a ) 

u 8a w o 

w 


(5.1.1.14) 


e) Calculate the vertical distance, z , from the body X-axis to the 

s 

centerline of the propeller slipstream at the longitudinal station of 
the quarter-chord of the mean aerodynamic chord of the immersed portion 
of the wing (Figure 5.2) by: 


z 

s 


__E_ 

57.3 





(5.1.1.15) 


where 

z^ is the vertical distance from the X-body axis to the thrust 
axis (Figure 5.2). 

f) Calculate the span of the immersed wing, b^/prop, by: 

b./prop = 2\4 ^ - (z - z )^ (5.1.1.16) 

1 p s w 

where 

(z - z ) is defined in Figure 5.1. 
s w ‘ 

g) Calculate the immersed wing area, S^/prop, by; 

S^/prop = (b^Zprop)c^ (5.1.1.17) 


h) 


Calculate the aspect ratio. 



of the immersed portion of the wing 


by; 


A. 

1 


b./prop 

X 


c . 


X 


(5.1.1.18) 
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i) Calciilate the increase in dj^namic pressure due to the propeller 

I 

slipstream on the immersed portion of the wing as follows: 



S (T '/prop) 
w c ' 



(5.1.1.19) 


j) Calculate the change in wing lift due to change in angle of attack 
resulting from propeller downwash, s , hy: 




- n(l + 


^eo 


a prop off 3 


S^/prop 


w 


(5.1.1.20) 


where 

(C^ ) represents the lift curve slope of the wing with 

a prop off 

the propeller removed obtained from Section 4.2 


The contribution of power to lift due to change in dynamic pressure 
on the immersed portion of the wing is obtained as follows: 

a) From Figure 5, 1.1. 5 obtain a value for the empirical constant 

as a function of the wing aspect ratio, A^, the immersed-wing aspect 
ratio, A^, and the propeller correlation parameter obtained from Equa- 
tion (5.1.1. 6) 

b) Calculate the increment of lift due to change in dynamic pressure 
on the immersed portion of the wing by: 




n K, 



w 


prop off 


(S^/prop) 

g 

w 


(5.1.1.21) 


where- 

(C_ ) is the lift coefficient of the wing obtained from 

L prop off " 

Section 4.2 as function of the wing angle of attack, a . 
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The lift of the airplane without horizontal- tail can be obtained 


with Equation (5.1.3) where (C^ ^prop off tail-off and pro-- 


wfn 


pellers removed lift coefficient obtained from Section 4,4,‘ 


The contribution of the horizontal tail to the lift of the airplane 
may be obtained as follows: 

a) From Figures 5. 1.1. 6 obtain a. value for the downwash increment due 

to power, ^ function of the propeller correlation para- 

meter, the downwash at the tail. T'rfrth- propeller (s) xemovediT.'(ej^)p^Q.^' o£f» 

obtained from Section' 4.9 and the .parameter z, /2R where z, is defined 

■ ^ - - .. ^T P ^T 

in Figure- 5,2, 

b) Calculate the "effective distance, z , parallel to the z-body axis, 

eff 

from the quarter chord of the horizontal-tail mean aerodynamic chord 
to the centerline of the propeller slipstream by; 


n J.J. 
eff 


= z - 


57.3 ^u ^p ^^h^prop off ^^^h^ power] - z^ 


where 


(5. 1.1. 2: 


is the distance from the quarter chord of the mean aerodynamic 
chord of the immersed portion of the wing to the quarter chord of the 
horizontal tail mean aerodynamic chord (see Figure 5.2), 

is the distance, parallel to the z-body axis, from the X-body 
axis to the quarter chord of the horizontal tail mean aerodynamic chord, 
obtained from Figure 5,2, 

c) From Figure 5.1,1. 7 obtain a value for the power-induced increment 

in dynamic-pressure ratio at the tail, Aq /q , as a function of the 

- ri ^ 

propeller correlation parameter, the ratio of horizontal tail immersed 

area and horizontal tail area, S. /S, , and the parameter z, /R . When 

n. n . n -- p 

X eff ^ 
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the increment obtained is negative, it can be assumed to be zero 
d) Calculate the angle of attack of the horizontal tail by: 


a, = a - (e, ) - (As, ) 

h h prop off h power 


(5.1.1.23) 


c) From Section 4.13 obtain a- value for the lift, (C„ t 

\(M) 

of ’the horizontal tail referenced to the horizontal tail area, S, , and a 
dynamic pressure ratio, q^/ q^, equal to one. 

f) Calculate the contribution of the horizontal tail to the lift of 
the airplane by: 

\(hf) \(hf) ^w 


(5.1.1.24) 


where 


(q, /q ) __ follows from Section 4.9, 

h «> prop off 


The total lift of the airplane can be obtained with Equation (5.1.2). 
5.1.2 Propeller Power Effects on Maximum Lift 

Over the linear lift-curve slope, increases in lift due to propeller 
power result from the factors discussed above. However, near or at max- 
imum lift an additional increase in lift occurs because the angle- of 
attack for stall increases with power. This effect depends primarily 
upon the ratio of the immersed wing area to the total area, S^/S^. The 
propeller power effect on maximiim lift can be obtained as follows: 
a) From Figure 5. 1.2.1 obtain a value for the empirical constant, K, 
as a. function of the ratio of immersed wing area to the total wing area. 
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b) Calculate the increment in taii-off lift due to power, (AC_) , 

^ * L power 

at the propeller off angle of attack for maximum lift: 




c) Calculate the increment of maximum lift due to propeller power by; 

AC, = K(AC,)’ (5. 1.2. 2) 

L L power 

max' 


The complete power-on tail-off lift curve is constructed as follows; 

a) Plot the linear portion of the power-on lift curve. 

b) Draw a horizontal line representing the power-on tail-off maximum 
lift coefficient. 

c) Translate the non-linear propeller-off portion of the lift curve to 

a tangency with a) and b) . This construction not only shapes the power- , 
on lift curve but also fixes the power-on stall angle. The construction 
of the power-on tail-off lift curve is also demonstrated in Eigure 5. 1.2. 2, 
The power-on lift curve of the complete airplane can be obtained 
by adding the. horizontal tail lift to the tail-off lift curve as sho'wn 
by Equation (5.1.2). 

5.1.3 Lift Characteristics of the ATLIT Airplane 

In Tables 5. 1.3,1 through 5. 1.3. 5 the calculations are summarized 
which lead to the lift curve of the airplane including power effects. ' 

The results for two power settings are shown in Figures 5, 1.3. 2 and 5. 1.3. 3 
and compared with experimental data from Reference 2. 

The predictions have been performed for a Reynolds number of 2.3 
million. However-, no power-on wind tunnel data were available of the 
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ATLIT in the "fully clean" configuration at a Reynolds number of 2.3 
million. Data were available for a Reynolds number of 3.5 million. As 
is shown in Reference 2, a change in. Reynolds number of this magnitude 
has an almost negligible effect on the lift of the airplane. 

The vertical position of the center of gravity of the ATLIT air- 
plane, used in this study, is situated at a distance of 0.333 m- (13.1 in.) 
above the wing root chord. 

The predicted lift curves show fair agreement with the wind tunnel 
data. However, when the predicted increment in lift due to power is‘ 
added to the experimental lift curve obtained with propellers removed, 
cowl flaps open, and engine inlets open, good agreement is obtained with 
the wind tunnel data of Reference 2. The increment in lift due to power 
can be written as follows; 


(AC •) = 

L power L 


- rc "i 

^ L^prop off 


(5. 1.3.1) 


where 

is the predicted lift coefficient of the airplane including 
power effects, obtained from Table 5. 1.3. 4. 

is the predicted lift coefficient of the airplane 
with propellers removed, obtained from Section 4.10. 

The effect of engine cooling on the lift of the airplane is briefly 
discussed in Section 4.10. 
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Table 5.1. 3.1; Lift due to direct action of the propeller forces. 


Symbol 

DeecripClon 

Reference 

Magnitude 


Incidence of thrust line to reference Xj^-axls, deg 

Figure 5.2 

0 

“t 

Angle of attack of thrust axis, deg 

Eq. (5. 1.1.1) 

a 

T^/prop 

Thrust coefficient 

- 

Variable 

n 

Number of propellers 

Figure 5.1 

2 


Lift component of the thrust 

Eq.(5.1.1.3) 

2(T^/prop) slna 








Table 5. 1,3.1: Continued, 


Symbol 

Description 

Reference 

Magnitude 

b 

P 

Width of propeller blade, m (ft) 


.127 (.417) at ,3R 

P 

.157 (.515) at .6R 

.108 (.354) at ,9R 
P 

R 

P 

Propeller radius, m (ft 

Table 2,1 

.97 (3.17) 

S 

Propeller normal force factor 

Eq. (5.1. 1.4) 

92.22 

8' 

Propeller blade angle, deg 


Variable 


Propeller normal force parameter, rad"^ 

Figure 5. 1.1,1 

Variable 

>p 

a ^ 

Propeller normal force derivative, rad ^ 

Eq, (5. 1.1.5) 

l*^l‘^^^SlJp^Kjj=80,7 

s 

w 

S„(T'/prop) , 

Reference wing area, (fc2) 

Table 2.1.1 

14.40 (155.0) 

8R 2 
P 

Propeller correlation parameter 

Eq. (5.1. 1.6) 

1.93(T*/prop) 

f 

Propeller Inflow parameter 

Figure 5. 1,1,2 

Variable 

Sp/prop 

Propeller dice area, (ft:^) 

Eq. (5.1.1. 7) 

2.93 (31.57) 

A 

w 

Wing aspect ratio 

Table 2.1,1 

10.32 

X ’ 

Distance of' propeller forward of c^/4, m (ft) 

Figure 5.2 

1,95 (6,41) 

_ 

Mean aerodynamic chord of Immersed wing area, 
m (ft) 

Figure 5.1 

1,34 (4.38) 

-Oe^Sa) 

Dpwash gradient at propeller 

Figure 5. 1.1. 3 

.155 


Wing incidence- at root, deg 

Table 2.1 

0.5 

a 

w 

Wing angle of attack, deg 

Eq. (5. 1.1.9) 

a. -b 0.5 

a 

Wing zero-lift angle of attack relative to 
wing chord, deg 

Table 4, 2. 4,1 

-2.89 

°T 

Angle of attack of thrust line, deg 

Table 4. 1.3.1 

a 

ct 

p 

Angle of attack of propeller plane, deg 

Eq.(5.1.1,8) . 

1.155a + 0.525 
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Table 5 . 1 , 3 . 1 : Continued 


T^Vprop 

8 R 2 
P 

Eq. (5.1.1. 6) 

£ 

Figure 5. 1.1. 2 

6', deg 
Appendix 6 

Figure 5. 1.1.1 

a ^ 

Eq. (5.1.1.5) 

0 

0 

0 

18 

.090 

.10 

,0458 

. .0882- 

.1.075 

18 

.090 

.10 

.0985 . 

.1900 

. 1.154 

18 

.090 

.10 


Table 5.1. 3.1; Concluded- 


o, deg 

(dCj_)^; 

Table 5.1. '3.1 


P 

•Eq. (5.1, 

1.10) 


T ’ 
c. 

“p. deg 

T ' 
c 


0.0915 


Eq.(5.1.1.8) 

0 

0.0915 

0.1970 

-4 . 

0 

-0.0064 

-0.0137 

-4.09 

-0.0029 

-0.0031 

-0.0033 

-2 

0 

-0.0032 

-0.0069 

-1.78 

-0.0013 

-0.0014 

-O.OOIS 

0 

0 

0 

0 

0.53 

0.0004 

0.0004 

0.0004 

2 

0 

0.0032 

0.0069 

2.84 

0.0020 

0.0022 

0.0024 

4 

0 

0.0064 

0.0137 

5.15 

0.0036 

0.0039 

0.0042 

15.9 

0 

0.0251 

0.0540 

18.90 

0.0129 

0.0139 

0.0149 ■ 


230 




















Table 5.1.3.2: Wing lift increments due to propeller slipstream 

effects • 



Propeller toTrelation- parameter Eq. (5.1.1.6) 1.93 (T^’ /prop) 

Factor for determining propeller downwasti Figure 5. 1.1. 4 Variable 

Factor for determining propeller downwash Figure 5. 1.1, 4 Variable 

Propeller normal force factor Table 5. 1.3. 3 0.10 

Downwash gradient behind propeller- Eq. (5.1.1.11) + O.IOC2 

Angle of attack of propeller plane, deg Eq. (5.1. 1.8) 1.155a+ 0.525 

Propeller- downwash behind propeller, deg Eq. (5.1.1.12) Variable 

Change in angle of attack of immersed xrfng, Eq. (5.1.1.13) — s /I. 55 

deg P 

Wing upwash at propeller plane, deg Eq. (5.1.1.14) .0.155 (a + 3.39)- 

Distance of propeller forward of Cj^/4, m (ft) Figure- 5.2 1.95 (6.41) 

Distance from X-body axis to thrust line, Figure 5.2 -0.128 (-0.417) 

m (ft) 

Defined in Figure 5.2 Eq. (5.1.1.15) Variable 

Distance from X-axis to quarter chord im- Figure 5.1 0.100. (0.329) 

mersed wing mean aerodynamic shord, m (ft) ' • 

Immersed span per propeller Eq. (5.1.1.16) Variable 

Mean aerodynamic chord of immersed wing Figure 5.1 1.34 (4.38) 

area, m (ft) 

Immersed area per propeller, m^ (ft^) Eq. (S.l.l.'l?) Variable 

Aspect ratio of immersed portion of wing Eq. (5.1,1.18) Variable 

Change in dynamic pressure ratio on immersed Eq. (5.1.1.19) Variable 

wing 

Lift curve slope of wing, rad~^ (deg~^) Table 4. 2. 4.1 5.090 (0.0888) 

Correlation parameter for added lift due Figure 5. 1.1. 5 Variable 

to power 


Lift coefficient of wing 


Figure 4. 2. 4.1 


Variable 




Table 5. 1,-3. 2: Continued 


T ’ 
C 

8r“2 

P 

^=1 

Figure 4. 1.1. 4 

^2 

Figure 4, 1.1, 4 

3e 

_E' 

3a 

Eq. (5.1.1.11) 

5= 

Eq. (5.1.1.19) 

0 

0 

0 

.0,25 

0.025 

0 

0.0915 

0.0832 

0.0783 

0.25 

0.1033 

0.2246 

0.1970 

0.1900 

0.1304 

0.25 

0.1554 

0.4836 
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Table 5. 1.3. 2: Concluded 


a, deg 

Op. deg 

Table 5. 1.3.1 
1 

deg 

Eq. (5,1.1.14) 

<lcg 

Eq. (5.1.1.12) 

(6a)s^, deg; Eq. (5.1.1.13) 

fti 

Eq. (5.1.1.15) 

T * 

c 

T ' 
c. 

T ' 
c 

0 

0.0915 

0,1970 

0 

0.0915 

0.1970 


0.0915 

0.1970 

-4 

o 

-0.095 

-0.1023 

-0.4217 

-0.6348 

0.0886 

0.3651 

0.5496 

0.0080 

-0.006 

-0.0295 

-2 

-1.78 

0,216 

-0.0445 

-0.1835 

-0.2763 

0.0385 

0,1589 

0.2392 

-0.2224 

-0.2380 

-0.2481 

0 

0.53 

0.526 

0.0133 

0.0546 

0.0823 

-0,0115 

-0.0473 

-0.0713 

-0.4743 

-0.4697 

-0,4666 

2 

2.84 

0.836 

0.0710 

0.2928 

0.4400 

-0.0615 

-0.2535 

-0.3816 

-0.7263 

-0.7018 

-0.6852 

4 

5.15 

1.146 

0.12S8 

0.5310 

0.7993 

-0.1115 

-0.4597 

-0.6920 

-0.9782 

-0.9338 

-0.9038 

15.9 

1 

18,90 

2.990 

0.4725 

1.9486 

2.9333 

-0.4091 

-1.6871 

-2.5397 

-2.4773 

-2.3122 

-2.2020 


ft 


T ’ 
c 

0.0915 


bj/ptop,' ft; Eq, (5.1.1,16) 


Sj^/prop, fc j Eq, (5.1.1.17) 


T ’ 


(ACj^)^ , Eq. (5.1.1.20) 


0 


-0.3210 

-0.5514 

-0.8033 

-1.0553 

-1.3072 

-2.8063 


-0.3350 

-0.5670 

-0.7987 

-1,0308 

-1.2628 

-2.6412 


0.1970 


-0.3585 

-0.5771 

-0.7956 

-1.0142 

-1.2328 

-2,5310 


6.3074 

6.2434 

6.1331 

5.9784 

5.7759 

2.9486 


0.0915 


6.3045 

6.2378 

6.1355, 

5.9954 

5.8152 

3.5060 


0.1970 


6.2993 

6.2341 

6.1371 

6.0068 

5.8409 

3.8173 


27.626 

27.346 

26.863 

26.185 

25.298 

12.915 


0.0915 


27.614 

27.322 

26.873 

26,260 

25.471 

15.356 


0.1970 


0 


27.591 

27.305 

26.880 

26.310 

25.583 

16.720 


0.0028 

0.0012 

-0.0004 

-0.0018 

-0.0032 

-0,0061 


0.0915 


0.0141 

0.0061 

-0.0018 

-0.0093 

-0.0164 

-0.0364 


0.1970 


0.0258 

0.0111 

-0.0033 

-0.0171 

-0.0301 

-0.0722 




Eq. (5.1.1.18) 

K, j Flguce 5.1.1. 5 

' L ^prop off 

(AC )^- ; Eq. 

(5.1.1.21) 

1 


T ' 



T ' 




T * 






c 




c 


0 

0,0913 

097ti 

0 

0,0915 

0.1970 

Figure 

0 

0,0915 

0.1970 

- 

1.439 

1.43 9 

- 

1.0 

0.96 

-0.0533 

0 

-0.0043 

-0.0088 

- 

1.424 

1.423 

- 

1.0 

0.96 

0.1243 

0 

0.D093 

0.0203 

- 

1.401 

1.401 

- 

1,0 

0.96 

0.3019 

0 

0.0235 

0.0486 

- 

1.369 

1.371 

- 

1.0 

0.96 

0.4795 

0 

0,0365 

0.0756 

- 

1.328 

1.334 • 

- 

1.0 

0.96 

0.6571 

0 

0.0485 

• 0.1007 


0.800 

0.872 

- 

0.98 

0.94 

1.494 

0 

0.0652 

0.1465 






























Table 5.1, 3.3: Tail-off lift characteristics with power-on of the 

ATLIT airplane,, 


a. deg 

Figure 4. 4. 3.1 


(ACj^)^; Table 5. 1.3,1 


0.0915 


0.1970 



, ; Table 5. 1.3.1 
P 


T ' 
c 


-0.0780 

0.1141 

0.3066 

0.4997 

0,6932 

1,5731 


0 

0.0032 

0.0064' 

0.0251 


-0.0137 

-0,0069 

0 

0.0069 

0.0137 

0.0540 


-0.0029 

-0,0013 

0.0004 

0,0020 

0.0036 

0,0129 


-0.0031 

-0.0014 

0.0004 

0.0022 

0.0039 

0.0139 


0.0004 

0,0024 

0,0042 

0.0149 


(AC,) - ; Table 
L Aq 
‘w 

5. 1.3. 2 

T ' 
c 

0 

0.0915 

0.1970 

0 

-0.0043 

-0.0088 

0 

0.0098 

0,0203 

0 

0.0235 

0.0486 

0 

0.0365 

0.0756 

0 

0.0482 

0.1007 

0 

0.0652 

0,1465 




Table 5. 1.3. 2 


0.0023 

0.0012 

-0.0004 

-0.0018 

-0.0032 

-0.0061 


0.0915 


-0.0093 

-0.0164 

-0.0364 


0.1970 


0.0258 

0.0111 

-0.0033 

-0.0171 

-0.0301 

-0,0722 


0.6936. 

1.5799 


; Eq. (5.1.3) 


0.0915 


-0.0777 

0.1254 

0.3287 

0.5323 

0.7356 

1.6409 


0.1970 


-0.0780 

0.1371 

0.3523 

0.5675 

0.7817 

1.7163 







































Table 5. 1.3. 4: Power effects on horizontal tail lift and total lift 

of the ATLIT airplane • 


Symbo], 

Description 

Reference 

Magnitude 

s 

V 

2 2 

Reference tjing area, m (ft ) 

.Table 2.1.1 

14.40 (155.0) 


2 2 

Horizontal tall area,, m. (ft ) 

Table 2.1.1 

3.60 (33.7) 


Vertical distance from thrust line to 
horizontal tail, m(ft) 

Figure 5.2 

-0.127 (-0.417) 

R 

P 

Propeller radius, m(ft) 

Table 2.1 

0.97 (3.17) 

z. /2R 
br P 

— 

— 

-0.066 

^^h^prop off 

Downwash at horizontal tail with propeller 
removed, deg 

Figure 4. 9. 3.1 

Variable 

A power 

Downwash increment at horizontal tail due 
to power, deg 

Figure 5. 1.1. 6 

Variable 

a — e -* E 

u P 

Inclination of slipstream centerline behind 
propeller relative to X-body axis 

Table 5, 1.3. 2 

Variable 


— 

Figure 5.1 

0.50 

z 

s 

Defined in Figure 5.2 

Table 5.1.3. 2 

Variable 

'h 

Vertical distance from X-*body axis to 
horizontal* tail, ta(ft) 

Figure 5.2 

-0.28 (-0.92) 


Distance along X-body axis from c./4' to 
Cj^/4, m(ft) 

Figure 5.2 

4.87 (15.99) 

^\ff 

Vertical distance from o, /4 to slipstream 
centerline, m(ft) 

Eq. (5.1.1.21) 

Variable 


Increment in dynamic-pressure ratio at 
horizontal tail 

Figure 5.1.1. 7 

Variable 

“h 

Angle of attack at horizontal tall 

Eq. (5.1.1.22) 

Variable 

^‘^h^%j^prop off 

Propeller-off dynamic pressure ratio at 
horizontal tail 

Figure 4. 9. 3. 2 

1.0 

^\(hf)^^h,% 

of tail referenced to Call area and a 
j. ^dynamic pressure ratio of 1,0 

Figure 4.13.3.1 

Variable 

wfa 

Tail-off with power on referenced to 

Table 5.1.3. 9 

Variable 
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Table 5. 1.3. 4: Concluded- 



















Table 5. 1.3. 5: Power effects on maximum lift. 


Symbol 

Description 

Reference 

Magnitude 

K 

Corr^atlon factor 

Figure 5. 1.2.1 

Variable* 

S./S 
i w 

Ratio o£ total immersed wing area to S 

at propellers-off (C^ ) 

L - max 
wfn 

Table 5a.3e2 

Variable 

L power 

Increment ±u lift due to power at propeller*" 

off (C^ ) 

L f max 
wfn 

Table 5. 1.3. 3 

Variable 


T ’ 
c 

w- \ w / 

Table 5,1*3. 2 

K 

Figure 5. 1.2.1 

(AC^)’ 

L power ■ 

Eq. (5. 1.2.1) 

laax 

Eq.(5.1.2;2) 

0 

0.1666 

1.065 

0.0068 

0.0072 

0.0915 

0.1981 

1.097 

0.0673 

0.0744 

0.1970 

0.2157 

1.110 

0.1432 

0.1590 
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1 


X 



A, = A- 

(Cr:) 


Vc,=i 


r<bi)el 

/I + 2 Aie\ 

— 1 
CNJ 

V * K ) 


(r^\ >e *e 


q=|(crj)^ 


1 + I 


le 


Si=bjCj 

Ai=-L 

q 


•Figure 5.1: Definition sketches for calculation of immersed areas. 
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Centerline of slipstreatn 






Figure 5.2: Definition sketch for calculation of propeller power 

effects 




p 'at 0.75 radius, deg 

Figure 5. 1.1.1: Propeller normal force- parameter (Reference 3) 
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■ 8Rp2 


Figure 5. 1,1, 2: Propeller inflow factor (Reference 3) 
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Figure 5. 1.1. 5: Correlation parameter for additional wing lift due 

to propeller power (Reference -3) 
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Single-engine airplanes 


8Rr 


Figure 5,1.1. 6: Increment in downwash due to propeller power 

(Reference 3) 
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^^^Wopsoff 



Sw(Tc/prop) 

8Rp2 


Figure 5. 1.1. 6: Concluded 

f 






5»X.X.7s Effect of piropelXet power on the dynamic pressure 
ratio at horizontaX taiX (Reference 3) 
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Figure 5.1.2,!: Correlation' factor for, maximum lift due to power 

(Reference 3) 



Figure 5. 1,2. 2: Construction of the power-on tail-off lift curve 



Figure 5. 1.3.1; Tail-off lift characteristics of the ATLIT airplane 
for different power settings. AC^ obtained from 

max 


Table 5.1. 3. 5 = 2.3 million) 

t'4 


2k9 







Figure 5. 1.3. 2: Comparison of predicted lift curve with full-scale 

wind tunnel data (T ’ = 0,0915, no stabilizer 
deflection) ^ 
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5,2 Propeller Power Effects on Pitching Moment 


Power effects of propellers introduce increments of pitching 
moment due to direct action of the propeller forces offset from the 
center of gravity and propeller induced slipstream effects on the 
wing, nacelles (or fuselage), and the horizontal tail. Although all 
the increments of lift due to power (Section 5.1) contribute to power 
induced increments of pitching moment, several additional contributions 
must be considered. These contributions include the propeller slip- 
stream dynamic pressure effect on G and nacelle (or fuselage) free 

. o 

moments . 

The pitching moment of the complete airplane can be considered as 
follows ; 


* ^^m^prop off + + (AC_)^ + (iC_), + (AC_)*- + (AC_) 




m Aq 


"m'S ■ ' m'Aq. 
p 


’ni (Ae. ) 

h power 


(5.2;i) 


where 


^^m^prop off propeller-off pitching moment obtained from 

Section 4.13. 

(AC )„ is the pitching moment due to offset of the thrust from 
mi. 

the center of gravity. 

(AC^)jj is the pitching moment due to offset of propeller normal 
P 

force from the center of gravity.- 

(AC^ )^- is the effect of propeller slipstream dynamic pressure 

™o S? 

increment on zero-lift pitching, 

(AC ) .- is the effect on pitching moment due to change in wing 
m 

lift resulting from change in propeller- induced dynamic pressure. 
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(AC^)e^ is the pitching moment due to change in wing angle of attack 
resulting from propeller slipstream. 


(AC^)g is the effect of propeller slipstream on body (nacelle or 
“ P 

fuselage) free moments . 


(AC^)^- is the pitching moment due to change in dynamic pressure 
acting on the horizontal tail, 

(AC ) /• A \ is the pitching moment' due to change in angle of 

h'’ power 

attack at the horizontal tail. 


In the following discussion of power effects on pitching moment 
the pitching moment of the complete airplane will be considered as fol- 
lows : 




P ““L "“p Vhf) 


(5. 2 ; 2) 


where 


(C ) -- is the propeller-off, tail-off pitching moment ob- 

^wfn IT o 

tained from Section 4.8, 

(AG^)^ is the net effect on pitching moment due to change in wing 
h 

lift resulting from propeller slipstream induced dynamic pressure and 

angle of attack changes on the wing, or: 

(AC ) ■= (AC ).- + (AC ) (5.2,3) 

m Wj^ ™ “ ^p 

C is the net pitching moment contribution of the horizontal 

™h(hf) 

tail for power-on conditions, which can be obtained as follows: 


C = (C ) + (AC ),_ 

"h(hf) "h(hf) 


(5.2.4) 
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The above mentioned increments in pitching moment can be determined 
by the following steps. 

The pitching moment increment, (AC )_ due to propeller thrust is 

m i 

obtained from: 

" n(T ’/prop) (5.2.5) 

mi c c 

w 

where 

n is the number of propellers. 

is the moment arm of the thrust relative to the center of grav- 
ity, obtained from Figure 5.2. 

c is the wing mean aerodynamic (geometric) chord, 
w 

T^’/prop is the thrust coefficient due to one propeller obtained 
from Section 5.1. 


The pitching moment increment, (AC^)^ , due to propeller normal 

“ P 

force can be determined as follows: 


(AC ) - (AC ) 

m N L N c cosa„ 

p p w T 


(5.2.6) 


where 


(ACj^)^ is the increment in lift due to the normal force of the 
P 

propeller, obtained from Section 5.1. 

Xp is the moment arm of the propeller force relative to the center 
of gravity, obtained from Figure 5.2. 

is the angle of attack of the thrust axis, obtained from Sec- 
tion 5.1. 
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The zero-lift pitching moment increment, (AC ) - , due to pro- 

m Aq_ 
o ^ 

peller slipstream effects 'on immersed portions of the wing-body or 
wing-nacelles at zero-lift condition can be obtained as follows: 

a) Calculate the zero-lift pitching moment coefficient of those com- 
ponents of the tail-off configuration that are not immersed by the 
propeller slipstream as follows: 



) 


area not 
immersed 


(C ) 

^o ^^prop off 


S -S. 
w 1 


w 


S -S. 
w 1 


(b -b.)c 
W 1 w 


(5.2.7) 


where 

(C ) is the C of the wing with the propellers removed, 

d W f- tm m 

o prop off o 

obtained from Section 4.5. 

b) Calculate the zero-lift pitching moment coefficient, (C ) . 

d i - — 

o prop off 

of those components of the tail-off configuration immersed by the pro- 
peller slipstream as follows: 


^^m ^wB ~ ^*"m ^area not 

o prop off o prop off o . 

Tnrmiat-eoH 


(5.2.8) 


where 

propeller-off C of the wing and nacelles, 
o prop off o 

(C^ ^wn ’ ^ multi-engine airplane and the propeller-off C 

o prop off % 

of the wing and fuselage, (C^ )^^^ , for a single-engine airplane, 

o prop off 

obtained from Section 4.6. 

c) Calculate the zero-lift pitching moment due to the change in 
slipstream dynamic pressure by: 




Ao S. c. 

— ^ — (C ) . 

- S - m X 

'^03 ° P^op off 


(5.2:9) 
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where 

is the increase in dynamic-pressure ratio of the immersed 
portion of the wing, obtained from Section 5.1. 


The pitching moment encrenent, (AC ) , due to change in the lift 


m w. 


of the wing resulting from power effects, is obtained as follows: 


(AC ) - - 

m w 

li 




W 


w 


(5.2.10) 


where 

is the distance from the aerodynamic center of the immersed 

wing area to the center of gravity, obtained from Figure 5.1. 

and (AC^ ) are obtained from Section 5.1. 
p 




The pitching moment increment, (AC ) , due to propeller slip- 

m D' 
p 

stream effects on the nacelle free moments (for multi-engine configu- 
rations) or on the fuselage free moments (for single-engine airplanes) , 
is accounted for by calculation similar to those in Section 4.8 which 
considered the free moments due to wing induced flows with the propellers 
removed. The following expression indicates the nacelle free moments 
increment : 


where 


(AC ) 
m n 




36.5 S c 
w w 



dx 


(5.2.11) 


aiid -e^ are propeller induced changes in flow inclination on 
the nacelle obtained from Section 5.1. 

4q^/q„ can be obtained from Section 5.1. 

y dx of the nacelle is obtained from Section 4.8- 

For single-engine airplanes the effect of power on the free moments 
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of the fuselage should be accounted for. The procedure is identical 
to the one above. 


The net pitching moment contribution of the horizontal tail for 

power-on conditions, C , is: 

■ “h(hf) 

^h 

= - — C (5.2.12) 

c h(hf) 
w 

where 

is the distance from the center of gravity to the quarter 
•chord of the mean aerodynamic chord of the horizontal tail, obtained 
from Figure 5.1, 

C, is the lift coefficient of the horizontal tail obtained 

h(hf) 

from Section 5.1. 

5.2.1 Pitching Moment Characteristics of the ATLIT Airplane 

Tables 5. 2. 1.1 to 5. 2. 1.4- account for' the pitching moment incre- 
ments due to the direct propeller forces and power induced slipstream 
effects on the wing and nacelles. These increments are summarized and 
added to the propeller-off, tail-off pitching moments in Table 5. 2.1.5 
to provide power-on, tail— off characteristics. These characteristics 
are added to the power-on horizontal tail contributions (in Table 5. 2. 1.5) 
to provide the pitching moment characteristics of the complete ATLIT- 
airplane. 

The results of Table 5. 2. 1,5 are plotted in Figures 5. 2. 1,1 and 

5. 2. 1.2 for T ' = 0.0915, while the results for T ’ = 0.1970 are 
c c 

shown in Figures 5. 2.1.3 and 5. 2.1.4. 

The- predictions have been performed for a Reynolds number of 2.3 


°h(hf) 
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million. No power-on wind tunnel .data, however, were available for 

the ATLIT in the "fully clean" configuration at a Re 3 molds number of 

2.3 million. Data were available for a Reynolds number of 3.5 million. 

In the case of the ATLIT, this increase in Reynolds ntimber will result 

in an increase in pitching moment, =: 0.03, in the linear lift 

'■^Re 

region. The predicted pitching moment, including the Reynolds number 
correction, is also shown in Figures 5, 2, 1.1 through 5. 2. 1.4. 

The predicted pitching moment curves show fair agreement with 
the power-on full-scale wind tunnel data. The discrepancy between the 
predicted and experimental curves is caused mainly by the pitching 
moment contribution of the engine cooling system. The prediction 
method does not take into account the pitching moment due to -the 
engine cooling, while from the wind tunnel data (s.ee Section 4.11) it 
follows that the contribution is significant. 

When the predicted increment in pitching moment due to power is 


added to the experimental pitching moment curve' obtained with propellers 
removed, cowl flaps open, and engine inlets open, good agreement is 
obtained with the pox-^er-on wind tunnel data of Reference 2, The incre- 
ment in pitching moment due to power -can be obtained as follows :■ 

•^^^m^power ^m ~ ^*"m^prop off (5. 2. 1.1) 

where 


^m predicted pitching moment coefficient of the airplane 

including power effects, obtained from Table 5. 2. 1.5. 

^*^m^prop off predicted pitching moment coefficient of the 

aitplane with propellers removed, obtained from Section 4,11. 
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Table 5. 2. 1.1: Pitching moment increment due to propeller forces 


Symbol 

Description 

Reference 

Magnitude 

n 

Number of propellers 

Figure 5.1 

2 


Wing moan aerodynamic chord) m (ft) 

Table 2.1.1 

1.225 (4.018) 


Distance from X-body axis to thrust line, m (ft) 

Figure 5.2 

-0.128 (-0,417) 

X 

p 

Distance from propeller plane to center of 
gravity, m (ft) 

Figure 5.2 

2.01 (6.59) 

“t 

Angle of attack of thrust axis 

Eq. (5, 1.1,1) 

a 

T^’/ptop 

Thrust coefficient per propeller 

- 

Variable 

P 

Normal force coefficient of the propellers 

Table 5. 1,3.1 

Variable 



(AC^) 

j, 5 Table 
P 

5. 1.3.1 


Eq. (5.2.5) 

P 

; Eq. (5,2 

.6) 

a, deg 

T ' 

c 

T ' 
c 

T * 
c 


0 

0.0915 

0,1970 

0 

0.0915 

0.1970 

0 

0.0915 

^0.1970 

-4 

-0.0029 

-0.0031 

-0.0033 

0 

-0.0095 

-0.0204 

-0.0047 

-0.0051 

-0,0054 

“2 

-0.0013 

-0.0014 

-0.0015 

0 

-0.0095 

-0.0204 

-0.0022 

-0.0023 

-0.0025 

0 

0.0004 

0.0004 

0.0004 

0 

-0.0095 

-0.0204 

0.0006 

0.0006 

0.0006 

2 

0.0020 

0.0022 

0.0024 

0 

-0.0095 

-0.0204 

0.0033 

0.0036 

0.0039 

4 

0.0036 

0.0039 

0.0042 

0 

-0.0095 

-0.0204 

0.0060 

0.0064 

0.0069 
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Table 5. 2. 1.2: Zero-lift pitching moment increment due to propeller 

power .• 


Symbol 

Description 

Reference 

Magnitude 

n 

Number of propellers 

Figure 5,1 

2 

"p 

Propeller radius, n (ft) 

Table 2.1 

0.97 (3.17) 


Reference wing area, (ft^) 

Table 2.1.1 

14.4 (155.0) 

b 

w 

Wing span, m (ft) 

Table 2,1.1 

12.19 (40.0) 


Wing mean aerodynamic chord, m (ft) 

Table 2.1.1 

1.225 (4.018) 


Total immersed wing area, (ft^) 

Table 5. 1.3. 2 

Variable 

n(bj,/prop) 

Total immersed span, m (ft) 

Table 5. 1.3, 2 

Variable 


Chord of immersed wing area, m (ft) 

Figure 5.1 

1.34 (4.38) 

o prop off 

Zero-lift pitching moment of wing, propellers off 

Table 4.5.1 

-0.0783 

m wn -- 

0 prop off 

Zero-lift pitching moment of wing plus nacelles, 
propellers off 

Table 4,6.1 

-0,0733 


Change in dynamic pressure ratio on immersed wing 

Table 5. 1.3. 2 

4.9098 (T^’/prop) 









Table 5. 2. 1.2: Concluded 




; Table 5, 

1.3.2 

b -b.j ft; Table 5 
w x’ 

.1.3.2 

(C ) Eq. (5.2.7) 

m area not 

° inmersed 

a, deg 

T * 
c 

T ' 
c 

X ' 
c 


0 

0.0915 

0.1970 

0 

0.0915 

0.1970 

0 


0,1970 


0.3565 

0.3563 

0.3560 

27.385 

27.391 

27.401 

-0.0457 

-0.0457 

-0.0457 


0.3529 

0,3525 

0.3523 

27,513 

27.521* 

2’?. 532 

-0.0460 

-0.0460 

-0.0460 

0 

0.3466 

0.3467 

0.3468 


27.729 

27.726 

-0.0465 

-0.0465 

-0.0465 


0.3379 

0,3388 

0.3395 

28.043 

28.009 

27.986 

-0.0472 

-0.0471 

-0.0471 


0,3264 

0.3287 

0.3301 

28.448 ' 

23.370 

28.318 

-0.0482 

-0.0480 

-0.0479 







(c ) ; Eq 

prop off 

(5.2.8) 

mm 

Bin 







T ' 



T ' 







c 



c 






0 

0.0915 

0.1976 

0 

0.0915 

0.1970 





-0.0326 

-0.0326 

-0.0326 

0 

-0.0028 

-0.0061 





-0.0323 

-0.0323 

-0.0323 

0 

-0,0028 

-0.0060 





-0.0318 

-0.0318 

-0.0318 

0 

-0,0027 

-0.0058 





-0.0311 

-0.0312 

-0.0312 

0 

-0.0026 

-0.0056 





-0.0301 

-0.0303 

-0.0304 

0 

-0.0024 

-0.0053 
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Table 5. 2. 1.3: Pitching moment increment due to power induced 

change in wing lift 


Symbol . 

Description 

Reference 

Magnitude 

, ■ — ^ 

X 

w 

Distance from c./4 to center of gravity, m (ft) 

Figure 5.2 

0.026 (0,086) 

% 

Uing mean aerodynamic chord, m (ft) 

Table, 2.1.1 

1.225 (4.018) 


Change- in wing lift due to power induced change in dynamic 
pressure on wing 

Table 5. 1.3. 2 

Variable 

P 

Change in wing lift due to power induced change in flow 
direction on wing 

Table 5. 1.3. 2 

Variable 




5. 1.3.2 

(4V 

; Table 
'P 

5. 1.3. 2 

(4C ) ; Eq. (5. 

m w^ 

2.10) 

a, deg 

T ' 
c 

T ’ 
c - 

T ’ 
c 


0 

0.0915 

0.1970 

0 

0.0915 

0.1970 

0 

0-0915 

0.1970 

-4 

0 

-0.0043' 

-0.0088 

0.0028 

0.0141 

0.0258 

-0.0001 

-0.0002 

-0.0004. 

-2 

0 

0,0098 

0.0203 

0.0012 

0.0061 

0.0111 

-0. 

-0.0003 

-0.0007 . 

0 

0 

0.0235 

0.0486 

-0.0004 

-0,0018 

-0.0033 

0. 

-0.0005 

-0.0010 

2 

0 

0.0365 

0.0756 

-0.0018 

-0,0093 

-0.01-71 

0. 

-0.0006 

-0.0013 

4 

0 

0.0485 

0.1007 

-0.0032 

-0.0164 

-0. 0301 

0.0001 

-0,0007 

-0.0015 




























Table 5. 2. 1.4: Pitching moiiient increment due to power effect on 

nacelle free moments 


Symbol 

Description 

Reference 

Magnitude 

n 

Mumber of nacelles 

Figure 5.1 

2 

S 

w 

2 2 

Reference wing area» m (ft ) 

Table 2.1.1 

14.4 (155.0) 

c 

w 

Wing mean aerodynamic chord, ra(ft) 

Table 2.1.1 

1.225 (4,013) 

-e 

u 

Upwash at propeller, deg 

Table 5. 1,3. 2 

Variable 

S 

Propeller induced downwash behind propeller, deg. 

Table 5. 1.3. 2 

Variable 


Change in dynamic pressure raClo on immersed wing 

Table 5. 1.3. 2 

4.9098 (T^'/prop) 

fw j' dx 

Integral of square of mean width of nacelle planfona 

Table 4. 8. 1,2 

1.15 (40.69) 

n 

segments of Ax length, m^(ft^) 


per nacelle 


o, deg 

c 

p 

,deg: Table 

— 
5. 1.3. 2 

Sy. deg 
Table 5. 1.3, 2 

(AC ) 
m 

^ ; Eq. (5 
P 

.2,11) 

T * 
c 

- 

T ^ 
c 


0 

0.0915 

0.1970 

0 

0.0915 

0.1970 

-4 

-0.1023 

-0.4217 

-0.6348 

0.0946 

o_ 

0,0014 

0,0029 

-2 

-0.0445 

-0.1835 

-0.2763 

-0.2155 

0.0009 

0.0017 

0.0026 

0 

0.0133 

0.0546 

0.0823 

-0.5255 

0.0018 

0.0021 

0.0024' 

n 

0,0710 

0.2928 

0,4408 

-0.8355 

0.0027 

0.0024 

0.0021 

Hi 

0.1233 

0.5310 

0.7993 

-1.1455 

0.0036 

0.0027 

0.0018 
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Table 5. 2. 1.5: Pitching moment characteristics with tail-off and 

tail-on with power on 


Table 5, 2.1.1 


Table 5. 2. 1.1 


).- i Table 5. 2. 1.2 


(iC ) ; Table 5. 2. 1.3 

ID W|_ 


-0.0095 -0.0204 
-0.0095 -0.0204 
-0.0095 -0.0204 
-0.0095 -0.0204 


-0.0095 -0.0204 


(AC^)^ ; Table 5. 2. 1.4 


-0.0051 -0.0054 

-0.0023 -0.0025 


0 


0 . 

0.0009 

0.0018 

0.0027 

0.0036 


0.0915 


0.0014 

0.0017 

0.0021 

0.0024 

0.0027 


.1970 


0.0029 

0.0026 

0.0024 

0.0021 

0.0018 


(C ) 
ra^j^'prop off 

Table 4. 8. 4.1 


-0.1256 

-0.0806 

-0.0365 

0.0067 

0.0490 


0 


-0.1304 

-0.0819 

-0.0341 

0.0127 

0.0507 


0.0006 

0.0039 

0.0069 


: i Eq. (5.2.2) 

"wfn 


T ' 
c 


0,0915 


-0,0028 -0.0061 
-0.0028 -0,0060 
-0.0027 -0,0058 

-0.0026 -0.0056 

-0.0024 -0,0053 




-0,0001 -0,0002 -0.0004 

0 -0,0003 -0.0007 

0 -0.0005 -0.0010 

0 -0.0006 -0.0013 


-0.0006 -0.0013 


0.0001 -0.0007 


0 . 

0.0455 


C j Table 5. 1.3,4 
HiChf) 


-0.0015 


J Eq. (5.2.12) 


0.1970 

0 

0,0915 

0.1970 

0 


0.1970 

-0.1550 

-0.0662 

-0.0682 

-0,0728 

0.2606 

0.2684 

0.2865 

-0.1076 

-0.0454 

-0,0504 

-0.0598 

0.1787 

0.1983 

0.2353 

-0.0607 

-0,0244 

-0.0346 

-0.0471 

0.0960 

0.1361 

0.1853 

-0.0146 

-0.0034 

-0.0165 

-0.0288 

0.0134 

0.0649 

0,1133 

0,0305 

0.0179 

0.0002 

-0.0151 

-0,0704 

-0.0008 

0.0594 


C^; Eq. (5.2.2) 


.1970 


0.1317 

0.1277 

0.1246 

0.0987 

0.0399 

































































265 



Figure 5. 2. 1.1: 


Comparison of predicted pitching moments with 
full-scale wind-tunnel data (T ' = 0.0915, no 
stabilizer deflection) ^ 
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Figure 5. 2. 1.4: Comparison of predicted pitching moments with 

experimental results (T^* = 0.1970, no stabilizer 
deflection) 



5.3 Propeller Power Effects on Drag 


The net drag change of the airplane due to propeller power results 

from: 

1) the component of the propeller thrust parallel to the X-stability 
axis 

2) the change in zero-lift drag due to slipstream dynamic pressure of 
those portions of the aircraft immersed in the propeller slipstream 

3) the change in induced drag due to the lift component of the direct 
propeller forces and the change in angle of attack of the immersed por- 
tions of the wing 

4) the change in cooling drag due to the power induced change in dyna- 
mic pressure acting on the immersed cooling system. 

For the subject airplane, where the propeller slipstream immerses 
the nacelle as well as a portion of the wing and horizontal tail, the 
drag with power on can be written as follows: 


Cr. " Cr, + ) + CAC^ + (AC„ ) + AC„ + (AC„) 

D • D _ ,, D T ' D w D h ' D 'n D. D cooling 

prop off o o <5 i ^ 

system 


(5.3.1) 


where 


^^D^prop off propeller-off drag of the complete airplane, 

obtained from Section 4.12 

is the component of the total thrust parallel to. the velocity 

vector, a positive thrust is equal to a negative drag contribution. 

(AC ) is the change in profile drag due to power effects on the 
u w 
0 


wing. 



(AC„ ), IS the change in profile drag due to power effects on the 
D h 
o 

horizontal tail. 

(AC^ ) is the change in zero-lift drag of the nacelles (including 
D n 
o 

nacelle-wing interference effects) . 

ACjj is the change in induced drag of the wing due to power effects, 
i 

(AC ) , . is the change in cooling system drag due to power. 

D cooling 

system 


The component of the total thrust parallel to the velocity vector, 
(AC^)„, can be calculated as follows: 

(ACj^)^ = -n(T^Vprop) cosa^ (5.3.2) 

where 

n indicates the number of propellers. 

T^’/prop is the thrust coefficient due to one propeller, obtained 
from Section 5.1. 

is the angle of attack of the thrust axis, obtained from Section 

5.1. 


The change in profile drag, (AC ) , due to power effects on the 

D W 
o 

wing can be obtained as follows: 


S . / prop Aq 

(“d \ T 

o o prop off w 


(5.3.3) 


where 

(C_ ) is the propeller-off zero-lift drag coefficient of 

o prop off 

the wing obtained from Section 4-. 12. 
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S^/prop is the immersed wing area per- propeller, obtained from 
Section 5,1, 

Aq^/q is the increment in dynamic pressure ratio, due to power, 
at the wing, obtained from Section 5.1. 


The change in profile drag, )j^» due to power effects on the 


horizontal tail is; 


o o prop off w “ 


(5.3,4) 


where 

is the propeller-off zero-lift drag coefficient of 

o prop off 

the horizontal tail, obtained from Section 5.1. 

S, /prop is the immersed horizontal tail area per propeller, ob- 
h. 

1 

tained from Figure 5.2. 

-Aq^^/q^ is the increment in dynamic pressure, due to power, at the 
horizontal tail, obtained from Section 5.1. 


The change in profile drag coefficient, (4C^ ) , due to power effects 

o 

on the nacelles with nacelle-wing interference effects included is: 


^n ^n(w) q 

o o prop off 


(5,3.5) 


where 


is the propeller-off zero-lift drag coefficient 
of the two nacelles with nacelle wing interference effects included, ob- 


\(w) 
o prop off 


tained from Section. 4. 12. 
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The induced drag increment, AG^ , of the wing due to propeller slip' 

i 

stream modification of the downwash over portions of the wing can be 
accounted for as follows: 

1) Calculate the lift coefficient, C^", as follows: 

+ (AC^)^/prop + (AC^)jj /prop (5.3.6) 
^^prop off ^ 


where 

(AG ) /prop is the lift component of the thrust per propeller, ob- 
L T 

tained from Section 5.1 on a per propeller basis. 

(AC ) /prop is the lift component of the propeller normal force 

per propeller, obtained from Section 5.1 on a per propeller basis. 

C is the lift coefficient of the wing alone with the 

L 

prop off 

propeller removed, h^ing stall angles extended to power— on stall 
angles (see Section 5.1.2), 

2) From Figure 5.3.1 obtain a value for the propeller drag factor, K^, 
as a function of the propeller correlation parameter, 

3) From Figure 5.3.2 obtain a value for the average propeller downwash, 

(e/e^), as a function of the propeller correlation parameter and the 

ratio, R /b . 

p w 

4) Calculate the effective propeller downwash angle, e, average over 
the entire wing by: 


3e 


where 


P P 


(5.3.7) 


8e /3a is obtained from Section 5.1. 

P P 

5) Calculate the' induced drag increment ratio on a per propeller basis 
by: 
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[(C ) /prop] /^t \ 

°l '' [ "oioo off j 

tt^A e 

1 j. w 


^ 1 
W j 

/' (iCj^)^Zprop + /prop\ 

''prop off \ " / 

180C, 

Li 

pi^op off J 

+ 

^ i 

V V / 




(5.3.8) 

where 

is the lift coefficient of the wing alone with propeller 
w 

prop off 

removed, having stall angles extended to power— on stall angles (see 
Section 5.1,2) . 

6) Calculate the induced drag increment, , of the wing due to pro- 

peller slipstream modification of the downwash over portions of the 
wing by; 


S. ^S.^w 

X X prop off 


[(C ) /prop] 

w power on 

(c~r ^ 

jj j/j 

i prop off 


- 1 


(5.3.9) 


where 


(Cj^ )^ is the propeller-off induced drag of, the wing, 

i prop off 

obtained from Section 4.12. 


The change in drag coefficient of the cooling system, (AC^,) , . , 

due to power induced changes ' in dynamic pressure behind the propeller 
acting on the cooling system immersed in the propeller slipstream is es- 


timated by; 


where 


cooling 

system 



cooling 

system 


) 


prop off 



(5.3.10) 


(C, 


D 


). 


cooling system 


■prop off 


is obtained from Section 4.12. 
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With Equation (5.3.1) the drag coefficient of the airplane inclu- 
ding power effects can be calculated. 

5.3.1 Drag Characteristics -of the ATLIT Airplane 

Calculations for power-on net-drag characteristics of the ATLIT 

airplane are summarized in Tables 5. 3. 1.1 to 5. 3. 1.4 as functions of 

the angle of attack, a, and the thrust coefficient, T^’. Table 5. 3. 1.1 

summarizes the zero— lift increments of drag due to power, Table 5. 3. 1.2 

the induced drag increments due to power, and Table 5. 3. 1.3 the change 

in cooling— system drag due to power. Table 5, 3. 1.4 s umm arizes all the 

power effects on drag and lists the power-on net drag. 

The results of Table 5.3.1. 4 are plotted in Figures 5. 3. 1.1 and 

5. 3. 1.2 and they show good agreement with the experimental results of 

Reference 2. The predicted results do not include a Reynolds number 

correction. The predictions have been performed for a Reynolds number 

of 2.3 million, while power-on wind-tunnel data for the ATLIT in the 

"fully clean" configuration were only available for a Reynolds number 

of 3.5 million. For the ATLIT an increase in Reynolds number of this 

magnitude will result in a slight decrease in drag, (ACj^)j^ = 0 to -0.0002. 

”Re 

Adding this correction to the predicted drag will result in an even 
better agreement with the wind tunnel data. 
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Table 5, 3. 1.1: .Zero-lift drag increments due to power 


Symbol 

Oescripfion 

n 

Number of propellers 

S 

w 

Reference wing area, (£t^) 


Change in dynamic pressure ratio at the wing 


due to power 

S^f prop 

2 2 

Immersed wing area per propeller, m (ft ) 

0 prop off 

Propeller-off zero-lift drag of wing 


Propeller— off zero— lift drag of horizontal 

o prop off 

tail referenced to tail area 


Change la dynamic pressure ratio at the 


'horizontal tall due to power 

S. /prop 

loimersed horizontal tail area per propeller. 

n. 

i 

m2 (ft2) 


Propeller-off zero-lift drag of both nacelles 

o ^ prop off 

with wing-nacelle Interference included, 


referenced to wing area 


Rafaratice 

Magnitude 

Figure 5.1 

2 

Table 2.1.1 

14.40 (155.0) 

Table 5. 1.3. 2 

'4.9098 (T^'/prop) 

Table 5. 1.3. 2 

h 

Variable 

Table 4.12.1,2 

0.00970 

Table 4.12.1.2 

. 0,00835 

■ Table 5. 1.3. 4 

Variable 

Table 5. 1.3. 4 

2.95 (9.68) 

Table 4.12.3.3 

0.01384 



S^/prop, ft^; Table 5. 1.3. 2 


(iC^ Eq. (5.3.3) 


Aqj^/q^; Table 5.1. 3.4 



27.626 27.614 27.591 
27.346 27.322 27.305 
26.863 26.873 26.380 
26.185 26.260 26.310 
25.298 25.471 25.583 


0.00003 

0.00004 

0.00004 

0.00005 

0.00006 



0.0915 

0.1970 

0,00078 

0.00167 

0.00077 

0.00165 

0.00076 

0.00163 

0.00074 

0.00159 

0.00072 

0.00155 



0.00011 

0.00013 

0.00013 

0.00014 

0.00017 


cv 

Eq. (5.3.5) 

T ’ 1 

c 1 

0 

0.0195 

0.1970 

0 

0.00311 

0.00669 

0 

0.00311 

0.00669 

0 

0.00311 

0.00669 

0 

0.00311 

0.00669 

0 

0.00311 

0.00669 


0.00392 

0.00392 

0.00391 

0.00390 

0.00389 


0.1970 




0.00847 

0.00847 

0,00845 

0.00842 

0.00841 
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Table 5. 3. 1.2: Induced drag increment due to power 


Symbol 

Description 

Reference 

Magnitude 

a 

Number of propellers 

Figure 5.1 

2 

b 

w 

Wing Span, m (ft) 

Table 2,1.1 

12.19 (40.0) 

“t 

Angle of attack of thrust axis, deg 

Figure 5.1 

a 

A 

w 

Wing aspect ratio 

Table 2.1.1 

io.32 

R 

P 

Propeller radius, m (ft) 

Table 2.1 

0.966 (3.17) 

^^^prop off 

Drag of airplane with propellers off 

- Table 4.12.8.1 

Variable 

0 

Zero-lift drag increment due to power 

Table 5. 3. 1.1 

Variable 

'^^prop off 

Induced wing drag with propellers off 

Table 4.12.4.2 

Variable 

prop off 

Lift coefficient of wing alone with 
propellers off and stall angles extended to 
power-on stall angles 

Table 4. 2.4.1 

Variable 

c ’• 

- 

Eq. (5.3.6) 

Variable 

3e /3a 
P P 

Race of change of propeller downwash with 
propeller angle of attack 

Table 5. 1.3. 2 

Variable 

i/Cp 

Averaged propeller downwash over wing span 
as a ratio of prop^ler downwash behind 
propeller 

Figure 5.3.2 

Variable 

G 

Effective propeller downwash averaged over 
wing, deg 

Eq. (5.3.7) 

Variable 


Propeller drag factor 

Figure 5.3.1 

Variable 



0 

0.0915 

0.1970 


0 

0.7056 

1.520 


R /b 
P w 

Figure 5.3.1 

3e /an 
P P 

Table 5.1. 3. 2 

0.0793 

4.0. 

0.0250 

0.0793 

3.76 

0.1033 

0.0793 

3.54 

0.1554 



e, deg 
Eq. (5.3.7) 


0 . 

0.01178a 

0,02611a 
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i prop off 


; Eq. (5,3,9) 


Table 4.12.4.2 


T * 
c 


0 

0,0915 

0.1970 

0,00037 

0.00005 

0.00095 

0.00239 

0.00075 

0.00002 

-0.00003 

-0.00002 

0,00329 

-0.00001 

-0.00001 

-0.00001 

0.00796 

-0.00006 

0.00034 

0.00087 

■0,01476 

-0.00013 

0.00099 ■' 

0.00255 









































Table 5. 3. 1.3: Change in cooling-system drag due to power 


Symbol 

De5crxpt3.on 

Reference 

Magnitude 

^ , ^prop off 
cooling 

system 

Cooling drag coefficient with propellers 
off 

Figure 4.12.7,1 

Variable 

s 

w 

Wing reference area, (ft^) 

Table 2,1.1 

14.40 (155.0) 

R 

P 

Propeller radius, n (ft) 

Table 2.1. 

0.966 (3.17) 


Change in dynamic pressure ratio at wing 
due to power 

Table 5. 1.3. 2 

4.9098 (T^'/prop) 


a, deg 

. . ^prop off 
cooling 

system 

Figure 4,12.7.1 

(AC ) - . ; Eq. 

D cooling 

system 

(5.3,10) 

T * 
c 

0 


.0.1970 

-4 

0.0083 

0 

0,00186' 

0.00401 

-2 

0.0067 

0 

0.00150 

0.00324 

0 

0.0055 

0 

0.00124 

0.00266 

2 

0.0045 

0 

0.00101 

0,00218 

4 

0.0040 

0 

0.00090 

0.00193 
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Table 5. 3. 1.4: Power-on drag of complete airplane , 


a, deg 

(ACjj)^; Eq. (5,3.2) 


; Table 5.3* 

■> 

1.1 


Table 5.3.1 

.2 


T ' 
c 

T * 

c 

T * 
c 


0 

0.0915 

0.1970 

0 

0.0915 

0.1970 

0 

0.0915 

0.1970 

-A 

0 

“0.09120 

-0.19652 

0 

0.00392 

0.00847 

0.00005 

0.00095 

0.00239 

-2 

0 

-0,091<!i4 

-0.19688 

0 

0.00392 

0.00847 

0.00002 

-0.00003 

-0.00002 

0 

0 

-0.09150 

-0.19700 

0 

0.00391 

0.00845 

-0.00001 

-0.00001 

-0.00001 

2 

0 

-0.09144 

-0.19688 

0 

0.00390 

0.00842 

-0.00006 

0.0003A 

0.00087 

4 

0 

-0.09128 

-0.19652 

0 

0.00389 

0.00841 

-0.00013 

0.00099 

0.00255 


cooling* 5, 3. 1,3 

ays tell) 

- power 


'^‘1* (5.3.1) 

D power on 

T ’ 
c 

T ’ 
C 

Table 4.12.8.1 

T * 
c 

0 

0.0915 

0.1970 

0 

0.0915 

0.1970 


0 

0.0915 

0.1970 

0 

0.00186 

0.00401 

0.00005 

-0.08455 

'-0.18165 

0.03724 

0.03729 

-0.04731 

-0.14441 

0 

0.00150 

0,00324 

0.00002 

-0.08605 

-0.18519 

0.03577 

0.03579 

-0.05028 

-0.14942 

0 

0.00124 

0.00266 

-0.00001 

-0.08636 

-0a8590 

0.03723 

0.03722 

-0.04913 

-0.14867 

0 

0.00101 

0.00218 

-o ; 00006 

-0.08619 

-0.18541 

0.04164 

0.04158 

-0.04455 

-0.14377 

0 

0.00090 

0.00193 

-0.00013 

-0.08550 

-0.18363 

0.04906 

0.04893 

-0.03644 

-0.13457 
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Figure 5.3.2: Average propeller downwash (Reference 3) 
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Figure 5. 3. 1.1: Comparison of calculated and wind-tunnel-determined 

drag at different power conditions - 
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Figure 5. 3. 1.2; Comparison of predicted drag with full-scale wind 
tunnel data at different .power conditions 
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CHAPTER 6 


CONCLUSIONS AND RECOMMENDATIONS 


In this report an analytical nethod is presented for predicting 
lift, pitching moment and drag of light, twin-engine, propeller -driven 
airplanes. The method is applied to the Advanced Technology rLight Twin- 
Engine (ATLIT) airplane and the calculated characteristics are compared 
with the full-scale wind tunnel data. 

The following conclusions can be made: 

1. The calculated lift curves show fair agreement with the wind tunnel 
results. However, the accuracy can be improved -as is indicated by 
the good agreement of the results obtained with the lifting sur- 
face programs of Reference 4 and 5 with the experimental data. 

2. The calculated pitching moment coefficients agree well with the 
wind tunnel results. 

3. In Section 4,12 the drag of the ATLIT airplane is predicted. The 
predicted drag shows good agreement with the wind tunnel data of 
Reference 2. 

4. The calculated lift and pitching moment for the airplane with 
deflected stabilizer show poor agreement with the experimental re- 
sults. 

5. In Chapter 5 the prediction of power-on aerodynamic characteristics 
is presented. The predicted power effects agree well with the wind 
tunnel results . 

The following recommendations are made to improve the accuracy of 
the predicted results: 
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1. In Section 4.2 the factor, has been presented. This factor 

describes the error in the wing lift curve slope, , obtained 

a 

with the Polhamus formula when compared with the value for 

a ■ 

obtained with the lifting surface method of Reference 5. A similar 

factor, K , is suggested for the wing angle of attack for zero-lift. 
”o 

Use of such a factor will greatly increase the lift coefficient pre- 
diction. 

2. The lifting line theory overestimates the wing maximum lift coef- 
ficient, . A study is suggested which would result in a cor- 

max 

rection factor, , as function of wing configuration and wing 

planform. 

3. In Reference 3 and 4 the dynamic pressure ratio at the horizontal 

tail, is assumed to be equal to one, as long as the tail is 

not situated in the wake. It is suggested to use q^^/q^ =1.0 only 
in the case of a T-tail, while q-j^/q„ assumed equal to 0.85 

for a fuselage mounted, stabilizer (ATLIT) and 0.95 for a fin 
mounted stabilizer. This change will result in a better agreement 
of the calculated "stabilizer deflected" results with the wind 
tunnel data. 

4. The wetted areas of the fuselage and nacelle have to be determined 
exactly to obtain an accurate prediction of the zero-lift drag coef- 
ficient of the airplane. Instead of using Figure 4.12.2.1 to pre- 
dict the wetted area of a body, the method described in Appendix F 
is recommended. 

5. The computer program of Reference 5 does not take engine nacelles 
into consideration. Inclusion of the nacelle will improve the 
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results obtained with this program considerably. This will result 
in predicted characteristics which show good agreement with the 
wind tunnel data and with the data obtained with the program of 
Reference 6, However, the program uses much less computer time 
than the program of Reference 6. 

6. The engine cooling system appears to affect the drag and 'pitching 
moment of the airplane considerably. A study is suggested to make 
it possible to incorporate this effect into the prediction of the 
aerodynamic characteristics . 
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APPENDIX A 


CALCULATION OF THE SHIFT IN ANGLE OF ATTACK 
FOR ZERO-LIFT PER UNIT WING TWIST IN THE 
CASE OF PARABOLIC TWIST 


According to Reference 8, the shift in angle of attack for zero- 
lift per degree of wing twist can be written as follows: 


^ I ^ C3£}dn (A.l) 

0 g 

where 

p == 2y/b is the non-dimensional spanwise station 
through C^ can be obtained from Figure 4. 2. 3. 4 
f is the lift distribution function obtained from Figure 4. 2. 3. 5. 
For parabolic twist the following expression can be written: 

= ap^ + bp + d (A. 2) 

At the wing root (p = 0) the ratio e/0 = 0, while at the wing tip 
(p - 1) e/0 = 1. Substituting this information into Equation (A. 2) 
leads to the following result: 

— = ap^ + (l-a)p (A. 3) 

The chord, c, at wing station, p, can be calculated by: 


c = (c^ - c^p + c^ 


(A. 4) 


The standard mean chord for straight tapered wings is: 

1+1 

c = c — 
g r 2 


(A. 5) 


Substituting Equations (A. 3) through (A. 5) into Equation (A.l) results 

/ 

in the following expression: 
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'La 


= -/ { + (l-a)T, }[ C2 -*■ C3f ] d. 


or: 


Ao 1 ' /2C \ 

-g^= -/ [^^^^Ci{an3+ (l-a)n2[ + ^^-H:3£jjan2 + (l-a)Tij+C2^{an2v/l-n^+ (l-a)TVi^ 


For unswept wings ~ 0) the lift distribution function, f, 

elliptical ; 


f = -Vi-n^' 

ir 

In this case Equation (A. 7) can be written as: 


^»-/p^|g^C 3 ^{an 3 +a^a)T,^) + J^{aT, 2 + a-a)n}+ (C2+ C3)^^ai,2Nii^+ (l-a)T^G,^}] 


From Equation (A. 9) follows: 

c 


.^£.. *1 (V S>7‘ 

^ -1 

+ i{iVi-n^+arcsimi)j + (C2+t3)^ (l-a)|--jS/(l-TI^)^'|J 
This expression -may be written as : 


. Ao 2C. 
0 


-Mil «-4 *t>+ (|- n>h «2+=3>Tl*<w-i> "•!! 


To determine the parameter, a, a third point of the twist 
bution line has to be known. In the case of the ATLIT: 

D f 0.395 e/9 = 0.1667 

From Equation (A. 3) follows: 

a = 0.9554 


(A. 6) 


|2 ^dn 

(A. 7) 


(A. 8) 


dn 

(A. 9) 


(A. 10) 


(A. 11) 
distri- 

(A. 12) 


(A.13) 
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The following data are known for the ATLIT: 


2irA 

c'^ cos A 
a 


c/4 


9.846 


X = 0.5 


(A. 14) 


From Figure 4. 2. 3. 4 follows: 

= 0.462 

= 0.179 (A. 15) 

C 3 = 0.359 

Using the information from Equations (A. 13), (A. 14) and (A. 15) in . 
Equation (A. 11) leads to the following result: 


Aa 

-^=-0.2705 ■■ . (A. 16) 

Because the sweep angle of the quarter-chord of the ATLIT is very 
small (A^y^<2“), the wing of the ATLIT may be assumed to be unswept. 
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APPENDIX B 


CALCULATION OF MAXIMUM LIFT COEFFICIENT FOR 
STRAIGHT TWISTED WINGS 

The method described in Section 4.2.3 is discussed in more 
detail in Reference 8. According to Section 4.2.3, the additional 
lift coefficient for C^ = 1 is: 


S 4,^. S* 


= r + -r — + 

Z 1 c/c„ IT 
a g 


c/c 


(B.l) 


and the basic lift coefficient is: 


^ ^Z ^ ^Z ^^4 
baa 


/ 4a \ 


cosA 


8 


(B.2) 


The ratio of the wing chord at spanwise wing station, p, to the mean 
geometric chord is , according to Appendix .A: 


c 

c 


X+1 


(Xn-h + 1) 


(B.3) 


In the case of linear twist the twist ratio can be written as: 


e 


(B.4) 


However, the ATLIT airplane has a parabolic twist distribution, as 
discussed in Appendix A. For the ATLIT the following twist ratio 
is valid : 


■|- = 0.9554ri^ + O'. 0446n 


The following data are known for the ATLIT wing: 
.2irA 


Cv 


= 9.846 


c, = 0.115 deg 


-1 


(Table 4.1.2) 


a 


(B.5) 
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= 1.59 (Table 4.1.2) 

max 

0 = -3 deg (Table 2.1) 

Aa 

— ~ = -0.2705 (Appendix A) 

a 

Then, from Figure 4. 2. 3. 4 follows; 

= 0.462 
= 0.179 
= 0.359 
= 0.560 


(B.6) 


The quarter chord sweep angle, ATLIT wing is nearly 

zero. In that case, the lift distribution function is; 


f = 




(B.7) 


and Equation (B.l) can be written as: 



- =1 + 


c/=g 



(B.8) 


Table B.l summarizes the calculations which lead to the value 
of the maximum wing lift coefficient. From this table follows that 
the minimum value of the ratio of (c - c ) to c at G = 1 is: 

Jj 

max b a 


- c„ 


max 


= 1.494 = (C^ ) 


max 


w 


(B.9) 


The minimum value of this ratio is considered to be the maximum lift 
coefficient of the wing. 
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Table B,l: Maximum wing lift coefficient prediction 


n 

c/c 

g 

Eq. (B.3) 

a 

Eq. (B.8) 

e/0 

Eq. (B.5) 

% 

Eq. (B.2) 

max b 
a 

0 

' 

1,333 

0.9758 

0 

0.0510 

1.577 

0.1 

1.267 

0.9999 

0.0140 

0.0496 

1.541 

0.2 

1.200 

1.0213 

0.0471 

0.0441 

1.514 

0.3 

1.133 

1.0386 

0.0994 

0.0343 

1.498 

0.35 

1.1 

1.0453 

0.1326 

0.0278 

1.494 

0.4 

1.067 

1.0506 

0.1707 

0.0203 

1.494 

0.45 

1.033 

0.0540 

0.2135 

0.0116 

1t498 

0.5’ 

1.0 

1.0552 

0.2612 

0.0019 

1.505 

0.6 

0.933 

1.0492 

0.3707 

-0.0203 

1.535 

0.7 

0.867 

1.0264 

0.4994 

-0.0454 

1.593 

0.8 

0.8 

0.9758 

0.6471 

-0.0710 

1.702 

0.9 

0.733 

0.8693 

0.8140 

-0.0913 

1.934 

1.0 

0.667 

0.4620 

1.0 

-0.0651 

3.582 
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APPENDIX C 


COMPARISON OF EXPERIMENTAL GA(W)-1 SECTION 
AIRFOIL CHARACTERISTICS WITH THIN AIRFOIL 
THEORY AND THICK AIRFOIL THEORY PREDICTIONS . 


In this study Reference 4 and Reference 5 have been used to predict 
the longitudinal aerodynamic characteristics of the ATLIT airplane. 
Reference 4 is based on the thin wing theory, while Reference 5 is 
based on the thick wing theory of Douglas Neumann. In this appendix 
a comparison will be shown between experimental results (Reference 9) 
and results obtained with the thin airfoil theory and the thick airfoil 
theory, respectively. 

In the thin airfoil theory the airfoil is represented by its mean 
camber line, z^. According to the thin airfoil theory, the lift coef- 
ficient is: 


Cji = 2r(a-a^)- 


(C.l) 


where a is the angle of attack with respect to the .airfoil chord in 


radians and is the angle of zero lift in radians: 


a = — / , ^ (1-cos 0) d0 
0 ir / dx ' ' 


(C.2) 


where : 


X = (1 - cos 0) 


(C.3) 


In the thin airfoil theory, the thickness has been found to have 

no effect on c and c- . In reality, the thickness does affect slightly 
x» HI 

both c.^ and c^. An improved thin airfoil theory is Weber’s method (Ref- 
erences 11 and 12) . This method shows that for thin airfoils and for 
small angles of attack the lift coefficient can be represented as follows 
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(C-4) 


c„ = 2ir(l + 0.8-)(a-a ) 

S, C O 

The thick airfoil theory takes the complete airfoil into account. 
The first method is the Douglas Neumann method (References 13 and 14) , 
which uses a surface source distribution. The second method replaces 
the airfoil surface by a vortex sheet instead of a source distribution 
(References 15 and 16). The resxilts of both methods, however, have been 
proven to be similar. 

In Table C.l the airfoil coordinates of the GACW)-1 airfoil are 
presented. The mean camber line of this airfoil can be represented 
by the following expression: 

10 

z = 2 3^ (C.5) 

c n 

n=0 


where : 

ao = 0. 

a^ - 0.22690069 
= -1.09038423 
a^ = -0.94538563 
a^ = 35.87784264 
a, = -176.70961777 

^ (C.6) 

a, - 442.91119511 

b 

a^ = -647.63095699 
ag = 558.01744471 
ag = -263.04762768 
a^Q= 52.38631204 

To calculate the angle of zero-lift, the slope of the mean camber line 
has to be determined: 
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(C.7) 


dz 


10 


c - (n-1) 

-z — =2 n a X- 
dx n 

n=l 

Substitution of Equation (C.3) into Equation (C.7) results in: 

9 


dz 

= 2. b cos^ 6 
dx n 

n=0 


(G.8) 


where : 


bg = 0.00734851 
b. = -0.00172714 


= 0.01065684 

b^ = 0.44789152 

b, = -0.27322050 
4 

b^ = -1.48035280 

b. = 0.30678088 
0 

b^ = 2.27192933 
bg = -0.03923674 
bg = -1.02317016 

Substitution of Equations (C.8) and (C.9) into Equation (C.2) leads to 
the following result: 


= -0.07856 rad = -4.5012 deg (C.IO) 

With Equations (C.l) and (C.IO) the lift coefficient of the GA(W)-1 
airfoil can be predicted and the result is plotted in Figure C.l. 

In Figure C.l the lift coefficient obtained with the thick airfoil 
theory method of Reference 17 is also shown. In the linear lift region 
the lift coefficient predictions of the thin airfoil theory and the thick 
airfoil theory are almost identical, and they show good agreement with 
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the experimental lift curve. The lift curve slope of the "thick airfoil" 
is steeper than the slope of the "thin airfoil" lift curve. This differ- 
ence is caused by the thickness effect as shown by Weber’s method. 
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Table C.l: NASA GA(W)-1 airfoil coordinates 


'camber line 


.02035 

-.03069 

.04165 


.05600 

.06561 

.07309 


.08413 


.10409 

.10500 

.10456 


.09917 

.09674 

.09374 


.08604 


.07639 


.06517 

.05913 


.04644 


.03313 

.02639 

.01965 

.01287 


-.00074 


-.00974 

-.01444 

-.02052 

-.02691 

-.03191 

-.03569 

-.04209 

-.04700 

-.05,037 

-.05426 

-.05700 

-.05926 

-.06265 


-.06517 

-.06483 

-.06344 

-.06091 

-.05683 

-.05396 

-.05061 

-.04678 

-.04265 

-.03830 

-.03383 

-.02930 

-.02461 

-.02030 

-.01587 

-.01191 

-.00852 

-.00565 


-.00248 

-.00257 

-.00396 


.001630 

.002955 

.005085 

.007370 

.008915 

.010155 

.011760 

.013045 

.014110 

.014935 

.015740 

.016415 

.017565 

.018605 

.019460 

.020085 

.020560 

.020890 

.021170 

.021390 

.021565 

.021675 


.021570 

.021280 

.020830 

.020280 

.019415 

.018520 

.017265 


.013740 

.011435 

.008585 

.005150 

.001040 

-.004285 
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APPENDIX D 


DETERMINATION 07 LIFT AND PITCHING MOMENT 
COEFFICIENT FOR THE ATLIT WITHOUT HORIZONTAL 
TAIL IN THE "FULLY CLEAN" CONFIGURATION 

During the ATLIT drag clean-up investigation several factors 

were established which affected the lift curve slope, C^ , and the 

a 

induced drag, Cjj , of the ATLIT. The two most important factors were; 
i 

1. leakage through the' spoiler-Fowler flap region along the 
wing trailing edge 

2. premature flow separation on the wing upper surface between 
the fuselage and the nacelles . 

The highest experimentally determined lift-curve slope was obtained 
with trailing-edge leakage eliminated and fillets installed at the 
wing-fuselage juncture. Therefore, in this study, any comparison with 
theory will be based on experimental data which includes these two 
fixes. However, these fixes were present on the ATLIT airplane only 
during the early phase of wind-tunnel testing. No data were obtained 
with horizontal tail removed until near the end of the test program. 

The effect of the horizontal tail on the airplane lift curve 
was established for the "power-off" condition for the airplane- "as built" 
at a Reynolds number of 2.3 i^llion. The curves are shown in Figure D.l 
and the incremental tail lift is established from these data. -Next 
the tail lift can be subtracted from the lift curve of the ATLIT air- 
plane in the "fully clean" configuration with the horizontal tail on. 

This procedure is shown in Figure D.2. 
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The same method can he used to determine the pitching moment 
coefficient of the airplane without horizontal tail in the "fully 
clean" configuration. In Figure D.3 the pitching moment increment 
due to the horizontal tail is determined. Figure D.4 shows the 
pitching moment curves of the complete airplane "fully clean" and 
also the pitching moment curve of the airplane without horizontal 
tail in the "fully clean" configuration. 
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Figure D.l: Determination of horizontal tail increment to lift 

(Airplane "as built," propellers stopped, Nj^^=2,3xl0 ) 
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Figure D.2: Determination of horizontal tail off lift 

(Airplane "fully clean," =2.3x10^) 

Re ' 


curve 
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Figure D.3: 


Determination of horizontal tail increment to the pitching 

moment (Airplane. , "as built," propellers stopped, N =2.3x10^) 

Re ' ' 
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Figure D.4: Determination of pitching moment coefficient with horizontal 
tail off (Airplane "fully clean," Nj^^=2.3xl0 ) 
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APPENDIX E 


DERIVATION OF THE DQWNWASH AT THE HORIZONTAL 
TAIL FROM THE FULL-SCALE WIND TUNNEL DATA 

In Reference 2 no experimental downwash results at tlie horizontal 

tail are shown. However, the average downwash at the horizontal tail, 

E, > as a function of the angle of attack, a, can be derived from the 
h 

graphs shown in Reference 2. 

In Appendix D the horizontal tail lift including tail-fuselage 
interactions, wing downwash and dynamic-pressure effects, has been de- 
rived from the experimental data as a function of the angle of attack. 
The lift of the horizontal tail, based on the reference wing area, can 
be written as follows: 

In Section 3.10 the lift of the horizontal tail has been calcula- 
ted (i^ =• 0) and the result was : 

0.0177C a - s. ) (E.2) 

hi(hf) ^ 

based on the reference wing area, S , a dynamic pressure ratio, q, / q , 

w n “ 

equal to unity arid valid up to = 12.2 degrees. 

With Figure D.l and Equation (E.2) the average downwash at the hor 
izontal tail, s^., can be calculated. The results of the calculations 
are summarized in Table E.l. 
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^ 

Table E.l: Average downwash over the horizontal tail of the ATLIT 
airplane 2,3 -million, no flap deflection) 


a, deg. 

%(hf) 
Figure D.l 

ct - deg. 

. 

Eq. (E.2) 

deg. 

-4 

-0,098 

-5.54 ■ 

1.54 

-2 

-0.075 

-4.24 

2.24 

0 

-0.052 

-2.94 

2.94 

-2 

-0.028 

-1.58 

3.58 

4 

-0.005 

‘-0,28 

4.28 

6 

0.018 

• 

1.02 

4.98 

8 

0.034 

1.92 

6.08 

10 

0.051 

2.88 

7.12 

12 ■ 

0.067 

3.79 

8,21 

14 

0.084 

• 4.75 

9.25 

16 

0.100 

■ 5.65 

10.35 

18 

0.120 

6.78 

11.22 

20 

0.140 

7.91 

12.09 
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APPENDIX P 


CALCULATION OF THE WETTED AEEA 
FOR FUSELAGE AND NACELLE 


In this appendix the wetted area o£ the fuselage and nacelle will 
be estimated in a manner different from the method described in Section 
4.12. The accuracy of the fuselage and nacelle wetted area estimation 
must be high to obtain good agreement between the predicted value of the 
airplane zero-lift drag coefficient and the experimental result. 

The wetted area of a body (fuselage or nacelle) will be calculated 
as follows; 

1. Determine the circumference, C , of the cross section of the body 
at a distance, x, from the nose. Repeat this n times, as is shown in 
Figure F.l. 

2. Plot the circumference, C^, as a function of the distance, x, from 
the nose, as is demonstrated in Figure F.2. 

3. Integrate the circumference, C , over the distance, x, from the nose 
to obtain the uncorrected wetted area of the body. This integration can 
be performed with Simpson's approximation as follows: 
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where C represents the circumference of the cross section at the center 
i 

of the interval Ax , . 

X 

4. Determine the area of overlap between body and lifting surface(s). 


The area of a wing section can be estimated as follows (Reference 8) : 


A = 0.68(t.c) = 0.68(t/c)c^ (F.2) 
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where c is the exposed root chord of the lifting surface, t/c is the 
thickness to chord ratio of the exposed root chord. 

5. Subtract the total area of overlap between body and lifting surfaces 
from the uncorrected wetted area to obtain the actual wetted area of the 
body. ^ 

The fuselage was broken into eight segments. The cross, ..sectional 
circumference was determined at nine stations: four nose sections, one 

windshield section, two cabin sections and two tail cone sections. The 
corrected wetted area of the fuselage appears to be; 

(S_ = 23.32m^ = 251. Oft^ (F.3) 

wet f 

The wetted area (including spinner) of one nacelle is : 

(S J = 5.55m^ = 59.7ft^ (F.4) 

wet n 

while the 8 inch extended nacelle has a wetted area: 

(S ) = 6.08ra^ = 65.4ft^ (F.5) 

wet n 
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Figure F.l: Stations at which cross-sectional circumference is 
determined 



Figure F,2: Circumference of fuselage crofis-section at different 
longitudinal stations 
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